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SUMMARY 


Recently, wing-body interference effects have assumed greater 
significance in aerodynamics than heretofore because of the trend 
toward lower aspect ratios and larger ratios of body diameter 
to wing span in high-speed aircraft and missiles. It has thus 
become necessary to effect refinements in the time-honored aero- 
dynamic assumptions concerning wing-body interference for both 
the subsonic and supersonic speed ranges. In this paper, a criti- 
cal survey of available theories, simplified analysis methods for 
use in design, and experimental results is given, including some 
new developments based on recent work at Cornell Aeronautical 
(1) an approximate subsonic 


Laboratory. Among these are: 


theory for low-aspect-ratio (4R < 4) wing-body combinations 
applicable to the class of wings with straight trailing edges, (2) 
a simple semiempirical design method for predicting subsonic 
characteristics of low-aspect-ratio wing-body combinations, (3) 
a calculation by variational methods of the characteristics of a 
series of high-aspect-ratio wing-body combinations according to 
lifting-line theory, (4) supersonic experimental results and anal- 
ysis for a series of wings with supersonic leading edges of aspect 
ratios ranging from 0.5 to 5.0, with particular emphasis on body- 
load distribution, (5) approximate methods for rapid calculation 
of body loading in supersonic flow for the configurations of (4 
above, and (6) general integral relations in the linearized theory 
of wing-body interference. It is also shown how a unified point 
of view in constructing approximate theories and design methods, 
and in interpreting and utilizing experimental data for all speed 
ranges can be achieved by the application of the general integral 
relations 


INTRODUCTION 


'y RECENT YEARS, the subject of wing-body inter- 
ference has assumed increased importance in aero- 
dynamics because of the design trend in transonic and 
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supersonic aircraft and missiles toward lower aspect 
ratios and lower ratios of wing span to body diameter. 
Many of the time-honored assumptions about wing 
body interference which served the aerodynamicist 
well in analyses of low-speed aircraft have been found 
to be woefully lacking when applied to some of the high 
speed aircraft now under development. For instance, 
the oft-used assumption that the slope of the lift curve 
of a wing-body combination can be closely estimated 
by calculating the slope of the lift curve for the wing 
alone, including the hypothetical extension of the wing 
which passes through the body is likely to give within 
5 per cent of the correct slope of the lift curve if the 
ratio of span to body diameter is 10 or greater, but it 
may lead to errors of the order of 20 per cent if this 
ratio is 3. The need for increased accuracy in the esti- 
mation of wing-body interference effects is not a pleas 
ant idea for the engineer to contemplate, since it sug 
gests the necessity for more complicated, laborious, and 
time-consuming calculations as part of any design anal 
ysis. It that the may 


not be as bad as it seems at first sight, since often 


appears, however, situation 
crude approximations and assumptions based on a 
sound physical understanding of the problem can lead 
to results satisfactory for many purposes. For ex 
ample, in estimating the total lift of a wing-body com 
bination, if interference effects account for 30 per cent 
of the lift, an error of 20 per cent in estimating the 
interference effects will lead to an error of 6 per cent 
in the answer, which is tolerable in many applications. 
There are few solutions of wing-body problems which 
are exact, even within the framework of the linearized 
wing theory. Most theories of wing-body interference 
are based upon approximations which involve assump- 


289 








290 JOURNAL OF THE 
tions in addition to those which are made in the various 
theories of wings alone. In many practical problems, 
it is found necessary to use these approximate theories 
for cases other than those for which they were derived, 
because of a lack of any theory more applicable to the 
problem at hand. To do this with any degree of suc- 
cess requires a knowledge of the basic assumptions and 
limitations of the available approximate theories and 
their relations to each other, as well as an understanding 
of the essential physical aspects of wing-body inter- 
ference. In the present paper, the available theories 
as well as certain recent extensions and new develop- 
ments are derived from a common basic formulation 
insofar as possible, so that the significant assumptions 
made in developing each theory and the relationships 
between them can be more clearly seen. A _ unified 
point of view is sought in the consideration of the vari- 
ous theories and approximate methods for both the 
subsonic and supersonic speed ranges. Relationships 
in the Trefftz plane and the reverse-flow theorems?! 
for wing-body combinations are shown to provide use- 
ful common bonds between almost all theories in both 
the subsonic and supersonic speed ranges. Because of 
limits on the length of this paper, only symmetrical lift 
problems are considered. While many of the theories 
and methods discussed apply to antisymmetrical lift 
problems as well as with minor changes, mention of such 
applications will usually be omitted. 

For subsonic flow, the spanwise lift distribution over 
a wing-body combination of minimum induced drag 
was first treated by Lennertz! for a body consisting of an 
infinitely long circular cylinder. His solution was gen- 
eralized by Pepper? to the case of wings in combination 
with an infinitely long body of any cross-sectional shape 
or with a multiplicity of such bodies. Lennertz also 
considered the problem of the longitudinal distribution 
of lift over the cylindrical body and found a distribution 
which was symmetrical about an assumed wing lifting- 
line and which was highly localized in the region of the 
wings. The case of a high-aspect-ratio wing of given 
plan form and angle-of-attack distribution mounted on 
an infinitely long cylinder of any shape was treated by 
Multhopp,* whose solution is based on_ lifting-line 
theory. Multhopp also indicated approximate correc- 
tions based on slender-body theory to allow for the 
effects of a noncylindrical body of finite length. Mul- 
thopp’s treatment of the lift on the body omitted that 
part of the body lift which would occur on an infinite 
cylinder of the same cross-sectional shape as the body 
in the region of the wings. This leads to discontinui- 
ties in the body loadings in the vicinity of the wings and 
requires rather arbitrary assumptions as to the down- 
wash to which the body should be considered subjected. 
An improved method for calculation of body loading for 
finite bodies has been developed and is presented herein. 
This method is based on slender-body theory for the 
body only and allows the calculation of the loading con- 


tinuously along the body. For an infinite cylindrical 
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body, Lennertz’s results are closely approximated. A 
solution for the lift of high-aspect-ratio wing-body com- 
binations based on a variational principle in the Prandtl 
lifting-line theory has been given by Flax and Treanor,! 
Slopes of the lift curve for a large range of wing-body 
configuration variables have been obtained and some 
are included herein. 

For the case of wings of aspect ratio less than about 
6, the Prandtl lifting-line theory predicts slopes of the 
lift curve for the wing alone which are known to be too 
high by an appreciable amount, and the same effect 
may be expected for wing-body combinations of aspect 
ratio below 6. Also, for swept wing-body combinations, 
the Prandtl lifting-line approach cannot be used. For 
wings alone, these cases are handled by the Weissinger 
theory, but attempts to set up a theory of the Weissinger 
type for wing-body combinations meet with difficulties 
in the determination of the downwash produced at the 
3/4 chord position of the wing in the presence of the 
body. Possible methods of attacking this problem are 
discussed, including one proposed for circular bodies 
by Zlotnick and Robinson® involving the use of a num- 
ber of horseshoe vortices on the wing, and their images 
in the body. It appears that while this problem can 
probably be solved very approximately by any one of 
a number of methods, the choice of the most suitable 
method will depend on the results of extensive correla- 
tion of theory with experiment. 

For slender-winged bodies of revolution having wings 
with straight trailing edges and cylindrical bodies be- 
hind the wings, Spreiter’ and Ward* have given closed 
solutions of simple form. Ward's results differ from 
Spreiter’s if the body is noncylindrical in the wing 
region. Miles has pointed out that this discrepancy 
results from the fact that Spreiter neglected quadratic 
terms in Bernoulli's equation which effect an interac- 
tion between the displacement flow around the body 
and the lifting flow around the wings. The slender 
body theory may be expected to be valid for exposed- 
wing aspect ratios such that the Jones theory® is valid 
for the wing alone. Although the range of these theo- 
ries is thus limited, the relatively simple form of the re- 
sults make them very useful in setting up approxima- 
tions, and in further extensions of the theory. Results 
of this type for winged bodies of elliptical section have 
been given by Nonweiler" and Kahane,'! while the case 
of a body of revolution separated from the wings by 
a gap has been solved by Mirels. By setting up the 
Spreiter-Ward formulation of the wing-body problem 
as a variational problem, it is shown herein that the lift 
distributions obtained for uniform angle of attack of the 
wing and body must always be the same as for the 
classical minimum induced drag case, as has always 
been observed in all solutions obtained. Further, the 
total lift of a slender wing-body combination can be 
readily computed from the results of the minimum in- 
duced drag problem with little effort. Another result 
which speeds up calculations is that the total lift of the 
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WING-BODY 
slender wing-body combination is the same as that of 
the slender wing alone which is obtained from the 
conformal transformation which maps the body on a 
vertical slit (Multhopp’s transformation). By the use 
of these two results, the total lift as well as the chord- 
wise distribution of lift can be obtained for a wide va- 
riety of cases with little effort, once the transformation 
appropriate to the body shape is known. 

A theory for low-aspect-ratio wings with straight 
trailing edges centrally mounted on an infinite circular 
cylinder has been given by Lawrence;!'* this theory 
is valid for aspect ratios well above those for which 
slender-configuration theory becomes inaccurate. A 
generalization of this theory to winged cylindrical 
bodies of any cross-sectional shape is given in this paper. 
The results of this theory always lead to an equivalent 
low-aspect-ratio wing problem, which may then be 
solved by methods already set up for such problems. '* 

The lift on afterbodies in the presence of wing trailing 
vortices has been studied by Lagerstrom and Graham, '* 
Moskowitz,'® and Adams and Sears.'® The work of 
Lagerstrom and Graham, as well as that of Moskowitz, 
deals with the forces induced on a long cylindrical after- 
body at an angle of attack as the wing vortices pass 
downstream at an ever-increasing vertical distance from 
the body. Lagerstrom and Graham found that there 
was a small additional positive lift on the body, while it 
can be shown that Moskowitz’s method results in a 
negative lift on the body equal in magnitude to all the 
positive lift induced on the body by the wing in the 
first place. It appears that Moskowitz is in error 
because of his neglect of the pressures induced on the 
body by the deformed vortex sheet as it passes over the 
top of the body. Adams and Sears studied the case of 
an afterbody of decreasing area in the presence of a 
wing vortex sheet; a small additional positive lift re- 
sults. The analyses of Adams and Sears, and Lager- 
strom and Graham may be applied as a first approxi- 
mation for slender afterbodies behind nonslender wing- 
body combinations. 

It has been found that rapid estimates of total lift of 
a winged body which is cylindrical in the region of the 
wings with constant angle of attack of the wings and 
body can be made by a relatively simple approximate 
method, by considering the exposed wings as the basic 
lifting element. The method was apparently first 
suggested by Ward,* who based it on the results of 
slender-configuration theory for winged bodics of revo- 
lution, but was unable to verify its validity for other 
cases. In the present work, it is shown on theoretical 
grounds that the method is approximately valid for 
both high- and low-aspect-ratio wings in combination 
with bodies of revolution, and it is shown how the 
method may be extended to cylindrical bodies of arbi- 
trary shape. There is, however, a need for more ex- 
perimental data on this subject. 

For supersonic flow, the earliest work is apparently 
that of Kirkby and Robinson," who analyzed the case 
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of a very-high-aspect-ratio wing mounted on a conical 
body by approximate methods. They did not deal 
with any action of the wing on the body, and calcu 
lated the interference effect as simply the increased 
lift on wings due to the increase in wing angle of attack 
produced by the upflow around the body. The next 
work in this field was that of Ferrari," 
sidered a rectangular wing of high aspect ratio mounted 
His method involves 


who con 


on a pointed body of revolution. 
the use of an iterative process which starts with the cal 
culation of the flow about the wing alone and the flow 
about the body alone, and proceeds to correct the solu 
tions in such a way that at each stage of the iteration 
process, the boundary condition on either the wing or 
the body is satisfied exactly, while the boundary condi 
tion on the other component is violated. Since the 
process can be shown to converge, carrying the proce 
dure through a sufficient number of steps would even 
tually result in a situation in which the boundary con 
ditions on both the wing and the body are satisfied to 
any desired numerical accuracy. Actually, the pro 
cedure becomes very laborious after the first few steps, 
and numerical with the theory have 
never been made beyond the first two steps. While in 
principle the iterative method will give the correct solu- 
tion, no matter what initial assumptions are made to 
start the process, the number of steps required to 


calculations 


achieve a given accuracy will obviously depend to a 
large degree on the strarting point of the iteration. 
The starting point of the Ferrari method was chosen so 
as to produce the smallest errors possible in the early 
steps, subject to the restraining condition that the cal 
culation of the early steps be kept as simple as possible. 
The results of the first iteration for wing and body show 
fairly good agreement with experiment. 

Morikawa” has attacked the problem of a rectangu 
lar wing mounted on a cylindrical body when both are 
at the same angle of attack. The region considered in 
his analysis is limited to points forward of the wing 
trailing edge. By using the Laplace transform in the 
x direction, he has succeeded in finding the Laplace 
transform of the exact solution to this problem in terms 
of a series involving the modified Bessel functions. . He 
was, however, unable to invert this transform to ob- 
tain the solution. He has given an approximate solu- 
tion to the problem which should be valid for small 
distances from the wing leading edge. The accuracy of 
this solution for practical configurations has not been 
determined and must depend upon experimental veri- 
fication. Morikawa’ also has given the solutions to 
numerous planar problems intended to simulate the 
wing-body problem; these solutions may be instructive 
for qualitative analyses and for setting up approximate 
methods for use in design. This will be discussed later. 

Browne, Friedman, and Hodes,** using the powerful 
methods available for conical-flow problems, have 
solved the supersonic wing-body problem for the case 


of a delta wing which starts at the apex of a conical 
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body. Their solution is exact within the limits of lin- 


earized theory. Although the configuration they have 
considered is not a particularly practical one, the re- 
sults are useful in checking the validity of those ap- 
proximate theories which can also be applied to this 
configuration. 

Lagerstrom and Van Dyke*! have made a very com- 
prehensive study of the supersonic wing-body inter- 
ference problem. Their arguments and results were, 
however, entirely qualitative and for the most part 
based on extension of the results obtained for planar 
problems. Of special interest in their work is a formu- 
lation of the wing-body problem, which was later used 
by Nielsen** in his work. In this formulation, which 
will be stated here for the case of the wing only at a con- 
stant angle of attack, a wing at angle of attack extend- 
ing through the body is taken as the basic wing. This 
wing is assumed to have supersonic leading edges. The 
wing flow then induces normal velocities across the body 
surface. If a potential can be found satisfying the 
linearized differential equation for supersonic flow and 
cancelling these normal velocities at the body surface 
without producing any velocities normal to the wing 
surface, the exact solution is the sum of this potential 
and the basic wing potential. Lagerstrom and Van 
Dyke have shown that such a solution can be con- 
structed by expanding the potential in a Fourier series 
in the meridional angle and considering only even terms, 
each of which contributes zero velocity normal to the 
wing. This formulation is valid only up to the wing 
trailing edge. 

Nielsen** has solved the supersonic problem of a rec- 
tangular wing at constant angle of attack centrally 
mounted on a circular cylindrical body at zero angle to 
the stream, using the Lagerstrom-Van Dyke formula- 
tion of the problem. Nielsen employed the Laplace 
transform in the stream direction to obtain his solu- 
tion. He was able to invert the transform of the solu- 
tion only in terms of a new set of functions, one for each 
term of the Fourier series. Only the first four of these 
functions were tabulated, so that only four terms of the 
for numerical 


Fourier series solution could be used 


work. The numerical results showed that the Fourier 
series solution converged very slowly in the vicinity of 
the wing leading edge and the solution needed to be 
faired into a result obtained by another approximation 
forward of x = 12 VM? — 1R. 
Nielsen and Pitts’ have given an additional presenta- 
tion of Nielsen’s theory with more emphasis on physical 
interpretation and application than is to be found in the 


More recently, 


original paper. 

Consideration of the practical side of the wing-body 
interference problem discloses that, as in the subsonic 
case, in many instances rather crude approximations 
will suffice to give satisfactory answers for many pur- 
poses. In considering any questions of accuracy in a 
practical way, due account should be taken of the limi- 


tations of the linearized theory itself. Second-order 
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effects on pressures cannot be readily computed for 
complex configurations, but are known to be of con- 
siderable magnitude. Also, the disturbances from wing 
to body and back do not propagate exactly along the 
Mach lines as assumed in the theory, and this, too, pro- 
duces some errors in the theoretical results. Recourse 
to experiment must be taken to assess the magnitude 
of such effects in many instances. Cramer” has given a 
comprehensive presentation of experimental pressure 
distributions at a Mach Number of 2 for a class of wing- 
body combinations consisting of bodies of revolution 
with centrally mounted wings having no subsonic edges. 
The aspect-ratio range covered was from 0.5 to 4.29, 
based on the exposed span of the models, and the body 
radius to wing chord ratio was approximately 0.7, 
Based on these test results, Cramer was able to show 
that fairly simple “‘rules of thumb’ augmented by the 
reverse-flow and Trefftz-plane theorems for wings on 
cylindrical bodies would suffice to give good estimates 
of the lift on the configurations tested. The results 
made it clear that the build-up of lift on the body behind 
the trailing edge contributes the greater part of the body 
lift, so that attention should be focused on the after- 
body for most configurations. The pitching moment 
of the configuration is also very strongly dependent on 
the afterbody lift distribution, and thus on the after- 
body length. No theory or method for estimating 
wing-body interference effects which fails to account for 
the afterbody length can thus be of general utility. 
Cramer has shown fairly good agreement between the 
body-lift distribution computed from the Ferrari 
theory and those obtained from experiment. A nor- 
malized body loading curve based on the Ferrari theory 
which has been calculated by Rae” is included herein. 
This curve can be used directly for all cases in which 
the wing tip Mach cones do not intersect the body; 
approximate methods of handling the effects of the wing 
tips when they do influence the body are discussed in 
the present paper. The application of slender-body 
theory as indicated by Ferrari's and expanded by 
Schrenk” offers a useful means of obtaining approximate 
evaluations of such effects. 

The interference pressures on the wing should theo- 
retically be obtainable by either the method of Fer- 
ari!®: 19 or that of The Ferrari iteration 
process has, however, not been carried far enough to ob- 


Nielsen.” 


tain the wing pressures to a high degree of accuracy. 
The Nielsen method is set up in such a manner as to 
favor the determination of wing pressures at the ex- 
pense of the lift distribution on the afterbody. Thus, 
the accuracy of Nielsen’s wing pressures should be 
greater than the Ferrari wing pressures for a given 
amount of numerical effort. As indicated previously, 
however, Nielsen’s method is limited in accuracy for 
distances downstream of a wing leading edge—body 
intersection of the order of the body radius. For wings 
of large chord, the Nielson method should give satis- 


factory accuracy. Even in this case, the wing span 
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WING-BODY I 
loadings obtained from the first iteration in the Ferrari 
method do not differ greatly from those of Nielsen or 
from experimental results. In studying experimental 
results for practical configurations, the departures from 
the theory due to second-order effects on the pressures 
are found to be significant. Also, the effects of the gap 
opened at the wing-body juncture in the case of mov- 
able wings are marked. *° 

For the case of wing and body at the same angle of 
attack, Ferrari's iteration the ap- 
proximate solution of Morikawa”*’ provide the only avail- 


solution'» '* and 


able theoretical data on wing pressures, since Nielsen's 
Both 
the first step in the Ferrari method and the Morikawa 


method has not yet been extended to this case. 


method estimate a fall in wing lift near the body, al- 
though Ferrari's result is in considerably better agree- 
ment with experiment. An analysis of experimental 
results for rectangular wings discloses that the wing 
pressures in the region between the side of the body and 
the Mach line from the wing-body intersection average 
about 90 per cent of the value at the Mach line, at least 
for distances up to 0.7 V \/* — 1 body radii from the 
wing leading edge. The theoretical value at the Mach 
line from the leading edge can easily be predicted from 
linearized wing theory (the result was first given by 
Kirkby and Robinson"). At present, this approach 
appears to offer the simplest, and at the same time, most 
practical method for predicting wing pressures for the 
case of the wing and body at the same angle of attack 
for the range of wing chord to body diameter con- 
sidered. Various approximate methods for roughly 
estimating supersonic wing-body interference effects 
have been given by several authors, including Laitone,”™ 
Lagerstrom and Van Dyke,*! Morikawa,** and Beane.*° 
These methods are based on the use of strip theory and 
analogies to planar models at one limit and on slender- 
configuration theory at the other limit. The results 
for intermediate cases have been inferred by means of 
plausible Such 
methods must be used with care, but when they are 
skillfully handled, and augmented even by meager ex- 
perimental results, they can be quite useful in design 


assumptions and approximations. 


and in the planning and interpretation of experiments. 


GENERAL THEORY AND CONCEPTS 


From the discussion in the introduction, it can be 
noted that a variety of methods applicable to special 
regimes of flight speed and configuration, and employ- 
ing different types of approximations, are available. 
In an attempt to unify the results obtained by different 
investigators and to clarify the relationships between 
the various theories, a general formulation of the theory 
will be developed, and the paths leading to the special 
theories will be indicated. The general formulation 
will be applicable to both supersonic and subsonic flows. 
By proceeding in this manner, the basic similarities and 
differences between the subsonic and supersonic theories 
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will be brought to the fore. In spite of the complexity 
of the wing-body problem, there are some general con 
cepts which are applicable to wide classes of configu 
rations and within the limitations of linearized theory 
to both subsonic and supersonic flows. These con 
cepts will prove to be useful in the discussion of many 
particular cases and will therefore also be considered in 


this section. 


General Formulation of the Wing-Body Problem 
Consider a wing-body configuration immersed in an 
inviscid uniform free stream flowing with the velocity 
U parallel to the positive x-axis. The angle between 
the free stream and every surface element of the wing- 
body configuration is assumed to be small, so that 
linearized theory is applicable. The wing is assumed 
to be so thin that it may be replaced by its mean camber 
surface, which is assumed to lie in the neighborhood of 
The local angle of attack of the wing 


The curvature of the body 


the plane s = 0. 
will be denoted by a(v,y). 
will be assumed to be finite everywhere except perhaps 
at the nose. Aside from this restriction, and the re 
quirements imposed by the applicability of linearized 
theory, the shape of the body is arbitrary. Effects due 
to displacements of the wing trailing vortex sheet will 
be neglected in this formulation, since such effects are 
difficult to formulate generally. Some special cases 
will be considered later. 

The (x,y,z) components of velocity may be expressed 
in terms of the perturbation velocity potential ¢ in the 
form (u,v,w) = (Uo + (O¢/0x), O¢/Ov, O¢/Oz). Neces- 
sary and sufficient conditions for specifying a unique 
potential function may be written in the form 


0?y/Oy"? + (07¢/02") + (1 — M*) (0%¢/0x") = 0 
in the flow field F (1) 
0¢/d0z = —U, acos (x,y) 
on the wing surface Ry (2 
O¢g/On = —U 5 cos (n,x) 
on the body surface Rg (3) 
g—>0 asx -> —a (4 
0¢/dx = 0 at wing trailing edges Ry, (9 


Eq. (1) is the well-known linearized differential equa- 
tion for the velocity potential in a flow at a free-stream 
Mach Number of 1/7. 
conditions that the flow through any element of the 
The symbol 0¢ On de- 
in the direction of the 


Eqs. (2) and (3) express the 
wing or body surface vanishes. 

notes the partial derivative of ¢ 
outward normal to the body surface. The symbol 
(n,x) denotes the cosine of the angle between the 
outward normal to the body surface and the positive 


Eq. (4) insures that the perturbation velocities 


cos 


X-axis. 
vanish far upstream and the Kutta condition at the 
wing trailing edge is formulized in Eq. (5). 

The aerodynamic forces acting on the surface ele- 
ments of the wing and body are related to the pertur- 
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bation velocities by the Bernoulli equation 


T ae ae , 
= k a dh ag 
bp xo e ox 
1 


Ge) +) +E) © 
ow oy oz/ I 

where /p is the pressure, y is the ratio of specific heats, 
and p»,Py are the density and pressure in the unper- 
turbed free stream. 

In view of the fact that the linearized differential Eq. 
(1) is used to compute the perturbation velocities, the 
use of the complete Bernoulli equation for the pressures 
may be questioned on the grounds of internal consist- 
ency. The objection may be stated as follows. 

The use of the complete Bernoulli equation for com- 
puting pressures from the perturbation velocities ob- 
tained by solving a linearized differential equation is 
certainly correct. However, an approximation to the 
Bernoulli equation which is consistent with the accu- 
racy of the linearized differential equation for the per- 
turbation velocities would yield equally accurate re- 
sults and would lead to a less complicated expression for 
the pressure. For the case of slender wing-body con- 
figurations in supersonic flow, Broderick**? and Light- 
hill** have shown that the following approximation to 
the Bernoulli equation is consistent in accuracy with 
the linearized differential equation. 


a j Uyd¢e l or)’ l oe)" 7 
Po p ie ) Or * y 5 + 2 = 


Adams and Sears'® have obtained the same result for 
slender configurations in subsonic flow. 

An approximate form of the Bernoulli equation con- 
sistent with the linearized differential equation for the 
potential and applicable to configurations which are 
not necessarily slender has not been developed. Thus, 
it appears that a safe procedure would be to use the 
exact Bernoulli Eq. (6) for nonslender configurations. 
It should be noted that in any particular case, there is 
no assurance that use of an approximate Bernoulli equa- 
tion will not empirically result in better agreement 
with experiment or nonlinear theory. In practice, the 
approximate Bernoulli Eq. (7) or the linearized Ber- 


noulli equation 


, OY 
Ox 


bo —- p=pl (8) 


have been found to yield satisfactory estimates of the 
pressure in many cases. 

In order to avoid the complications introduced by 
using the nonlinear Bernoulli Eqs. (6) and (7), the re- 
sults of this paper will be restricted mainly to cases 
where the linearized Bernoulli Eq. (8) yields accurate 
results. The general wing-body problem in which the 
body is cylindrical in the neighborhood of the wing-body 
intersection may be treated in this fashion if only the 
difference in pressure between the upper and lower sur- 
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faces of the body or wing is computed. It is clear that 
the results obtained herein may be extended by making 
use of the exact Bernoulli equation where necessary. 

Eqs. (1) through (5) and Eq. (8) represent a mathe 
matical formulation of the linearized wing-body prob 
lem in subsonic or supersonic flow. The methods em 
ployed by various investigators will be summarized with 
particular emphasis on the range of applicability of the 
results and the approximations used in simplifying the 
problem. 
Iterative Procedure 

The most general procedure for solving the wing-body 
problem, at least in principle, is the iterative method 
first proposed by Ferrari.'* It will be shown that this 
method 1s applicable to all configurations for which the 
wing-alone and body-alone problem can be solved at 
both subsonic and supersonic speeds. The conver- 
gence of this process may be established in the subsonic 
case by appeal to the Schwarz “‘alternierendes Ver- 
fahren’’** and in the supersonic case by the existence 
proof for Volterra-type integral equations.” 

The potential ¢ is written as a sum of the form 


where each term satisfies the differential Eq. (1) and 
the condition, Eq. (4), that the perturbations vanish far 
upstream. It will be shown that the boundary condi- 
tions, Eqs. (2) and (3), on the wing and body surface, 
are specified 


respectively, are satisfied if gy" and gy" 


by the following conditions: 


Ogp™/On = —Upcos (n, x), Ogea™/On = 
—(dO¢gn On. Bom FO. od Von Rg (10) 
O¢n “ Oz =—_ =_ l oa iS, V), Ovn ‘ z= 
— (Og, D/Oe), R= ie... 64 Von Ry (11 
at £ e612... Von Ry, (12 
Ox 


It will be noted that gg, and gw? satisfy boundary 
conditions only on the body and wing, respectively, 
and hence require solution of the body-alone and wing- 
alone problems. It is assumed that a method for solv- 
ing the body-alone and wing-alone problems is avail- 
able. It is evident from Eqs. (10), (11), and (12 
that the potential functions gy" may be calcu- 
lated in order of increasing k starting with k = 0. 
Upon substituting the series, Eq. (9), for ¢ in the 
boundary conditions, Eqs. (2) and (3), and making use 


sine tO 
» ¥B 


of the defining Eqs. (10) and (11) for gw and gx", 
there result the formulas 
0¢/0zs = —Uy a(x, vy) + (O¢gp' )/dz) Rn 
13) 
O0g/On = —U> cos (n,x) + (OeHw ) /On) Re 
14) 
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Upon comparing Eqs. (13) and (14) with the bound- 
ary conditions, Eqs. (2) and (3), respectively, it is seen 
for ¢ satisfies all 
conditions of the »roblem if Og, ) Oz and O¢n on 
vanish on the wing and body surfaces, respectively. 
Oz and Ogy On approach 


that the series expression, Eq. (9), 


Itcan be shown that Ogyw . 
zero aS N= 2. 

The above iteration process will be useful for nu- 
merical calculations only if each iteration is reasonably 
simple and if the process converges rapidly. In prin- 
ciple, it is immaterial whether the wing be considered as 
extending through the body or ending at the body, pro- 
vided the same definition of the wing region is used 
throughout the analysis. The rapidity of convergence 
and numerical labor, however, are strongly dependent 
on the definition of the wing region and the method of 
starting the calculation. In certain cases, it 1s pos- 
sible to use only a finite number of steps in the iter- 
ation process, so that the method reduces to simply a 
superposition of a small number of component poten- 
tials. This can be done in the subsonic case if the wing 
‘is treated by lifting-line theory, as shown by Multhopp,* 
and in the slender-configuration case. It is also pos- 
sible in certain supersonic problems, as shown by Lager- 
strom and Van Dyke.*! All of these special cases will 


be discussed later. 


Trefftz-Plane Relationships 

The simplest relations of a general type may be ob- 
tained by considering the flow far downstream of the 
wing in the Trefftz plane (x approaching « In this 
case, O¢/Ox and 0°¢/Ox* approach zero and Eq. (1) be- 
comes 
O-¢g 
Oz? 


0" 
oy" 


(15) 


This equation is valid for both subsonic and supersonic 
flows, since the last term of Eq. (1) vanishes in either 
case. For supersonic flows, however, the region of 
validity of Eq. (15) is restricted to points near the 
traces in the Trefftz plane of the configuration under 
consideration. The treatment of inter- 
ference problems in the Trefftz plane has been discussed 
The problem of basic impor- 


wing-body 


in detail by Lennertz.! 
tance in wing-body interference is that of a wing-body 
combination consisting of wings of finite size mounted 
on an infinite cylindrical body with its generators paral- 
lel to the x-axis. 

The perturbation pressure on any surface element 
of wing or body, according to the linearized theory, is 
given by 

_ OY 


p= —pl ( ie 


(16) 
Theref« re, 
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WING -BODY COMBINATION IN THE 
TREFFTZ PLANE 


TRACE OF wine 
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©ig./ 


since Thus, the integrated pressure 
which has been experienced by any element of wing or 
body is completely determined by the value of ¢ in the 


The total lift on a wing-bodyv combina 


g(—o) = 0. 


Trefftz plane. 
tion is then given by 


(1S 


= pl, f yg cos (n,z) ds 
where cos (7,2) is the cosine of the angle between the 
outward normal to the surface of the body or the wing 
traces and the positive z-axis and the integral is taken 
over the contour enclosing the wings and body as illus 
trated in Fig. 1. If the lift distribution over the wings 
is given, so is yg over the trace of the wings in the Trefftz 
Further, the velocity Og On normal to the body 


plane. 
Thus, the problem of deter 


surface must be zero. 
mining the distribution of ¢ over the body, and hence 
the integrated body pressure, is reduced to a straight- 
forward problem in two-dimensional potential theory. 
The circumferential distribution of integrated body 
pressure on the body is, however, of less interest than 
the total body lift. This has been found in reference 
31 by the use of Green’s*' theorem, to be 


Lz = a plvg cos (n,z) ds = 


$ wings? Uo #u — gif(s) ds (19) 
where ¢, is the value of g on the wing upper surface and 
yg, is the value on the wing lower surface; on the wings 
the outward normal to the upper surface is used to de- 
termine cos (,z). The function f(s) is the velocity 
normal to the wing surface which would be produced at 
the point s on the wing by a unit vertical transverse 
velocity of the body alone in the negative z direction. 
Since (¢g, — ¢,) is equal to the circulation on the wing, 


the body lift may be written as 


Le = ply Spies M(s)f(s) ds (20 
or 
Lp _ . l( s)f(s ds ] } 
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RATIO OF INDUCED LIFT ON AN INFINITE CIRCULAR CYLINDRICAL 
BODY FOR VARIOUS WING-SPAN LOAD DISTRIBUTIONS 
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Here s is taken as the spanwise coordinate over the 
traces of the wings, and /(s) is the force normal to the 
wing surface per unit span. 

To make Eq. (21) even more concrete, it will be illus- 
trated by application to the common case of a plane 
wing centrally mounted on a circular cylinder. In this 
case, the well-known potential for the two-dimensional 
flow around the cylinder alone is 


yg = (R?*/r) cos @ 
and s = rat é = 90 
f(r) = (l/r) (0¢/08) |6 = 99 
Thus, 
f(y) = K?/r* 


and 


Lp = ply Swings T(r) (R®/r*) dr (22 


A similar result has been obtained for the antisym- 
metrical case, but it will not be considered here. 

It should be observed that if the lift is computed as 
the rate of transport of vertical momentum across the 


Trefftz plane per unit time, it is given by 
L=—-SS ply (0¢/0s) dy dz (23) 
By integrating this expression, it can be reduced to 
L= f ple cos (n,z) ds 


which is identical to Eq. (18). The derivation of Eq. 
(1S) from Eq. (23) rather than from Eq. (17) has, how- 
ever, deeper significance than is at first apparent. 
As Ward has shown in the case of slender-body theory, 
Eq. (1S) takes into account more than the linearized 
Bernoulli equation for pressures used in deriving Eq. 
(17); in fact, it accounts for, among other things, the 
quadratic approximation to the Bernoulli equation as 
given in Eq. (7). In certain problems, the momentum 
theorem basis for Eq. (18) is essential. 
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Eq. (21) provides a rapid means for estimating the 
lift induced on an infinite cylindrical body of any cross 
section due to any wing lift distribution. As applied 
to the case of wings centrally mounted on a circular 
body, Eq. (22), it was first obtained by Lennertz.' 4 
result equivalent to Eq. (22) has recently been stated 
by Zlotnick and Robinson‘ in terms of a representation 
of the wing lifting elements by horseshoe vortices 
Eqs. (21) and (22) show a remarkable property in that 
the ratio of total lift induced on the body to the total 
wing lift is relatively insensitive to the details of the 
wing loading for many practical cases. To illustrate 
this, the ratio of induced body lift to wing lift, Zz, Ly, 
based on Eq. (22) has been plotted versus the ratio of 
body diameter to wing span Rb) in Fig. 2 for the fol 


lowing three cases: 


(1) /(y) = const. 

, 
(2) (vy) = const. V [b + (R?/b)]? — [vy + (R?*/y)] 
(3) l(y) = const. (1 + (R?/y?)) 


The second of these cases is the distribution of lift for 
minimum induced drag according to Lennertz.' It is 
also the distribution obtained from slender-body theory 
by Spreiter.’. The third case is the wing loading ob- 
tained from strip theory for a wing operating in the up- 
wash of a circular body at the same angle of attack. It 
can be seen from Fig. 2 that the ratio, Lz;/Lw, does not 
vary much with these changes in wing-lift distribution. 
This result proves to be of considerable value in con- 


structing approximate theories. 
The Reverse-Flow Theorem and Its Consequences 


A reverse-flow theorem for wings of arbitrary plan 
form mounted on infinite cylindrical bodies has been 
derived in reference 31. A statement of the theorem 
has also been given by Bleakney and Morikawa.* For 
wings mounted on a body at zero angle of attack, the 


theorem states that 


/ I A pot dA = / | Apiwe dA 24 
a u a u 


where A/» is the normal force per unit area on the wings, 
w» is the specified velocity normal to the wing surface 
in direct flow, Ap; and @, are corresponding quantities 
for a case in which the flow direction is reversed, and 
integration is over the wing surfaces. w». and @, are 
defined so as to be of the same sign if the corresponding 
angles of attack in the respective direct and reversed 
flows are of the same sign. @, and w. may in general 
be arbitrarily prescribed boundary conditions on the 
wing surface. As one example of the application of Eq. 
(24), the case of a rectangular wing in the x,y plane 
Then if Ap; is 


the pressure difference across the wings for the case of 


mounted on a body will be considered. 


the wing at constant unit angle of attack to the »x,) 
plane and ae(x,y) is an arbitrary camber and twist dis- 
tribution, the total lift due to ao(x,y) is given bv 
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Loy = S Siw Apo dA = JS Sy Apr a(x, y) dA 


In other words, Ap; is the influence function for wing 
lift. In this ease, no reversal of flow is really required, 
since the rectangular wing presents the same plan form 
to the flow for either direction of flow. The reversal 
must, however, be taken into account in specifying 
the chordwise distribution of Ap. 

If it is desired to take into account slightly deformed 
body surfaces in the vicinity of the wings (such as occur 
in certain stages of iteration or superposition processes 
but with the body still at zero angle of attack at +, 
the appropriate form of the reverse-flow theorem is more 


general in form, namely 


/ / p2 we, dA : / / Pi ww dA (26 
. JW+B « SJW+B 


where p is the pressure acting on an element of surface 
and w, is the value of Og On, the fluid velocity normal 
to the surface; the integration is now to be taken over 
all surfaces of the wing and bodies, the integration con 
tours in transverse planes being similar to those shown 
in Fig. 1. An application of Eq. (26) will be given in 
the discussion of the supersonic wing-body interference 
problem. 

If the body is at uniform angle of attack az, in the 
direct-flow case and @ ,, in the reversed-flow case, the 


theorem takes the form 


cy QR + / / Ap wei dA = 
J JW U, 
Lp apn, + / / Ape = dA (27) 
et it l 


where @, and w» are again arbitrary and Lz, and Lp, are 
the lifts on the body in the direct and reversed-flow 
cases, respectively. If the configuration is rotated as 
a unit about the y axis so that the body axis makes a 
unit positive angle of attack with the x,y plane, then 
dp 1, and @, = cos (v,z). As in the Trefftz-plane 
theorem, cos (7,2) on the wings in this case is the cosine 
of the angle between the outward normal to the upper 


wing surface and the z-axis. Then 


a = / / Ape cos (1,2) dA 
J JR 
~ / J Ap, —.dA (28) 
Sa Uo 


simply the total lift on the 


W. and 


The terms on the left are 
wing-body combination in case 2 for arbitrary 


az, SO that 


Lu b> Lp-ep + Jf. Ap, 7 dA 


In this instance, 4p, and Lz, together form the influence 
function for the total lift of a wing-body combination 


(29) 


N 
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with arbitrary wing camber and twist, and arbitrary 
body angle of attack. 

l and @; 


By setting w. = O on the wings, dz 


0, (w,/ Uy) = cos (n,z) in Eq. (28), the result 1s 


Lp / / Ap cos (n,z) dA 30 
° Si 


This result has a simple interpretation, if it is recog 
nized that the term on the right 1s the lift contributed 


by the wings, so that 
Lp Li 3] 


Lx», is the lift induced on the body by the wings each at 


angle of attack cos (7,2) in case 1, while Ly, is the lift 
induced on the wings, which are at zero angle of attack, 
by the body in case 2. If the wings are in the x,y 
reduces to the simple 


the 


plane, cos (7,2) l, and Eq. (51 
statement that the lift induced on the body by 
wings at unit angle of attack is equal to the lift induced 
on the wings by the body at unit angle of attack in re 
verse flow. 

The use of the reverse-flow relations becomes par 
ticularly simple if the wing has a spanwise axis of sym 
metry, or if one is dealing with slender-configuration 
theory or lifting-line theory, because in these cases the 
reverse-flow solutions are not really necessary to use the 
relations. These cases will be considered later. 

Another useful result which can be obtained from the 
reverse-flow theorem is a variational principle which is 
of value in setting up approximate theories. The pro 
cedure used has been explained in some detail in refer 
ence 37. The essential result is the variational equation 


x ‘ Og» 
0 Pi = prw ~ po® dA = 0 32) 
- W+B on 


where the barred quantities correspond to reverse-flow 
conditions as before. Thus, for first-order variations 1n 
p, and fp» from the exact values corresponding to the 
specified boundary conditions, this expression deviates 
from zero only by higher-order terms. The exact solu 


tions correspond to (O¢)/ 0”) @,, and (Og2/On) = 


W,,. Ogs/On in this equation is to be regarded as the 
correct value corresponding to ps. @,, and w,, are to 
be regarded as the values required by the boundary 


conditions of the problem; they do not vary. If pe is 
exact, (Og /On) = wW,,; if pe varies from the exact value 
of pressure, so does Og./On vary from w,,. If py 1s 
taken to be an approximate solution for @,, appro- 
priate to a constant unit angle of attack in reverse flow, 
and if the solution in direct flow is taken to be py = 
cfo(x,v,z) on the surface of the wing-body combination 
where c is an arbitrary constant which must be deter- 
mined and fs(x,y,z) is some plausible distribution of 
pressure on the boundary in the direct flow case, it 
can be from the 
theorem that the value of total lift obtained from 


shown Eq. (32) and reverse-flow 
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— _— 
— w,,) pj)dA = 0 (3: 
B+Ww\on 


also has the property of being stationary with vari 
In this 


“~~ 


ations in fp; and p. from the correct values. 
sense, P; is a proper weighting factor to use in satisfying 
the boundary conditions on the average as is often done 
in approximate solutions. If /; is exact, then Eq. (33) 
gives the exact answer for the total lift, regardless of 
what is assumed for the function fs representing po. 
This is readily seen, since if the value of Og./On calcu- 
lated from the assumed distribution /: is given by 


Og2,/OnN = C¥2(X,V,2 


Eq. (33) gives 


(SP Lmoul( ff, 


The total lift is given by 


Lwe = € / / fo cos (n,z) dA (35) 
JW+B 
so that 


Lwe = (/ I Wr PI aA) (/ / gop, aA )x 
W+B / JW+B 
/ / fo cos (n,z) dA (36) 
JW+B 


At this point it should be noted that Ll cos (”,z) = 
@,1, since case | in reverse flow corresponds to unit angle 
Also, because 


esp ) (34) 


of attack of the wing-body combination. 
of the reverse-flow theorem as given in Eq. (26) 


/ / fo = dA = / / £2 Pi dA (37) 
Jw+Bs Us, JW+B 


since f2 and go are, respectively, the correct p. and 
Og./On for some case in direct flow, even though not 


the one of interest. Therefore, 


Lwe = JI Wy pi dA 
W+B 


But by the reverse-flow theorem Eq. (26) 


/ / W,,pidA = 
W+B 
/ / Uy cos (n,z)py dA (39) 
JW+B 


where p» is the exact solution for the specified w,,. 
Therefore, /; is the proper influence function for lift of 
the wing-body combination and Eq. (38) thus gives 
precisely the correct lift. 

This procedure can be generalized to permit the use 
of several weighting factors in combination with a series 
of terms for f.. As shown in reference 37, proper 
weighting factors to make lift, pitching moment, and 
rolling moment stationary can be chosen. For the 
special cases lifting-line theory and slender-body theory, 
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the use of reversed flows becomes unnecessary and 
simpler variational methods can be derived. 
discussed in the section dealing with these theories 


SUBSONIC WING-Bopy INTERFERENCE 
Multhopp’s Formulation 

The subsonic problem of a lifting surface mounted 
on an infinite cylindrical body of any cross-sectional 
shape will now be formulated, along the lines indicated 
by Multhopp,* and the solutions of this problem for 
certain limiting cases indicated. The results of such an 
analysis can, with modifications such as those given 
by Multhopp, be applied to give approximate results 
for finite bodies as well. 

Consider a thin wing mounted on an infinite cylin- 
drical body, whose generators, which lie approximately 
in the direction of the x-axis, make a small angle a, 
with the x,y plane. The free-stream velocity ly is 
taken to coincide with the x direction, and the wing 
surface is assumed to lie approximately in the x,y plane. 
The angle between the wing surface and the x,y plane at 
any point is denoted by a(x,y). The entire configura- 
tion is assumed symmetrical about the y,z plane. It is 
convenient to consider the flow to be composed of two 
parts: first, the flow around the isolated body at angle 
az and second, the flow produced by the wings in the 
presence of the first flow. It may be noted here that 
the first flow is simply a two-dimensional flow past 
the body cross section. The second flow must be de- 
termined in such a way as not to violate the boundary 
conditions on the body. The velocity potential, gp, 
describing the isolated body flow field is determined 
by the condition 


V*o = 0 in the flow field (40 


Og /On = —U va cos (n,2) on the body surface (41 


where cos (7,2) is the cosine of the angle between the 
outward normal to the body surface and the 2-axis. 
The velocity potential, ¢:, describing the second flow 


is determined by the conditions 


A*9, = 0 in the flow field 
O¢) : ° ° 
> = — UU a(x,y) — 0¢g/0z on the wing surface 
(42 
O¢1 - e 
= 0 on the body surface (45 
On 
gi = O atx = -—« 14 


In addition, the pressure difference between the upper 
and lower surfaces at the wing trailing edge must vanish 
according to the Kutta condition. The linearized per- 
turbation pressure is given by 


—p = pU(d0¢/O0x) (45) 


so that at the trailing edge 
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(Ogi, /OxX) — (Ogy,/0x) = O 16 


where the subscripts « and / are used to identify upper 
and lower surfaces, respectively. 

The 
tions on the solution to the combined problem, since 


conditions on gp and ¢ give the proper condi- 
¢ = go + ¢1 Satisfies the equations previously given in 
the general formulation, Eqs. (1-5). 

The solution for gp can readily be obtained if the body 
cross section can be mapped conformally onto a circle 
or a Straight line, which can easily be done in a wide 
variety of cases. The problem of finding ¢; is consider- 
ably more complicated and, in general, can be done 
With further 


approximations, however, it is possible to find ¢g in 


only for some simple limiting cases. 
many cases of practical interest. The g; problem is 
essentially that of a wing mounted on an infinite body 
at zero of attack with an angle of attack equal to 
+ (1/U%) (Og For convenience, define 


iS). 


ar,y 


l O¢go(V) 
l QR Oz zs=0 


f(y) = (47) 


The function f(y) is then the vertical velocity produced 
along the wings by the isolated body in two-dimensional 
flow for unit vertical velocity of the body. The down- 
wash, ww, on the plane z = 0 produced by the lift on 


the wings is, according to lifting-surface theory, '® 


1 oO aa 4 dé dn 
ww(x, vy) = u1(E,n) X 
tor Ov Rwy — 1 
/ : 
f ¥ if — £° + y= 978 
41+ 5 (48) 
l x—¢€ 
where “4; = (Og\u/Ox) (Og,,/Ox) and w; = (O¢)/O0z). 


It should be noted that the lift per unit area on the 
wing is given by /(x,y) = plot, which is directly pro- 
portional to u;. The downwash produced by the lift 
on the body generally involves more complex consider- 
ations which will not be gone into here, because the re- 
sult is usually so involved that no calculations can be 
made with it. Since, however, a Green’s function for 
the problem exists, it is possible to set down an oper- 
ator A(x,v; £7) such that the additional downwash, 
Wiz, produced by an element of wing lift in the presence 


of the body is 


WB = / J dé dy u,(£,n)A(x,y; &,n) (49) 
Ru 


Wiw + Wp. 


ary condition on the wing, Eq. (42), 


= 7 an 
didn X 
tr Oy J Ru 
Vix — £2? + (y— )*| 
i+ oogneeanint tna 
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/ | dé dynuy(é,n) A (x,V; &0) 
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Thus, w, = In accordance with the bound- 


Uja(x,y) + Uvanf(y) = 


ui(&,n) J 
beats ( 


(50) 
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From this integral equation in ™, the wing lift distribu- 
tion can, in principle, be solved for any given case, sub 
ject to the Kutta condition, w; = 0 at the wing trailing 
edge. 

One particular case in which the determination of the 
correction kernel A becomes simple occurs when x ap 
proaches © ; in other words, when the downwash in the 
Trefftz plane is required. In this case, Eq. (50) be 
comes 


. 


1 oO / "dy / ane ; - 
Ww, = uy(&,n) dé + 
27 OV J-by — 9 Sx0(7 
/ I. 


where x) and x; are the x coordinates of the leading and 


dé dnuy(E.n)K(@©,y; En D1) 


trailing edges, respectively, while ) is the semispan of 
the wing. In integrating from —b) to #b, it must be 
understood that the region of the body itself is excluded. 
The lift per unit span may be defined as 


pl / 


and the circulation as 


l(n) = ui(&,n) dé 


I(n) TE i at 
‘ey * —. = uy(é,n) dé 
pl J xol(n 


Then, 


|) & 2” ie 
WwW) = dn so 
2r Oy J-byY — 7 


/ I dé dnuy(é,n) K( ©, vy; &) (52) 
J Ru 


Integrating the first term by parts and then differen- 


tiating under the sign gives 


l ie dV dn dn 
whic 
27 'y~*7 
/ I dé dnuy(é,n) K( ©, vy; En (53) 
J Ri 


The exact expression for w, in the Trefftz plane is, how- 
ever, knownin another form from the work of Multhopp.* 
If 7 = g(t) is the conformal transformation which takes 
the body in y,z plane (¢ = y + iz) into the z axis in the 
¥,2 plane (7 = 7 + 72), then the wing-body combination 
in the ¢ plane is mapped into a wing-alone in the / plane; 
the vertical slit along the = axis representing the body 
is necessarily a streamline, since the flow 
perpendicular to the 2 axis vanishes by 
Thus, the boundary condition on the body 
In the transformed plane then 


] f (dP /d4) di 
Qn. b y —_ n 


The downwash velocity in the physical plane is given 


component 
symmetry. 


is satisfied. 


(54) 


“ 


ol 
| 


by 
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WwW, = @, R.P. (di/dt) (55) 


where R.P. denotes “‘real part of,’’ o 


; *b LIK / js) Js 
| (5) l / (dI'/dn) dn (56) 
dt] 2r 5 ¥-7 

It should be noted that this expression for w in terms of 
@ is an approximation based on neglecting the sidewash 
velocity on the transformed wing. Multhopp points 
out that this effect should usually be small 
wing with a body symmetrical about the plane of the 
A check of some results obtained 


for a plane 


wing the error is zero. 
with this approximation for a body unsymmetrical about 
the wing plane is given later. 

Since the discontinuity in potential (¢, — ¢,), and 
thus the value of [, are not altered in a conformal 


transformation, 
r(n) = Tla(n)] (57) 
and 
(d0' /dy) dq = (dV /dn) dn (5S) 


Substituting this relation into Eq. (56) 


di\ 1 “> (dV /dn) dy 
w, = R.P. ) / (59) 
= 2r 6 V7 


and 4 are explicit functions of y and y, respec- 
The velocities of the flow around the cylinder 


where ¥ 
tively. 
alone may also be given in terms of the transformation. 


In fact, the upwash function, /(y), is given by 


f(y) = R.P. (di/dt) — 1 (60) 
Thus, 
R.P. (di/dt) = 1 + f(y) (61) 


and 


wi(o, y) = (62) 


[1 + f(y) ] he dn) dn 
rt FoF 


27 


Comparing this result with Eq. (55), it is clear that 


the body effect on the downwash is given by 


Rw 
| ar j{1 + fy) ] : 4 
= = ni dn 
2a b dy - 7 = nh 
For lifting elements located at substantially a single 


chordwise location, say x = 0, the downwash at that 
location is one-half the value in the Trefftz plane, so 


that 
I “? (dT /dn) d 
tar sii 6 ¥~7F 


Wip( ©, ¥) 


1 (O0,¥) = [ 


where, as before, it is understood that the region of the 
body is excluded from the range of integration. Also, 
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*b oan 
Wi B(O,Vv) = : * | = ¢ * —_ _ . up 
; bdnl ¥- 4H y—a! 


iy J dé, dnuy(é,n) A(O0, v; 0, 69 


It is convenient at this point to proceed to the discus- 
sion of special cases for which the complete formulation 
and solution of the wing-body problem is possible. 


Slender Body Theory 
The slender-configuration theories of Munk*® and 
Jones’ as extended to slender wing-body combinations 
by Spreiter’ and Ward® offer a particularly simple case 
with which to begin. In these theories, the assumption 
is implicitly made? that as far as the calculation of in- 
duced downwash is concerned, each point on the wing 
is in the Trefftz plane of all upstream points, and in- 
finitely far ahead of all downstream points. Mathe- 
matically, this results from the approximation in Eq. 
(48) 
Vv (4 ey a n) a” x ‘a 


x—é€é x 


Ser 


Ste 


which would be expected to be valid if vy — n is always 


small compared to |v — &, 1.e., for slender configura- 


tions. In this case, it is convenient to define 


7 


v(x,V) = / U,(&,V) dé (67 
J x \ 


This is 1/pl) times the lift on a strip from the leading 


edge up to the point x. The downwash is given by 


[1 tio) fe "(¥ on) ' 
dy (OS 


Da 
-7 


w(x,Vv) = 


which may also be written as 


1o fe : . 
W(X,V) = = — di 69 
Zr Oy J—-b ¥ — ij 


This must satisfy the boundary condition, Eq. (42), so 


that 
+ f(y)] Le (Ov On) , 
~~ an 


70 


a(x,v) + apzf(y) | ” (Ow /On) - 
U = - dyn (71 
La Fy 2rjJ-b ¥—7 


The most convenient way to solve this integral equa- 
to the transformed co- 


Upa(x,y) + U QR, (vy) = 


or 


tion is to convert it entirely 


ordinates, jand 7. This gives 


fa[xv(¥)] + anfly(*) ]) 7 l | ee On) dq 
{ 1 + f[y(¥) ] § w+ F-—% 


-) 
i= 


U" 


t The explicit formulation usually given is quite different from 


this one, but is equivalent 





sta 


the 


whi 
of t 
tot 
pla: 


lift 
on 1 


Nov 


and 


SO 
Stat 
laps 
ont 


fon 


ISCUS- 


ation 


and 
tions 
case 
ption 
f in- 
wing 
d in- 
ithe- 

Eq. 


(66 


Vays 


‘ura- 


(67 


ding 


(OS 


(69 


70 


rom 





WING-BODY 
This equation may then be inverted directly by well- 


known methods to give as the solution 


oY — 2 "> falxn(a) |} VB — 4 - 
E é = — dq (73 
Ov 7 V b - — y - } ! ae 7\n | ‘ = — 


A simple and important case occurs when a(xv,V) 1s con- 


stant as for an uncambered untwisted wing, and when 
the wing and body are at the same angle of attack, ap, so 


that 
a V,V = AR a 74 


”) then becomes 


| f° (O%/d%) da ” 
( a 2 3 (40) 
27 ( ‘ile 


[his corresponds to the case of a slender wing alone of 
The solution for this 


~! 


Eq. 


span 26 at angle of attack ap. 
case is Well known and is 


V(ix,y) = 2Uyay Vb? — ¥ (76) 
It should be noted that 6 and ¥ are functions of x. The 
total lift in this case is given by 
L = wpa max. (77) 
or 
L — . 
= 2a max . d) 
q 


where Dna, is the value of 4 at the widest cross section 
of the wing. It will next be shown that Z is in fact the 
total lift on the wing-body combination in the physical 
plane. 

According to the Trefftz plane relation, Eq. (21), the 
lift on the wing-body combinations in terms of the lift 
on the wings is given by 


L= pl; / r(y) [1 + fiw) | dy (79) 


b 


where the body is excluded from the range of integra 
tion. The lift on the wings in the transformed plane is 
given by 


*h 
L = pl / r'(¥) dy (SO) 

h 
Now, from the transformation equation and Eq. (61) 
d¥/dy = R.P. (di/dt) = 1 + f(y) (SI) 

and in accordance with Eq. (57), 
r(y) = Flv) ] (S2) 


so that L = L. 


Stated: If in conformal transformation 


Thus, the following rule may be 
which col- 
lapses the body of a slender wing-body combination 
onto a vertical slit in the plane of symmetry, the trans- 
formed wing span is 24, the lift of the wing-body com- 
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bination is equal to the lift of an isolated slender wing 
of span 26. 

Applying this to the simplest case of a wing of span 
2b mounted on the centerline of a circular cylinder, the 


conformal mapping function is 


t— (R 83 
and 
R.P.(dt, ¢ + (R* when : QO (M4 
Also, 
R2/ 
and 
h, b — (R2/b 


Substituting this into Eq. (7S) gives, for the lift of the 


wing-body combination, 


L = rpl*ao|b 


(R y)] SO 


which is in agreement with the results of Spreiter. Sub 
stituting for 6 and ¥in Eq. (76) leads to 


Y(x,V 2U ya V [b (R?/b) | y R?y) |? (86 


which is equivalent to the results of Spreiter and Len 
nertz for the spanwise lift distribution. It should be 
noted that this expression is identical to 


2UvaV [|b + (R?/b)] y + (R?/d)] S7 


W(x,y 


which is more commonly given. 
If the wing is not centrally mounted, but is located a 
distance below or above the midplane of the body, 


the coordinate of the tip is 


t=b+ th 
and 
a (R2h) /(b? + hr? 


From Eq. (7S), the lift of the configuration is then 
L g = 2rrayb?} | - {R (b> + h li SS 


Still another interesting case in slender-body theory 
may easily be solved, namely that of a wing centrally 
mounted on a body which is an elliptic cylinder. It is 
convenient here to use the complex variable u z+ 


! 


iy rather than ¢/ = y + 7s. For this case, the trans 


formation, which maps the body on to a vertical slit, is 


A? + B?] 


(S9 


n= [1/(A — B)] [Au BWV u 


where B is the axis of the ellipse in the y direction and 
A is the axis of the ellipse in the z direction. The trans 
formed semispan of the wing is obtained by substituting 
u = 1b into this expression as 

bh = [1/(A — B)] [Ab 


BWV A* — B? + b*] (90) 


According to Eq. (7S), the lift of the wing-body com 


bination is 
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L/q = 2xa0[(Ab — BV A? — B? + 6?)/(A — B)}? 
(91) 


This result is in agreement with those of Nonweiler'’ 
who used Spreiter's’ approach and Kahane'! who used 
methods similar to Ward’s.* 

The virtue of the rule given here is that the total lift 
on any wing-body combination for which the necessary 
conformal transformation can be found may be calcu 
lated without detailed calculation of the pressure dis- 
tribution. 

The relationship between the solutions obtained from 
slender wing-body theory and those obtained from the 
classical theory of wings of minimum induced drag also 
merits some attention. It has often been observed that 
the results are intimately related. 
Jones theory for a slender wing at constant angle of 
attack gives an elliptic distribution of lift which is the 
the Spreiter 


For example, the 


distribution for minimum induced drag; 
theory for a slender wing centrally mounted on a circu 
lar body gives a distribution of lift identical to the 
Lennertz solution for minimum induced drag for that 
case. The relationship between the two theories is 
easily seen by setting up the slender-body theory in 
variational form as in Eq. (32). The variational equa 
tion is, for the wing at constant angle of attack, 


*b *b 
if l'(y)w(y) dy — 2 (vy) Uvao X 
b b 


[1 + fy) dy, = (0 (92) 
where w(y) is given by Eq. (56) and integration is over 
Carrying out the variational process 
The minimum in- 


the wings only. 
leads back to the integral Eq. (75). 
duced drag problem as set up by Lennertz and, in more 
general form, by Pepper,’ requires that the induced drag 
be made a minimum subject to the constraining condi- 
tion of constant lift. Since an infinite body at zero 
angle of attack can contribute no drag, the induced 
drag arises from the wings only and is given by 


U. {° 
D= . - / l'(y) wy) dy (93) 
21 6 —b 


The lift of the wing-body combination is from Eq. (21) 


*b 
py | I(y) [1 + f(y) ] dy (94) 
—b 


where, as before, the integration is over the wings only 
The variational problem of finding the minimum in- 


Ny = 


duced drag reduces to 


b 


l 
5 (vy) w(iy)dy +A 


12U,,J_, r(y) X 


{ i 
[1 + f(y)] dye = 0 (95) 
where A is an undetermined Lagrange multiplier which 
must be introduced to stipulate the condition of con- 
stancy of lift. The variational equations of the two 
problems are thus identical, except for thé constant 


TICAL SCIENCES MAY, 1954 
multipliers of the second integrals in Eqs. (92) and 
(95). It follows then from the form of these equations 
that the lift and downwash distributions can differ at 
most by an arbitrary multiplicative constant. This 
arbitrariness stems from the fact that the solution of the 
minimum-induced-drag problem cannot be expected 
to give a value for the lift in terms of boundary condi- 
tions on the wing, but only a relation between lift and 
drag for the particular distribution which gives mini- 
mum drag. It is now possible to state the foliowing 
rule: For any slender wing-body combination having 
a evlindrical body with constant angle of attack of the 
wings and with an equal angle of attack of the body, the 
distribution of lift over the wings is identically the one 
giving minimum induced drag for a wing mounted on an 
infinite evlindrical body of the same shape. 

Examining the results in additional detail, it may be 
seen from Eq. (95) that the Trefftz-plane downwash 
distribution for minimum induced drag is always 


w(y) = woll + f(y)] OG 


is an arbitrary constant. In order to satisfy 
the slender wing-body 


The lift in 


where w, 
the boundary 
theory, zw; must equal Lyay (see Eq. (56)). 


condition in 


the minimum-induced-drag problem can be determined 
in terms of w) from momentum considerations (see 


reference 39), as 
L = pl oto >, 97) 


where Y is determined from a momentum integral in 
the Trefftz plane. Thus, the results of the theory of 
wings of minimum induced drag may be applied with 
ease, since Y is the quantity which must be determined 
to obtain the induced drag in terms of the lift. For in- 
stance, in the case of a wing centrally mounted on a 
circular body, which was considered previously, the 


theory of minimum induced drag gives*® 
>> = xl[b — (R2/b)]? (98) 


From this 


L = rpUPao[b — (R?/b) }? (99) 


which is in agreement with Eq. (85). 

As a further application of the relation between the 
solution for minimum induced drag and the slender- 
body solution, consider again the case of a wing mounted 
on a circular body a distance d below the central plane 
of the body. As pointed out previously, the solutions 
obtained in this paper for such cases are approximate 
because the transformed wing has some dihedral and 
the mean sidewash at any point is not zero. Pepper 
has given the exact solution for the minimum induced 
drag problem in this case. Some of his results for > 
are shown in Fig. 3, and compared with © from Eq. (SS). 
It may be seen that the error is small. 

A somewhat more complicated problem of some in- 
terest may be solved in this same way, namely, the case 
of a wing-body combination with a gap of span s on 
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WING-BODY 


each side between the wing and the body. Consider 
first the case of a wing alone with gap (the total gap in 
this case is 2s), a case which has been treated by Ble 


viss and Struble.*” The value of Y for this case is** 


= K’ — 2k’ 
= b? ) + 5? 
~ r| ( kK’ | 


where A’ and £’ are elliptic functions of the first and 


(100) 


second kind, respectively, of argument k’ = V 1— (s, 6)*. 
The lift of the slender wing with gap is then 


(K' — 2E! 
L= restau] 6( : = ) + | (101) 
kK 


The case of a gap on a slender wing centrally mounted 


on an infinite circular cylinder may readily be obtained 


by making use of the earlier result that the lift of the 
wing-body combination is the same as the lift on the 


wing alone which is obtained by conformally trans- 


forming the body to a slit on the vertical axis. The 


transformation is again 

i= t— (R*/i) 
Thus, 
bh = b — (R?/b) 
j= 5 — (5/5) 


Substituting ) and § for b and s in Eq. (10) gives for the 
lift on the wing-body combination with gap 


se *) K’ — E’ 
_= Uv?ao\( bo — = + 
(s _ e)' (102) 
where the argument of K’ and FE’ isk = V1 — (5/5)? 


This result was first obtained by Harold Mirels of Cor- 
nell University in a paper as yet unpublished, using 
methods analogous to Ward’s.* 

The reverse-flow theorem leads to particularly simple 
results for slender configurations; in fact, no flow re- 
versal at all is required. Results identical to those of 
the reverse-flow theorem for this case can be obtained 
directly by considering two distributions y,(x,y) and 
¥o(x,y) corresponding to the angle distributions a;(x,y) 
and «as(v,vy), respectively. Then by Eqs. (70) and 
(69), the integral equations for y7 in the two cases are 


Uvar(x,y) + Usan,(x,y) f(y) = 
1 oO 6 : 
——— dj 
2m OV bv— 7 


+ Barred quantities in these equations correspond to con- 


(103) 


and 


formally transformed values 
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COMPARISON OF APPROXIMATE AND EXACT 

THEORETICAL LIFTS FOR SLENDER WING-BODY 
COMBINATION UNSYMMETRICAL 
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Fig. 3 


Uyae(x,y) + Uoap,(x,y) f(y) = 


l a) z y2 li 
or dyJ-17 —-4H ” 


Multiplying the first equation by y2 and the second by 
y, and integrating over the span gives 


. 


(104) 


*b 
Uo / Yr [ay + ap, f(y) ] dy == 
oY b 


- is o f° Lil) 
/ yo(v) / = : djidy (105) 
2a b OV . b y - n 
and 
*h 
Vy, J V1 [as a ap, f(y) |] dy = 
b 
| 7%, Pa “ vo(n) 
(y) dj dy (106) 
x LH ~ f=. silk 


where, as usual, the body is omitted from the range of 
integration. The terms on the left-hand side may 
easily be shown to be equal by integrating parts, ex- 
changing the order of integration, and taking account of 
Eqs. (68) and (69). 


-. s 
/ Yo [a + ap, f(y)|] dy = / WX 
- hb 


[a2 + anof(y)| dy 


Thus, 


(107) 


This is valid for any two angle-of-attack distributions 
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ay(x,v) and ao(x,v). A particular case which will be 
required in mei developments in the case a2 = agy = 
1, so that yY2 is the circulation distribution for unit angle 


of attack of wing and body given by Eq. (87). Then 


" 
| 2U,V 8? —- v [ar + ap, f(y)] dy = 
b 
*h 
/ vill + f(y) ] dy 
. b 


The integral on the right-hand side is the total 1 pl’ 
times the lift forward of the location under consider- 


(10S) 


ation, L(v), by Eq. (79). Thus 
(109) 
b 


*h 
L(x) = 2p: Vi? Vv? [ar + ap, f(y) | dy 


This expression for the special case of slender wings 
centrally mounted on a circular body with zero body 
attack has been Morikawa and 


angle of given by 


Puckett. ! 


Low-Aspect-Ratio Wing-Body Combinations 
Slender-wing theory, although attractive in its sim 
plicity, proves to be accurate only for wings of very 
small aspect ratio; for instance, the error in slope of 
the lift curve of a triangular wing is already of the order 
of i0 per cent at an aspect ratio of unity. The error in 
the predicted pitching moment of rectangular wings is 
even more serious. In reference 13, an approximate 
theory was developed which gave satisfactory accuracy 
for many purposes up to an aspect ratio of four. The 
theory has been extended to wing-body combinations 
in reference 12. The essential approximation in this 


analysis is the approximation in Eq. (48). 
eo Fy oe 4 


+ V(x —¢)?+ 7) 


x—é& 


L tix 
J (110) 
| 


where @ is an average local spanwise distance which is a 
function of x only. This represents a somewhat more 
accurate approximation than Eq. (66), which is the ap- 
proximation used in the Jones theory. The approxi- 
mation of Eq. (110) approaches the Jones theory ap- 
is large compared to the 


x —é& 


proximation when é 
A good value for 8(*) has been found to be 


span, B. 
Thus, from 


the local exposed semispan of the wing. 
Eq. (48), the downwash due to the direct effect of the 


wing lift is 


dé d 
Wiy(X,¥) 12 f., : = u(&,n) X 
v — ¢) +V(x— £)? + B? 
i+, cng ee 
x—€& 


The term in braces is a function of x only so that 
w, may be written as 
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ee & 
|) (.V,4 — / de — 
tor J x { 
c— & + V Ww — &)* + B* 
x 
(x _ &) { 


O / (e 141(&,n) mas 
dn 112 
Oy be» YY — 
The term braces describes the variation of down- 
wash with distance downstream of the location of a 
lifting element. If it is assumed that the body con- 
tributions to the downwash behave in the same way 
as the direct contributions of the wing lifting elements 
as far as this variation is concerned, then from Eqs. 
(112) and (69 
+ apf(y) = 


a(xv,V) W(V,V) = 


0 f*? uls& a), 
= an 115) 
Oy é } ee n 


To reduce this integral equation in two variables to an 


integral equation in the spanwise variable only, the 


method of spanwise integration with weighting fac- 
tors!® * 4° is employed. The pressure distribution for 
a constant angle-of-attack distribution on the wing in 
reverse flow, used as a weighting factor, has the prop- 
erty of providing a stationary (in the sense of the calcu- 
variations) integral for total lift, as indicated 


In the present case, the approximate span- 


lus of 
previously. 
wise pressure distribution in reverse flow is taken to be 
the result obtained from slender-body theory, Eq. (76) 
for the wing and body at constant angle of attack. 

po(x,y) = const. Vb? — (114) 


Multiplying Eq. (113) by this factor and integrating 


over wing span gives 


mi) 
/ Uy [a(x,y) + azf(y)] Ve? — FP dy = 
b 
l 7 : %=— < B- 
/ iy J ded + x 
Lor v >| 2.x — &) ( 
To = @ ‘ dj 
vo? — a | : 11(&,7) 
oy ¥—7ij 
If the order of integration is changed, this becomes 
*b 
j 
b 


| Px x —¢ V(x ge + 8 
| diji+" 5 tt Lx 
lr x | 2(x el &) f 


Str 
aa 
él 
| 
Str 
+ 
Ww 


(115) 


0 [a(x, vy) + apf(y)] b? — y dy = 


*b . = 
=m aes O b d : 
/ dy Vb? — ¥° / : uy(&,n) (116) 
b Oy J-bY — fj 
By virtue of the reciprocal relation Eq. (107), since 
V b? — F* is the distribution of p for unit angle of attack 
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WING-BODY 


of wing and body, the last double integral is reducible to 


avery simple form, namely, 
| ; T° n° Oo t dn 
/ dy Vor - y uy(&, = 
Ie ; Ov/ -by¥ — fj 
/ uy(&jn) 11 + f(m)] dn (117 


It is apparent from Eq. (79) that the integral obtained 


AS 
oT) 


on the right-hand side of the last equation is simply the 
circulation per unit chord on the wing-body combina 
tion It is convenient to define this quantity as 


i. 7” si 
g(x) = (x,V + I(y)] ay 
OX. ini 


(11S 


Also, the term on the left-hand side of Eq. (116) may 


be defined as 


R(x = / la X,Y + api v IVA - y dy 
b 


. 


(119 


This latter quantity, as has been shown in Eq. (10S 
is simply the circulation on the wing-body combination 
up to a distance x from the origin according to slender 


configuration theory. Eq. (116) then becomes 


or 


This one-dimensional integral equation is identical in 
form to the one derived for low-aspect-ratio wings alone 
in reference 13. The only differences are that k(x 
in this case is the Spreiter-Ward value for the wing-body 
combination, rather than the Jones value for a wing 
alone, and 8 is approximated as the exposed span of the 
wing rather than as the span of an isolated wing. 
Methods for solution of Eq. (121) for circulation per 
unit chord g(*) by collocation have been given in refer- 


ence 13. 


Lifting-Line Theory for Wing-Body Combinations 

The theory of wing-body interference for high-aspect- 
ratio wings mounted on infinite cylindrical bodies ac- 
cording to Prandtl lifting-line theory has been defi- 
nitively presented by Multhopp.* It will be desirable 
here to briefly indicate how the lifting-line approxima- 
tion may be arrived at from lifting-surface theory when 
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a body is present. Reissner** has shown that for wings 


alone, the approximation 


>) 


V (x — €)* + (y n y n 122 


in Eq. (48) for wi, leads to the Prandtl result. With 


this approximation, 


1 oO dé dn | y—n) 
Wiy- = i+ - u,(&,n 123 
ws OV . JRwV— 7 { \ - £ { 


If the wing is considered to be rectangular for simplic- 


¢ 


ity, this may be written 


l i / ag dn Ou ay n ' 124 
bit fr J JSRwV — n On ' \ ct 7 


or 


"Or on | “m(E,V) | ; 
iy = dyn + dg 125 
ly bv - n Ix 2 \ £ 


TT ew 


where c is the chord and 


(yn) = / u,(&,n) dé 126) 


The first integral in Eq. (125) will be recognized as the 
downwash of a trailing vortex sheet, of strength dl dn, 
at the front edge of the sheet, while the second is the 
downwash of a two-dimensional vortex sheet of strength 
ui(x,v) along the chord. The second term does not 
depend on lifting elements anywhere except at vy. 
Thus, to the lifting-line approximation, the presence 
or absence of a body should make no difference in this 
second term on the wing. It may be observed that the 
first integral which gives the effect of the trailing vortex 
sheet is one-half the Trefftz-plane downwash due to 
lifting elements on the wing itself. The corresponding 
value for a wing plus infinite cylindrical body is given 
by Eq. (64). Substituting the latter expression for the 
downwash in place of the first integral in Eq. (125) and 
setting the downwash equal to the value required by 
the boundary condition, Eq. (42), gives the following 
integral equation 
, l 
Uofa(x,y + ap f(y) ] = 1+ (vy > 4 
Lor 


“(dT dn) d l eS 4,(&,V 
/ a +. / “_ dt (127 
sb =F 27 2s=—¢ 


Multiplying by the lift distribution in two-dimensional 
reversed flow as the weighting factor and integrating 


over the chord gives 


ve I( 2)T x 
Uo [a(x, vy) + an f(wl@/ dx = 
2 icin) — 2 
7h , : 
] , XC (dl /d ) dn (y 
{1 + f(y) / = : = (128) 
ior ° 2 i F—F 2 
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If ais a function of y only, this becomes 
: l 
Uvla(y) + azf(y)] = [1 + f(y)] x 
4a 
"b Lan ) : 
(dl dn) dn I'(y) 
enw ae (129) 
. 3 =—9 1C 


This is the integral equation for a wing-body combina- 
tion according to lifting-line theory. It may be solved 
by any method which can be employed for the lifting- 
line equation in the case of a wing alone. Multhopp* 
has given solutions for several particular cases. 

In order for the lifting-line formulation of the wing- 
body problem to be valid, not only must the wing aspect 
ratio be large, but, also the ratio of body diameter to 
wing chord must ke large. This is because there is a 
rapidly varying angle-of-attack distribution due to 
flow around the body in the vicinity of the body, and 
in such regions, lifting-line theory is not accurate. 
Multhopp points out that, actually, the induced down- 
wash at the 3/4 chord position of the wing (the cen- 
troid of the chordwise distribution of lift in reversed 
flow, and therefore the centroid of the weighting func- 
tion) will be greater than one-half the Trefftz-plane 
downwash. Thus, the predicted lifts near the body 
will be too high. Further, the flow field of a lifting line 
in the presence of the body is actually readily determi- 
nable only in the Trefftz plane, Eq. (62), and at the line, 


MAY, 1954 


Eq. (64). At other longitudinal locations, it is very 
difficult to obtain. 

A variational method which allows relatively rapid 
determination of the total lift of a wing-body combi- 
nation according to lifting-line theory has been given in 


reference 4. The variational equation is 


"bore *} 
6 / _ dy — 2U | r ary) + aR f(y) | dy as 
a bac — b 
*h 
/ rw ay | = 0 (130 
b 


Oa 


where, as usual, the region of the body is excluded 
from the range of integration and w is given by Eq. 
(64). 


gives 


Upon performing the variational operations, this 


PA i 


oc 


+ w= } a(y) + ap f(y) f Ui 13] 


which is identical to Eq. (129) when ap is taken as 27, 
The total lift obtained in this way, for constant a, must 
be a maximum as shown in reference 4. In other 
words, any approximate answer will be less than the 
exact theoretical lift. The variational equation has 
been solved by the Rayleigh-Ritz method using a two- 
term approximation for the lift, namely, 


| R2\2 R2\2 
cc, = 8 V(° — >) — (y — _ 4 
[ : 1) 4 44 (y — (R?/y)) 198 
A, — A3) A; 32 
: * (b — (R*/b)) 


where A; and A; are arbitrary constants to be deter- 
mined from the variational equation. If A; is taken as 
zero in Eq. (132), it may be observed that the cc; dis- 
tribution reduces to that of Lennertz and Spreiter, Ea, 
(S86). Calculations by this method have been made 
for a range of wings of various aspect ratio, taper ratio, 
and wing span to body diameter ratio in reference 4. 
Typical results are shown in Figs. 4 and 8. 

By manipulations on the integral equation, Eq. (129), 
similar to those employed to obtain Eq. (107), the fol- 


lowing result is obtained 


b 
/ I’; [a2(y) + QR, f(y) ] dy = 
b 
~h 
/ I's [ai(y) + ap, f(y) ] dy (1383) 
= 


where ay and az, produce circulation I», and a and 
ap, produce circulation T';. This is again in accord with 
the reversed-flow theorem except that, for lifting-line 
theory, flow reversal is not necessary. 

The variational principle is sometimes useful in an- 
other form, particularly if a rapid estimate of lift is re- 
quired. In this case, use is made of the fact that the 
expression 





lre 


typ 
sufl 
dov 


and 


func 
prov 
the: 
lar 


This 
Wei 
fact 
tion 
frst 
thos 
the 
finite 


I ipid 
ombi- 


ven in 


luded 
Eq. 
, this 


S 22. 
must 
»ther 
1 the 

has 
two- 


-ter- 
nas 
dis- 
Ea. 
lade 
tio, 
e 4. 


29), 


fol- 


ind 
ith 
ine 


in- 


the 
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*} . 


, or 
lwdy + / dy 

l b e b dol 
U ins / 
I ia(v) + apl(y ldy 
tT Pe y) 1 dys 


js stationary and is a minimum. Thus, if an assumed 
l distribution such as Lennertz’s and the corresponding 
@ distribution are substituted on the right-hand side 
the term on the left-hand side will differ from the exact 
value by terms of higher order than the error in the 
assumed lift distribution. For a(y) = ag = unity, the 
term on the left-hand side is equal to the reciprocal of 


pl) times the total lift on the wing plus body. 


Theories of the Weissinger Type 

One of the shortcomings of lifting-line theory men- 
tioned previously, namely, that the downwash at the 
3 4 chord point is actually greater than one-half the 
lrefftz plane value assumed in lifting-line theory, may 
be overcome by the use of theories of the Weissinger 
type. The other shortcoming arising from the lack of a 
sufficiently tractable form for the correction to the 
downwash due to the presence of the body w,, except 
at the lifting element and in the Trefftz plane still re- 
mains, however, and will necessarily limit the accuracy 
of any theory. 

The Weissinger theory can be constructed from the 
lifting-surface equation by making the approximation 


in Eq. (48) 


V (c2/4) + (y — 9)? (135) 
and then integrating the equation weighted by the 
function V (c/2 + x)/(c/2 — x) over the chord. The 
procedure is similar to that used to obtain lifting-line 
theory, and will not be repeated here. For a rectangu- 
lar wing alone, the result obtained is 


: r'(y) 1 f° (dV /dn) d 
Usa(y) = — + J ae ih 
: ™C j7 J—-b Vy 
l {¢ dn) dy 
47 b jee 
j 1(y—n)? 2\y — nh) 
(en 


c* Cc 


(136 


This corresponds to the well-known equation of the 
Weissinger ‘ZL Method.” It is equivalent to the satis- 
faction of the boundary condition at the 3 4 chord posi- 
tion with the lifting line at the 1/4 chord position. The 
frst two terms on the right-hand side are identical to 
those appearing in the Prandtl lifting-line theory, while 
the third term represents an additional correction for 


{nite aspect ratio. 
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The form of Eq. (136) as originally given by Weiss 
inger is obtained by combining the terms on the right 
hand side in a different way to give 


; | [ (dT dy dy l / 
l a(y) = + 4 
2r J-b y n tr J -! 


(dl dn) dn j +( n)~| a 
cag: Eee” gag Ee 


Here, the first term on the right-hand side is the down 
wash in the Trefftz plane and the second contains all 
corrections due to the finite length of the vortex sheet 

Returning now to consideration of wing-body com 
binations, it is clear that the main difficulties in the ap 
plication of this type of theory arise in connection with 
the last term of Eg. (136). All the other terms have 
already been treated in the discussion of lifting-line 
theory. The case of a circular body may be handled 
approximately by making use of the well-known 
centrally mounted wings, 


method of images. For 


this is particularly simple. In the use of this method, 
the lifting line is considered to pass through the body, 
but boundary conditions are satisfied only on the ex 
posed wing, since the images themselves effect the 
satisfaction of the boundary conditions on the body to 
some degree. In fact, the body boundary conditions 
are exactly satisfied at the chordwise location of the 
lifting line and in the Trefftz plane. The values of 
dV /dyn on the elements of lifting line inside the body are 
determined by the values on the appropriate point of 
the wing as 


r’(R?/y) = I’'(y y¥<K 138) 


With the image approximation, Eq. 
written for a wing-body combination as 


ry) 1 {°° (d¥/dn) d 
Uvfa(y) + azf(y)] = —— + / ville a 


1 tor. b 7 
WB " ” 


I " (d0 dn) j Hy — 9) 2 | 
J. ingylt 4 19 ~ wf 


rf 
im“ Rp Y— 7 


where the subscript IV’,B on the integrals over the span 
indicates that the range of integration includes the sur- 
face of the wing passing through the body in contradis- 
tinction to what has been done in most other cases in 
this paper. Eq. (139) is valid over the wing only, 
since I’’ in the body region is specified’ by Eq. (138) 
For wings which are on circular bodies but which are 
not centrally mounted, the method of images may still 
be applied, but the equations become more compli- 
cated, due to the fact that the image vortices do not lie 
on the extension of the plane of the wing through the 
body (the image of a vortex must lie on the line con- 
necting the vortex location with the center of the 
body). Incidentally, it should be noted that for circu- 
lar bodies the method of images would give results in 
slender-body theory and in lifting-line theory com- 
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pletely equivalent to those given above based on con- 
formal transformation. 

In the case of a straight wing mounted on a body of 
noncircular section, the approximate procedures which 
may be adopted are not so intuitively obvious as in 
but there are at least two 
The first approach 


the case of a circular body, 
possibilities which can be explored. 
would be to adapt Eq. (136) so that it conforms to lift- 
ing-line theory at the limit of high aspect ratio and to 
slender-configuration theory at the limit of low aspect 
This can be done by writing Eq. (156) to include 
(129) and 


ratio. 


the effects of the body (by analogy with Eqs. 


(y) + f(y) 
Y) 4 [1 (3 fol f 
[1 aly (dT cara 


¥-—%) 2). - 
V+ “3 7-H Aly (140) 


where the range of integration does not include the re- 
gion of the body and the bars indicate the same trans- 
The introduction 


(71)) as 


+ apl(y)] = 


Uvla(y) 


(dV /dn) d 
tf tf n ( n ihe 
7 =F 


formed coordinates as in Eq. (129). 
of the transformed coordinates ¥ and 7 into the last 
which represents three-dimen- 
however, lead 


term of this equation, 
sional flow effects, is arbitrary; 
to the correct slender-configuration result in the limit 
It is also easily seen that as c > 0, Eq. 
(129) for lifting-line theory. If 


it does, 


asc—> 
(140) approaches Eq. 
the form of the integral equation given by Weissinger, 
Eq. (137), is used, a possible ener is 


[1 + f(y) ’ (dV /dn) d 
i 


+ apf(y) } 
7 


U la(y) 
{1 + f(y) ] > (aT dn) dn § 5 = 9) 
ar wo 

(141) 


where the range of integration excludes the region of 
and the bars again transformed 
This equation approaches the equation 
but does not 
It has 


the body indicate 


coordinates. 
of slender-configuration theory asc > ©, 
approach the lifting-line equation as c > 0. 
the same interesting property as slender-body theory: 
for a = ag = const., the lift of the wing-body combina- 
tion is the same as the lift (as given by the Weissinger 
theory) of a wing alone having the same chord, and the 
transformed span. Whether Eq. (140) or Eq. (141) 
gives more satisfactory results in the practical design 
range of wing parameters is not known and must await 
detailed comparison with experiment. 

For swept wings, the direct application of the Weiss- 
inger procedure seems impossible, even where the image 
approximation is used, since the images of points near 
the tip may lie aft of the wing root trailing edge. For 
this case, Zlotnick and Robinson® have employed a sys- 
tem of horseshoe vortices and their images to represent 
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the lifting elements of the wing-body combination 


the boundary condition, Eq. (42), is then satisfied at 
the 3/4 chord position on the wing. 


trarily 


They also arbj- 
multiplied the downwash due to the bound 
+ R 


the presence of the body. 


vortices by (1 y*) in an attempt to account for 
This, as they have pointed 
out, amounts to assuming that the effect of the body 
on the flow field produced by the bound vortices is th 
as the effect of the body in modifying a uniform 
Zlotnick and 


Robinson® also make some reference to work in this field 


same 
two-dimensional downflow from infinity. 


in a paper by Zlotnick and Diederich which is as yet 
not available to the general public. 

For swept wings on bodies of noncircular cross sec 
tion, it appears necessary at present to resort to quite 
crude approximation in calculating the spanwise load 
distribution. It is obviously possible to construct an 


analogue to the horseshoe-vortex and image method 
of Zlotnick and Robinson by distributing the lift on the 
noncircular body arising from a given wing horseshoe 


vortex, as given by Eq. (21) laterally over the body 
(this can be done so the total lift is distributed exactly 
if desired, on the basis of a Trefftz-plane solution as dis- 
cussed previously). The body lift distribution could 
then be approximately represented by horseshoe vor- 
tices inside the body contour; the longitudinal location 
of the front edge of the horseshoe must, of course, always 
coincide with that of the wing horseshoe vortex which 
induces it. Since the labor involved in this procedure 
is considerable, it may in many cases be desirable to 
resort to even cruder approximations such as replacing 
the body by a plane wall at the side of the body, or 
replacing the body by the extension of the wing through 
the body set at zero angle of attack. In each of these 
two latter approaches, the Weissinger procedure in 
its original form could be used subject to the modified 


boundary condition on the wing 
a= aly) +- apf (y) 


The relative merits of these various methods can 
only be evaluated by extensive correlation of theory and 
experiments which are not yet available. For many 
practical purposes, it may well prove that the simplest 


approximations are adequate. 


The Moment Produced by the Body Lift in Subsonic Flow 

The Trefftz-plane theorem for lift may easily be aug- 
mented by a similar result for pitching moment. For 
infinite bodies of circular cross section, this result has 
already been given by Lennertz! and in extended form 
by Zlotnick and Robinson.’ Consider a lifting strip of 
width dx carrying a given spanwise distribution of lift 
p(x,v) dx. The potential due to this lifting strip will 
be g(x,y,z) and in the Trefftz plane it will be g~(y4,3). 
If the same lifting strip with the given spanwise dis- 
tribution is placed in reversed flow, the potential, de- 


noted by g, is 
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e = —[¢(x,y¥,2) — ¢g~(y,3z) ] (142) 


since the two flows are uniquely determined by the 
specified distribution of lift and differ only in the re- 
quirement that the potential vanish at x = — © for the 
direct flow and at x = +o for the reversed flow. The 


linearized pressure at a point on the body is given by 
p = pl o(0¢/Odx) 
in direct flow and by 
pb = —pl(0¢/0x) 


in reversed flow. Thus, for any point on the body p = 


p, and it may be concluded that the pressure distribu- 
tion on the body produced by the spanwise lifting strip 
is longitudinally symmetrical about the location of the 
strip. The element of body lift produced by the strip 
is, from the Trefftz-plane theorem, Eq. (21) 


dL = JS Ap(x,y) f(y) dy 


The moment of the body lift about x = 0 is then given 


un= f f xf(y) Ap(x,y) ds 
Ww 


where Ap is the pressure difference per unit area, and 
the integration is over the wing surface. 


by 


(143) 


Distribution of Lift Along an Infinite Body 


In order to apply the various theories discussed thus 
far to bodies which are finite in length and are non- 
cylindrical, it is expedient to examine first the longi- 
tudinal distribution of the lift induced by a wing lifting 
element on an infinite body. This will be done on the 
assumption that the body alone is slender, as was done 
by Multhopp? in his treatment of this problem. Fora 
cylindrical body, this assumption of slenderness really 
has nothing to do with the shape of the body, but is 
actually an assumption that the wing-induced upwash 
at the body varies slowly enough in stream direction so 
that slender-body theory may be This 
assumption will actually be violated near the wing for 
lifting elements in close lateral proximity to the body, 
but great accuracy cannot be obtained in this region 
Further, it will 


employed. 


without considerable complication. 
be seen that fair agreement with more accurate methods 
of analysis and with experiment can be obtained even 
in this region. 

Consider a cylindrical body in the presence of two 
symmetrically disposed lifting elements carrying lifts of 
amount pl\I'Ay at a lateral distance y, from the axis 
of the body. For convenience, they will also be as- 
sumed to lie in the x,y plane through the body axis 
parallel to the direction of flow. The upwash induced 
along the axis of the body may be obtained from the 
expressions given by von Karman and Burgers** for the 
induced-flow field of an airfoil of infinitesimal span as 
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rAy (x? + y?) (r — x) + 2°7 
— 


w= . . (144) 
4ar ri(r — x)* 
Similarly, the longitudinal velocity is 
u = (TAy/4m) (2/r*) (145) 


where the origin of coordinates is at the lifting element 


and 


r= Vx? 4+ y? + 2? (146) 


Taking account of the fact that there are two such lift- 
ing elements, the induced upwash at the axis of the body 
(z = O) is given by 


rAy x 
Ss . + 
<7y" Vai 


The total lift on a circular body due to the upwash ac- 
cording to slender-body theory is given by 


(147) 


Le, = / (xpR?Uww) dx (148) 
Thus, 
Le, = reR?Uow (149) 
since 
w = 0 
Finally, 
Le, = (pUo. Ay) : (150) 


y? 
Since there are two lifting elements carrying lift 
(pl I Ay), this is just one-half the amount of lift which 
is induced on the body according to Eq. (22). The 
other half must arise from the pressure field of the wing 
lifting elements acting on the body surface. This lift 


is given by 


Le, = / f plu cos (n,z) ds dx (151) 


where cos (7,2) is the z direction cosine of the outward 
normal to the body surface and ds is the element of 
The value of u is positive above 
If the first term 


body circumference. 
the plane z = 0 and negative below. 
of a Taylor's series expansion in 2 is used to represent 


u, Wis given by 
rAy Zz 
or (\ x + y2)3 


“u= (152) 


where account has been taken of both lifting elements. 


The expression for Ly, then becomes with z = R cos @ 
"© TR®Ay pU es 
La, = cos? 6 dé dx 
de (Ve + 2) Jo 
or 
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This then accounts exactly for the other half of the lift 
on the body. Its distribution along the body is given, 
tu the slender-body approximation, by 


dLz . FR*Ay | 7 
-_ = ply r 2 ( Led) 
dx 2 (Vet) 
It is now necessary to go back and find the longi- 


tudinal distribution of Lz,, the lift induced by the up- 
wash. This is given by 


dLp,/dx = mrpR?Uy (dw/dx) (155) 


dw/dx is obtained by differentiating Eq. (147), which 


gives 
dw TAy I 
dx 2 (Vx? + y?)' 
Thus, 
dLp, pl iT AyR? l -2 
= ; : (156) 
dx 2 (v x? + y?)' 


This is identical to the expression obtained for dLz,/dx, 
so that it may be concluded that not only are the total 
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amounts of both components of lift equal, but that 
they are identically distributed longitudinally. 

The significant thing about the longitudinal distrj- 
bution of body lift is that it diminishes very rapidly 
with increasing distance from the lifting element which 
produces it. As may be seen from Eqs. (154) and (156 
for large x it decreases as 1 /x*. 

Lennertz! has made calculations for the distribution 
of lift on a cylindrical body in the presence of a wing 
without restriction to slender-body theory. Lennertz’s 
results are shown in Fig. 5 for the body lift induced by 
a uniformly loaded lifting line with R 6 = 0.1 and by a 
lifting line carrying the lift distribution for minimum 


induced drag with R/b = 1/2. 


For comparison, the 
induced body lift for constant wing circulation accord- 
ing to Eqs. (154) and (156) is shown for R/b = 1,2 and 
U1. 


ent method are only in the longitudinal distribution: 


It is emphasized here that any errors in the pres- 


the total body lift is exact. The body lift is, according 


to Eqs. (154) and (156), for constant I, 


dls _ jurr 1% dy 
dx R (rv “2 + y2)3 
dlp — pUS j b R 
dx @/RLVxs 452 Vet + al 
(157) 
This form is indeterminate at x = 0, but it is not diffi- 


cult to show that the value of dL, /dx there is 


(dLp/dx),.9 = (pUol'/2) [1 — (R?*/d?)] 


Considering the elementary nature of the methods em- 
ployed here, the agreement between Lennertz’s values 
and those from Eq. (157) is surprisingly good. 


Distribution of Lift Along a Finite Body 


Applying the slender-body theory to determine the 
lift distribution on a body of revolution of finite length 
and varying radius gives 


dLy/dx = f pUuR cos 6 do + 
[(d/dx) (rpR?Uww) | 


(158 
If u is represented by (Ou /0z)oz, as before, this becomes 


dLy/dx = f pU(OU/dz)oR? cos? 6 dd + 


(d/dx) (rwpR?Upw) (159 


It may be noted here that because of the irrotationality 


Ow /Or. 
dLp/dx = rpR?U)(Ow/Ox) + (d/dx) (rpR?Uww 


conditions 0u/O0z = Thus, 


The expressions for w and Ow /Ox on the axis of the body 
due to any distribution of lift over the wing span may 
be obtained from Eq. (147). 

Comprehensive experimental data on the distribution 
of lift over prolate spheroidal bodies in combination 
with rectangular wings obtained by Moller have been 
presented by Schlichting** in his comprehensive mono- 
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graph on wing-body interference. These data have 
been compared by Schlichting with the results of an 
analysis by Multhopp’s method.* Multhopp’s method, 
which is mainly intended to serve as a means of pre 
dicting the moment about the wing | 4 chord point 
due to wing-body interference, is based on a number of 
approximations to the flow field around the body. 
First, it is assumed, as has been done here, that the 
load distribution on the body may be calculated from 
the slender-body theory of Munk.** Second, it 1s as 
sumed (although not stated in this way) that the total 
lift which would be produced on an infinite cylindrical 
body by the wing lifting line acts on the noninfinite 
body as well; for an infinite body this lift has no mo- 
t chord point and it is assumed that 


An attempt 


ment about the | 
this is the case for a finite body as well. 
is therefore made to exclude this part of the body force 
the This is 


first, the forces on the body pro 


distribution from moment calculations. 
done in two ways 
duced by the direct action of the pressures induced by 
the wing (dL, dx) are ignored, and second, the induced 
flow transverse to the body axis is approximated by a 
flow field which omits any contribution from (dL,, dx) 
in the immediate vicinity of the wing itself. This flow 
field is obtained by calculating the actual upwash on 
the axis of the body ahead of the wing from the wing 
span loading, assuming the flow is parallel to the wing 
along the wing chord (zero upwash if the wing chord 1s 
parallel to the body axis), and assuming that behind 
the wing the upwash changes linearly from a value of 
zero at the trailing edge to a value of a{l — (de da)| 
at the tail, where ¢ is the downwash angle at the tail 
calculated from the wing loading. The load distribu 
tion along the body resulting from this special induced 
flow field is then calculated from Munk’s slender 
body theory.* 

A comparison of the load distribution along the body 
computed according to the present method with experi 
ment and with the load distribution according to Mul 
thopp'’s method as given by Schlichting**® is shown in 
Fig. 6 for a mid-wing configuration having a rectangular 
wing of aspect ratio 5. In calculating the theoretical 
upwash for this case, the wing-lift distribution was ap 
proximated by one of constant average circulation. 
The same is usually done in the application of Mul- 
thopp’s method. The slender-body theory calculations 
by the present method were also corrected at the ends 
of the body by substituting sin 7 cos 7 for dR dx, where 
ris the angle whose tangent is dR/dx. This correction 
accounts for the inclination of forces due to normal 
pressures at the ends; it conforms to the exact solution 
for an ellipsoid at constant angle of attack. Also, the 
forces on the afterbody were computed with due allow- 
ance for the closing of the vortex sheet around the body 
according to the method of Adams and Sears’ which is 
discussed in the following section. 

It may be seen from Fig. 6 that the results obtained 


lrom the present theory are in fairly good agreement 


N 
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with experiment. In fairness to Multhopp’s method, 
it must be stated that it serves very well the purpose 
for which it was intended, namely, the prediction of the 
shift in neutral point of a wing-body combination 
caused by interference. Schlichting** shows that theo- 
retical predictions of the neutral point shift by Mul- 
thopp’s theory are in excellent agreement with experi- 
ment for a wide range of wing and body parameters. 
It is felt, however, that the present method is worthy 
of consideration because it possesses logical consistency 
which is lacking in Multhopp’s method and gives a bet- 
ter approximation to the actual load distribution along 
the body, while at the same time it requires but little 
more labor than Multhopp’s method. Moreover, the 
present method can be directly applied to swept wings 
merely by computing the upwash along the axis due to 
while Multhopp’s cannot be applied 


a swept wing, 


without further modification. 
Loading on Afterbodies 

In slender-configuration theory, the basic assumptions 
lead inevitably to the result that the lift induced on 
the body occurs at the same longitudinal location as the 
wing lifting element inducing it. Thus, if the after 
body of a slender wing-body combination is of constant 
cross section, no lift can be induced on it behind the 
wing trailing edge unless displacement and rolling up 
of the vortex sheet are taken into account. For an 
afterbody of revolution which closes to a point behind 
a slender wing-body combination having a straight 


wing trailing edge, Adams and Sears" have given an 
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expression for afterbody lift which takes into account 
the shifting inward of the wing trailing vortices as they 
follow the streamlines around an afterbody of decreas- 
ing cross-sectional area. Vertical displacements and 
rolling up of the vortex sheet are neglected, so there is 
an implied restriction to short bodies and light loadings. 
Their result for the position, ye, of a vortex originally at 
y, at the wing trailing edge is 


R.2 — R,;? (160) 


yo? — yy? = 


where R, is the radius of the body at 1 and R» is the 
radius at the location 2 at which the new position of the 
vortex is desired. Their general expression for the 
afterbody lift between points 1 and 2 was obtained 
from the momentum integral in the Trefftz plane, Eq. 
(18). (Note again that in slender-body theory, every 
plane is treated as a Trefftz plane insofar as the flow 
ahead of the plane is concerned.) In the present nota- 


tion, the afterbody lift as given by Adams and Sears is 


) R2 
= 2 r(iy) (1+ —)dy - 
R Vo° 
by . Ry? 
r(iy) (1+ —] dy (161) 
Ri Vi" 


This is 


ALi» 
ply 


where yo and }. are obtained from Eq. (160). 
simply the difference in momentum of the flow between 
constant circulation 


stations 2 and |. The case of 
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(which corresponds to rolled-up vortices) is of interest 
and can be easily dealt with. For constant circulation, 
Eq. (161) may be integrated to give 


ALy/pUy = 21} [bo — (R2?/by)] — 
[br — (Ri?/by)]f (162 
so that 
ALy = 2pUcT'[(be — bi) — (Re®/be) + (Ri2/b:)] (168 
and from Eq. (160 
bo = VR? + b? — RP ‘4 


If the body closes so that Ry, goes to zero, the afterbody 
lift becomes 


ALi = 2pl |S b,? = R,? — b, . 2 (R;? b,) | 165 
If 6, R is large, the radical may be expanded in series 
and only terms up to order (Rj/6;)* retained. Doing 
this leads to 

ALi = pu ol (R:? b;) 166 


which in the case of large b;/R amounts to a small posi- 
tive lift. For comparison, the lift on the body at station 


| is for constant circulation 
Lis — 2plT[R, —_ (R,? b) | 167 


The ratio of afterbody lift to body interference lift for 
constant circulation based on Eqs. (160) and (165) is 
shown in Fig. 7. It can be seen that the ratio grows 
large for R; b; approaching unity. 

Another source of lift on long afterbodies is the dis- 
placement of the wing vortex sheet upward from the 
body as the flow moves downstream. The problem 
has been studied in detail by Lagerstrom and Graham" 
on the basis of slender-body theory. Moskowitz! 
has estimated the total loss of lift behind 
wing-body configuration with cylindrical afterbody on 
In this case, 


a slender 


the assumption of a very long afterbody. 
the wing trailing vortices are eventually displaced to a 
great enough distance from the body that the wing- 
induced upwash at the body base becomes effectively 
zero. Rather sizable losses of lift on the afterbody were 
predicted in contrast to the results of Lagerstrom and 
Graham. In fact, it can be shown that Moskowitz’s 
calculations give an afterbody negative lift which is 
exactly equal to all the positive body interference lift 
induced ahead of the wing trailing edge. It appears 
that these results are in error, since only the effect of 
change in upwash at the body, dLz dx, was taken into 
account in Moskowitz’s calculations. Actually, if the 
trailing vortices are deflected from a position parallel 
to the body axis, a longitudinal component of lift 1s 
induced directly by pressures acting on the body. As 
was shown by Lagerstrom and Graham, the momentum 
equations may be applied to this case with simple re- 
sults if rolled-up concentrated vortices are considered. 
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Far downstream, the rate of flow of momentum (ig- 
noring the uniform flow past the cylinder) is simply the 
momentum of a vortex pair plI'2b. where 2. is the 
distance between the two rolled-up vortices. The rate 
of flow of momentum just behind the trailing edge of 
the wing is the same as calculated for station | in Eq. 
162 Thus, the lift on the afterbody, which is the 


difference of these two values 1s 


Li = 2pUc jbo — [br — (Ri*/d)); (168 


Lie = 2pUT [(b2 — b;) + (R17/b (169) 
If interactions between the vortices are neglected, it is 
easv to show from the two-dimensional flow around a 


circle that 
bo = db, — (R,7/d; (170 


Substituting this into Eq. (169), it is seen that to this 
approximation the afterbody lift is zero. 


and Graham have shown that a small] positive lift does 


Lagerstrom 


occur if the interactions of the vortices are consid- 


ered. 


Simplified Method for Estimating Lift of a Wing-Body 
Combination 
Ward’ has suggested that the result obtained from 
slender-configuration theory for the lift of a winged body 
of revolution in terms of the lift, L*, which the exposed 
wing would experience if mounted on a plane wall may 
be applicable to more general wing plan forms than the 


derivation would imply. His result is 


Lw, = Lw*{1 + (R/d)]? (171) 
A similar result is suggested by calculations according 
to lifting-line theory‘ in Fig. 8. It is interesting to note 
that whereas the lift of the wing-body combination in 
terms of the lift of the wing including its extension 
through the body shows a marked dependence on as- 
pect ratio, as shown in Fig. 4, aspect ratio dependence is 
largely suppressed when the wing-body characteristics 
ire referred to those of the exposed wing in Fig. 8. 
Also plotted in Fig. + are Ward's results and several ex- 
perimental results** * * for wings of both high and 
low aspect ratio. The agreement between the two 
theories and experiment is fair, considering the wide 
range of configurations covered. For the experi- 
mental cases, the lift of the exposed wings alone was in 
all cases estimated from the Weissinger theory, since 
whenever wing tests are made to supplement wing-body 
tests, current experimental practice is to test not the 
exposed wing, but a hypothetical wing, including the 
extension of the exposed wing through the body. 

From the Trefftz-plane relation, Eq. (21), it is ap- 
parent how the method of approximation can be ex- 
tended to arbitrary body shape. The average angle of 
attack experienced by the wing is given by 
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where /(y) is the same body upwash function introduced 
previously and a is the distance from the plane of sym- 
metry to the wing-body intersection. If C),,.* is the 
slope of the lift curve for the exposed wing, then the 
slope of the lift curve of the wing in the presence of the 


body will be approximately 


& 1+ oa)C 173) 


ay La 


The lift induced on the body can be estimated from the 
Trefftz-plane relation, Eq. (21); assuming the lift on the 


wing is distributed as |1 + /(v)] 


/ 1 + f(y) [fly) dy 


The slope of the lift curve for the wing-body combina- 
tion is then 


. (175) 


au 


where C),,,, 1s based on the exposed wing area. If, in 
the calculation of the induced lift on the body, the 
wing-lift distribution is assumed constant, 4 = o, and 


the equation for C),,,,, becomes 


(176) 








314 JOURNAL OF THE 


For a body of revolution, this would lead exactly to 
Ward's equation; for most cases, however, Eq. (175) 
seems to give slightly better results. 
improving the present method at the expense of addi- 
tional labor are obvious; better approximations for the 
wing span-load distribution in the presence of body up- 
wash may be made, the induced body lift could be 
simply calculated from Eq. (21) once this improved 
loading were available, and finally, corrections for 
finite-body length could be estimated where necessary. 


SUPERSONIC W1ING-Bopy INTERFERENCE 


The aerodynamics of supersonic aircraft have been 
subjected to intensive study for a relatively short time, 
but already there has been a great deal of work on the 
wing-body interference problem. The various pro- 
cedures which have been employed in attempts to ob- 
tain exact linearized solutions of this problem will first 
be reviewed in historical order. This will then be fol- 
lowed by a discussion of the practical applications in the 
light of experimental results. Slender-body theory, 
which is applicable to sufficiently slender configurations 
in supersonic flow, will not be discussed further, since 
the discussion already given in connection with the sub- 


sonic problem is directly applicable. 


Ferrari Method 

Ferrari'* '* considered the problem of a flat rectangu- 
lar wing attached to a body of revolution in supersonic 
flow where the wing and body are at a common angle of 
attack. He employed an iterative process in his solu- 
tion; the basis for this process has already been given 
under ‘‘General Theory and Concepts.’ The same no- 
tation will be employed here. If the wing is assumed 
to be extended through the body, the projection on the 
z = O plane of the region common to both will be de- 
noted by Rws. The four-step iteration procedure em- 
ployed by Ferrari may be summarized as follows: 


Step 1: 
Ogw™/dz = — Uva on Rw + Rwez 
Step 2: 
Ogp” /On = — Uva cos 6 on Rpg 
Step 3: (177) 
Ogw'? /Oz = —OgR/OZ on Ry 
Oow )/dz = 0 on Ry B 
Step 4: 
Ogz'/On = —Ogp/On on Rp 


The first two steps of the iteration procedure may be 
carried out by making use of the Ackeret*® solution for 
a two-dimensional wing and the Tsien®® procedure for 
a body at angle of attack. In computing the third 
step, Ferrari makes use of a Fourier expansion of the 
spanwise lift distribution on the wing which may be 
regarded as essentially a representation of a doublet 


distribution. 
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In the fourth step, the velocity potential of the flow 
around the body is taken to be of the form 


_ m _ 
g= 7. On(x,r)r” cos m6 78) 


in polar coordinates with r = 0 on the body axis, where 


the ¢,,'s are given by 


, *O 
> {| O n 
Yn =? ( fm(x — ar cosh uw) du 
ry Or cosh v/a 


9) 


where the /,,’s may be regarded as multipole distribu- 
tions along the axis and a = V JJ? — 1. Cramer“ has 
shown that a form of the solution much more conven- 


ient for computation is 


mom *0 
(—1)"a 
On = : pS 
r v/a 


- cosh 


(x — ar cosh wu) cosh mu du (180 


The f,,’s may be found numerically by the methods 
originally given by von Karman and Moore for /).* 
From the form of Eq. (178), it is clear that only ¢, 
contributes to the lift distribution along the body, 
since the linearized pressures are proportional to 0g Ox 
and these are multiplied by cos @ to obtain the lift. 


The Browne-Friedman-Hodes Conical Solution 


A conical wing-body configuration is defined as one 
which can be generated by a line rotating about a fixed 


9 


point. Browne, Friedman, and Hodes”? have applied 
the powerful methods of conical flow and conformal 
transformation to this problem. Analytic solutions for 
delta wings with both subsonic and supersonic leading 
edges were obtained in this fashion. Although the 
conical configuration is of small practical importance, 
this work contains the only solution of a wing-body 
interference problem with subsonic leading edges for 
nonslender wings. The results of reference 22 have 
beea used to estimate the regime of applicability of the 
slender-configuration theory of Spreiter’ and Ward.° It 
is clear that the conical configuration will be of great 
interest as a check of the more general theories for inter- 
ference between delta wings and bodies which lie in the 


future. 


The Morikawa Procedure 


Morikawa” has considered the problem of a rectangu- 
lar wing attached to a cylindrical body of circular cross 
section at a common angle of attack « in a supersonic 
flow field. Upon specializing the results of the general 
formulation previously given to an infinite circular 
cylinder at a Mach Number of 7/2, the problem may 
be formulated in terms of a potential gy as follows: 


O-gw/Ox? = (O° gw Oy") + 
(O° gw /Oz”) in flow field F (181) 


Ogw/Or = 0 onr = 1 (182) 
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| onx> tt: r> 1 
(1S3) 


(1/7 Ogw/08) = — Uv [1 + (l/r- 


Use of the particular Mach Number, v 2, in the formu- 
lation of the problem does not restrict the generality of 
the result, since the Prandtl-Glauert transformation 
may be used to obtain results for an arbitrary Mach 
Number. 

The Laplace transform is applied to the differential 
equation and the boundary conditions. With the defi- 


nition 


eo 


V(7,6,5) = Lew (7,00) = / dse-* gw(7,0,x) 
0 


(184) 


the boundary-value problem is reduced to the form 


sy = 0*y/dy" + O0*y/02? 
= O°y/Or? + |(1/r) (Ow/Or)] + 


[(1/r?) (0°y/00?)] (185) 

oy/or = 0 onr = 1 (186) 
(1/r) (Ow /00) = —|[Up(a/s)] [1 + C/r*?)] on r = 1 
¢é=0 

(187) 


(1/r) (Ow/00) = [Uc(a/s)] [1 + (i/r*)J onr = 1 
6= xf 


(188) 


The solution obtained by Morikawa for the boundary- 
value problem in ¥ may be written explicitly in the fol- 


lowing form 


Uo ‘ ’ 
V(r,0,s) = = | dp (1 + ,) \Kolse: ) + 
TS l p- 


4U be l T2n'(s) 
Ko(sp-)} — = | dp\{1+ ) > &st 
7s l p~ n=0 Ko,'(s) 
Ko, (sr) Kon(sp) cos 2n@ (189) 


where J>,, Ko,, are the modified Bessel Functions, e, is 
equal to one-half when n = 0 and equal to | otherwise, 
and p+, p— are defined as follows: 


p+ = Vr? + 2rp cos 6 + p*| 


(190) 
= Vr? — 2rp cos é + p?s 


The potential gw may be computed from y by the in- 
verse Laplace transform 
a 
gw(r,0,x) = L7! Yi7,6,s) = : | dse™* W(r,6,5) 
<7 ys 
(191) 
where y is any real number greater than zero. Thus, 
the wing-body problem is reduced to two integrations 
and a summation and may be considered solved in prin- 
ciple. In practice, however, the analytical difficulties 
are so great that almost no progress has been made in 
extracting numerical results from the Morikawa formu- 


lation in its exact form. Morikawa, however, has ap 
proximated the series in Eq. (189) by a simpler function 
which he constructed semiempirically so as to fit condi- 
tions very near the leading edge. The range of va- 
lidity of this approximation has not been determined 


mathematically. 


Nielsen-Matteson Method 


Nielsen and Matteson” considered the problem of an 
untapered wing of double-wedge profile and 60 


attached to a cylindrical body of circular cross section 


sweep 


in supersonic flow with wing and body at zero angle of 
attack. 
inscribed in the cylinder. 
used to cancel the normal velocities induced on the 
It is clear 


The body is replaced by a set of plane areas 
A collocation procedure is 


body by the disturbances due to the wing. 
that the sources distributed over the body in this pro- 
cedure will cause normal velocities on the wing. Thus, 
this procedure may be looked upon as the first steps of 
an interative process. It is evident that the same 
procedure may be used to treat lifting cases if all edges 
are supersonic. The difference between the Ferrari 
and Nielsen-Matteson methods appears to be solely 
in the method of solving the body-alone problem which 


arises in the analysis. 
Nielsen Procedure 


The wing-body problem for a rectangular wing and 
infinite cylindrical body of circular cross section in 
supersonic flow has been treated by Nielsen for the 
case where the body is at zero angle of attack. Nielsen 
employed the special formulation of the problem sug- 
gested by Lagerstrom and Van Dyke”! in which a basic 
wing passing through the body gives one component of 
the solution. The flow set up by this wing violates the 
boundary condition of zero velocity normal to the body 
boundary. A second component solution must then 
be added to the first to cancel the flow across the body. 
This, thus far, is the same as the Ferrari problem. At 
this point, however, the Lagerstrom-Van Dyke formu- 
lation requires that flow normal to the body surface be 
cancelled without at the same time producing any new 
velocities normal to the wing surface. This can be 
done, for supersonic leading edges, by representing the 
solution by a Fourier cosine series in the meridional 
angle (6 = 0,7 on the wing) with even terms only. 
Thus, in effect, the iteration procedure is terminated 
in two steps. 

The potential ¢ is expressed in the form 


¢ = ow’ + ¢p (192) 

where gy’ satisfies the boundary condition 
Ogyw’/0z = —Uoa Rw + Rwe (193) 

and ¢ » satisfies the boundary conditions 
Ogn/Or = f(x,0) — (Ogy®/Or) onr = 1, x >O (194) 
Ovp/00 = 0 on@é=0,7, x>0 (195) 
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and the differential equation 
O°¢p/0x? = (0°pp/Ov") + (0792/02?) 


The wing is assumed to extend through the body at 
angle of attack a so that the Ackeret solution may be 
used to obtain gw®. The boundary-value problem 


is solved by taking the Laplace transform of the differ- 


ential equation and the boundary conditions. With 
the definitions 
V(r,0.s) = / dxe” *“pp(r,0,x) = Lop (196) 
Q 
F(0,s) = ii dx~** {(6,*%) = Lf(0,x) (197) 
0 
the transformed boundary-value problem may be 
written as follows 
sty = (O°y/or?) + [(1/r) (Ow/Or)] + 
[(1/r?) (02/007) ] (198) 
Oy /or = F(6,s) onr = 1 (199) 
oy/o9= 0 6=0,7 (200) 


The solution to the above boundary-value problem 
may be written in the form 


e,F’,(s) cos 2n0K.,,(s,7) 


¥(r,0,s) = > — (201) 
n=0,1 Ss Ko,'(s) 
where the following definitions are used 
ve ° x 
Pils) = | drF(7,s) cos 2nr (202) 
T JO 
I 9 
&=s5 n= 0O (203) 
én = 1 n>O 


It will be noted that the F,,(s) are the Fourier cosine 
series coefficients of the function F(@,s). The potential 
¢ is found by taking the inverse transform of the above 


expression for y or 


Ko, (sr) F,(s) 


=n 


¢n(7,0,x) = dle, cos 2ndLe-! (204) 
n=0 


sKo,‘(s) 


Nielsen has evaluated the inverse Laplace transform 
appearing in the above formula in terms of known 
functions and a new function .VW/s,(x) defined below 


Mo,(x) = £7! [1/se’Ke,'(s) ] (205) 
The following formulas may be obtained by elemen- 
tary applications of tables of Laplace transforms (say 


Van der Pol and Bremmer) 
Ko, (sr) F,,(s) 2\** = a(2n) ee 
} 1 7 = ° d 
§ tha, (5) r a(4n)Jo 


Qn—} 
- is + 1 — §* — ry” 


Str 

vS 
to 
Str 


—— (206) 
of" 


where g»,,(&) is defined by the formula 
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Pon(x) = J dé fon’ (E)Mon(x — &) 207) 
0 


and r(m) = n! 
Nielsen has tabulated the Mzo,(x) for 
n = 0,1,2,3 and x from 0 to 4 with intervals of approxi- 


functions 


mately 0.2. 

It should be noted that since the solution of the 
second component problem gives no flow across the 
wing plane, it is not valid behind the trailing edge of 
the wing. An entirely new solution must be worked 
out for the body behind the trailing edge based on the 
downwash field of the wing and body ahead. In the 
Nielsen solution, by contrast with the Ferrari solution, 
all terms in the Fourier series for potential contribute 
to the body lift, not just the first term. 


THE STATUS OF SUPERSONIC WING-Bopby INTERFERENCE 
THEORY—COMPARISON WITH EXPERIMENT 


Two essentially different procedures have been em- 
ployed in calculating supersonic wing-body interference 
the method of iteration 

The iter- 
Ferrari! '% 


for nonslender configurations: 
and the method of separation of variables. 
ation procedure, which is employed by 
and Nielsen and Matteson,” reduces the wing-body 
interaction problem to a set of general wing-alone and 
body-alone problems. The wing-alone and body-alone 
problems are handled exactly in a mathematical sense 
by Ferrari, and Nielsen and Matteson. In general, 
the iteration procedure is not entirely satisfactory for 
two reasons. Although the first stages of the iteration 
process may be computed rapidly, the numerical calcu- 
lations of the later stages become very tedious. Fur- 
thermore, there is some question as to the rapidity of 
convergence in regions which are strongly affected by 
wing-body interferences. In spite of these deficiencies, 
the iteration procedure is the only method which can 
be used in its present form to compute the character- 
istics of general wing-body configurations with subsonic 
leading edges. It will be noted, however, that such a 
calculation has not been carried out and appears to re- 
quire considerable labor. 

The method of separation of variables has been em- 
ployed by Browne, Friedman and Hodes” for the cont- 
cal case and Morikawa”? and Nielsen**’ for configura- 
tions with supersonic wing edges. Separation of the 
wave equation in conical coordinates is used by Browne, 
et al., to obtain a relatively simple solution of the coni- 
cal wing-body problem. Although a conical wing-body 
configuration is rarely employed in the design of air- 
frames, the solution has been used to establish approxi- 
mately the regime of applicability of slender-configu- 
ration theory and may be used as a test case in evalu- 
ating future theories. 

The procedures employed by Morikawa” and Niel- 
sen** are quite similar. By writing the potential as a 
sum of two terms, Morikawa reduces the wing-body 
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WING-BODY 


problem to a problem where the normal velocity is 
Nielsen 


makes use of a similar device to reduce the configuration 


specified on the wing, and is zero on the body. 


problem to one where the normal velocity is specified 
on the body and is zero on the wing. The Laplace 
transform is employed in both theories to reduce the 
problem to a quasi two-dimensional boundary-value 
problem. Morikawa determines the Green's function 
appropriate to the problem as a series expansion. The 
solution to the problem is then expressed as the inverse 
Laplace transform of an integral over the wing alone. 
Morikawa does not evaluate the inverse Laplace 
transform of the Green’s function, which would re 
quire the calculation of a hitherto untabulated func 
tion .\J/s,(1 


wing-body problem as the inverse Laplace transform of 
The Green's function 


Nielsen expresses the solution of the 


an integral over the body alone. 
derived by Morikawa may be used to obtain the same 
result. Since the Green’s functions for the Morikawa 
and Nielsen procedures are essentially the same, the 
same untabulated inverse Laplace transform which arose 
in the Morikawa theory appears in his result. Nielsen 
has carried out a numerical evaluation of the function 
Mo, (x) for n = 0,1,2,3. It is clear that four terms of 
the series obtained in the Nielsen and Morikawa 
theories may be evaluated from the tabulation of .\/2,. 
Although both theories are applicable to all configura- 
tions with supersonic leading edges, they have been 
applied only to a rectangular wing. Morikawa’s 
treatment of this problem is based on a semiempirical 


approximation to the Green’s function. 


In general, the method of separation of variables 
leads to a solution in series form, so that the problem 
of rapidity of convergence arises as in the iteration 
process. The accuracy of a four-term series turns out 
to be poor near the leading edge of a wing-body com- 
bination. In fact, Nielsen found it necessary to supple 
ment his series solution in the region from the leading 
edge up to a distance 1.2 V 1/° — 1 R downstream of 
the leading edge (where R is the body radius) by an 
approximate solution based on a quasi-plane-wall argu- 
ment. The distance 1.2 V M? — 
to the wing chord for many configurations, even for 
MJ = 14. The numerical labor required to obtain even 
one term of the series is comparable with the first stage 


| R is comparable 


of the iteration procedure. 

It is not possible at present to draw definite con- 
clusions as to the best methods for computing super- 
sonic wing-body interference. However, it does appear 
that, for wing-body combinations of large wing chord 
compared to the body radius, which are not followed 
by afterbodies (tail surfaces), the Nielsen procedure 
would offer advantages. On the other hand, for con- 
figurations having shorter chords and followed by long 
afterbodies (the more usual wing case) the Ferrari 
method offers the maximum amount of useful informa- 


tion for a given amount of labor. 
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A study of the available experimental data has shown 
that, while available theories fall far short of perfection 
in their present form, they do provide sufficiently ac 
curate results for most design purposes when applied 
to cases close to those for which they were derived 
Fig. 9, based on experimental data from reference 50, 
shows the pressures measured on the flat upper surface 
of a rectangular wing centrally mounted on a body of 
revolution which was essentially of constant cross sec 
tion downstream of the wing leading edg« The tests 
were made with wing and body both at a common angle 
of attack of 4 
are all referred to the theoretical predictions ot 
209 his 


For reference purposes, the pressures 
a for 
mula given by Kirkby and Robinson, Eq 
formula represents the exact solution according to 
linearized wing theory for the pressures in the region 
outside the root and tip Mach lines when the wing ts 
operating in the upwash field of an infinite cvlindet 
Beskin’s approximation For points theoretically 
inside the root Mach line, Kirkby and 
values on the Mach line were used as reference values 
It can be seen that Kirkby and Robinson's result 1s 1n 
very good agreement with experiment outside the root 
Mach line, being no more than 5 per cent high at the 


he »binson S 


worst. 
ories of Morikawa and Ferrari (first iteration 


The pressure distributions predicted by the the- 
are also 
shown in Fig. 9 for comparison. Morikawa’s results 
for the region outside the root Mach line (which, of 
course, is not influenced by the body except by upwash 
Ferrari's re 
When 


the region inside the root Mach line is examined, how 


are identical to Kirkby and Robinson's 
sults for the same region are somewhat high. 
ever, the situation is reversed. Ferrari's result shows 
fairly good agreement and Morikawa’s is too high. It 
should be emphasized that the test of Morikawa’s 
semiempirical approximation occurs only inside the root 
Mach line. His solution outside the Mach line was 
obtained separately in a manner similar to Kirkby and 
Robinson, and has nothing to do with the solution in- 
side. Ferrari's solutions inside and outside the root 
Mach line are, however, related since he used a single 
Fourier series to represent the potential over the span. 

The experimental load distributions along the body 
from reference 52 for a high-aspect-ratio wing mounted 
on a long cylindrical body are shown in Fig. 10. The 
body was long enough so that the wing-tip Mach cones 
did not strike the body until very near the base of the 
body. The cases of common angle of attack of the 
wing and body and wing incidence only are included. 
For the wing incidence study, one test was made with 
the opening at the wing-body juncture created by the 
wing angle unsealed, and the other with it sealed. The 
theory of Ferrari is the only one available for afterbody 
loading. The results of the first iteration according 
to this theory are shown in Fig. 
In the first step of the iteration process, 


10 for comparison with 
experiment. 
as presently constituted, there is no difference in the 
computed body loading between the incidence and 
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angle-of-attack cases. This could be improved, but 
only at the expense of additional labor. 
that the theoretical loading lies between the measured 
loadings for the wing incidence and angle-of-attack 
cases, being closest to the case of wing incidence with 


It can be seen 


sealed gap. 

No comparison of theory with experiment for the 
wing pressures in the wing-incidence case was possible 
because the only wing pressure measurements in such 
cases were made with large gaps at the wing-body junc- 


ture. 


Approximate Methods in Supersonic Wing-Body 
Interference 
The theoretical 

supersonic wing-body interference calculations become 


methods which are available for 
extremely laborious and even intractable when applied 
to many of the practical problems which arise in de- 
sign. Therefore, less exact and more simple methods 
for analysis of wing-body combinations have had to be 
devised. Of course, the results of more exact analyses, 
particularly if augmented by adequate experiments, 
serve as major guideposts in any attempt to find ap- 
proximate methods which can be relied on, at least for 
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Fig. 9 
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DISTANCE FROM ROOT 


some types of configurations. It must be kept in mind 
that higher-order effects are always present to modify 
at least some of the results which may be obtained even 
from the most exact linear theory. For instance, if one 
has to determine with great accuracy the hinge mo- 
ments on a movable control surface which is influenced 
by the body, it will do him little good to seek to deter- 
mine the pressure distribution from linearized theory 
with great accuracy. If the effects of gaps and viscos- 
ity do not influence the results, second order potential- 
flow effects most certainly will, unless control-surface 
thickness and angle of deflection are held down to im- 
practical limits. On the other hand, such things as the 
overall lift and pitching moment of reasonably conven- 
tional configurations can probably be predicted with 
satisfactory accuracy even by approximate solutions 
of the linearized equations. At present, only for the 
case of wings centrally mounted on a circular body of 
revolution are there sufficient theoretical and experi- 
mental data to permit any sort of general conclusions to 
be drawn. 

The earliest attempt at an approximate theory of 
supersonic wing-body interference was apparently that 
of Kirkby and Robinson.” They considered the case 
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WING-BODY I 


of a rectangular wing and a conical body at a common 
angle of attack. 


body upwash around the conical body and to compute 


Their procedure was to evaluate the 


the wing lift on the basis of two-dimensional strip theory 
Ackeret) *” 
chordwise upwash equal to the mid-chord value at each 
No attempt was made 
Their 


on the additional assumption of constant 


spanwise location on the wing. 
to evaluate the effect of the wing on the body. 
assumption of the applicability of strip theory on the 
wing would make the method applicable to wings of 
chord small compared to R V J/* — 1. 
of investigating the validity of the strip assumption, 


In the course 


Kirkby and Robinson have given in an appendix the 


simple useful solution, exact within the linearized 
theory, for a wing of infinite span subjected to the up- 
wash field of an infinite cylinder (in other words, 
Beskin'’s approximation®! when applied to finite-length 
bodies). This is the same body-induced angle of attack 


used by Multhopp? for finite lergth bodies, namely, 


a, = ap (R? y*) (208) 
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where az is the body angle of attack Their solution, 


for the difference between upper and lower surface 
pressure coefficients, is 
; j an Reapy( VM | ' 
AC, = 4 + 209 
VY/ie — 1 M? — 1 “a he 


where ay is the wing angle of attack, ag is the body angle 
of attack, and x is measured from the wing leading edge, 
while y is measured spanwise from the body axis. As 
may be seen from Fig. 9, this result is in good agreement 
with experiment in its region of validity 

Returning to Kirkby and Robinson's basic method, 
their use of strip theory would result in a wing-lift co 
efficient, based on exposed area, 


la R 
fi (1 4 4 
VM? - 1 D 


This neglects wing-tip effects 
for the 


210 


where / is the semispan. 
which may be corrected separately, at least 


basic wing lift. On the other hand, Cramer’s ‘Rules 
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of Thumb’’” based on experiment lead to the result 
that, for wings of large span followed by a sufficiently 
long afterbody, the total lift on the wing-body is 


: La 9 R 91 
Cr B =z a ] + ~ } (Z1 1) 
VM? — 1 4 


The lift due to the body alone must, of course, be added 
to this result. 
divided between wing and body, but is somewhat 
Thus, Kirkby and Robinson's 


The interference effect is not equally 


greater on the body. 
formula is more appropriate for the case of no after- 
body. In their theory, they take no account of load- 
ing on the body up to the wing trailing edge, but they 
also do not account for the fall-off of hit in the wing- 
body juncture. As Cramer has shown, these effects 
tend to compensate, so that Eq. (210) for total lift 
is probably not too bad. Qualitative theoretical 
arguments for the validity of this compensation have 
been advanced by Lagerstrom and Van Dyke?! on the 
basis of analogies to planar lift preservation theorems. 
As will be shown later, the reverse-flow theorem for 
wing-body combinations provides a more rigorous basis 


for this compensation. 
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If the wing chord is very small compared to the body 
radius and the wing span is large compared to the body 
radius, then strip theory should more nearly apply on 
the wing, since the region affected by the body would 
be a small part of the wing area. In this case, the lift 
coefficient along the span would be closely given by 


ta i? 

SS eer 

VM? - 1 oe 
The wing-lift coefficient would then be given by Eq. 


(210). From the Trefftz-plane theorem, Eq. (22), the 
lift of a long body in the presence of the wing is given by 


Le Lan 7 R?\R? 
= 2 (: + —~ |— ay (233) 
gq VM2?—1 JR y? Jy 


212 


Cy) = 


Integrating this, and finding the lift coefficient based 


on exposed wing area, gives 


ha R E 1R 1/R 
+ +r ( ) (214) 
/M2?—150L3 30 d \¢ 
Adding this to Eq. (210) gives 
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This can be regarded as the limiting value of the lift 
obtainable from a rectangular wing-body combination. 
It is undoubtedly too high for most practical cases, 
and something between Eqs. (210) and (211) will be 
more appropriate. 

The prediction of the pitching moments of wing-body 
combinations is a somewhat more difficult problem 
than the determination of the lift, and thus, the ap- 
proximate methods are not quite so simple. The calcu- 
lation of the first iteration for the body loading in the 
Ferrari method in which the wing-tip Mach cones are 
assumed not to strike the body, consists of two parts: 
first, the computation of the loading due to the pres- 
sure field of the wing acting directly on the body, and 
second, the loading due to the flow set up around the 
body to cancel the velocities produced by the wing. 
On a sufficiently Jong body, only the first part of the 
loading produces lift, since the second part of the load- 
ing is due to equal positive and negative angle-of- 
attack changes at the body due to the leading and 
trailing edge Mach waves. The second loading pro- 


duces a couple, however. By going to the limiting 
case of vanishing body radius, Ferrari has computed 


the moment of the couple as 


(216) 


AM, = mpl'y*R°ca 


where the minus sign denotes a nose-down moment, « 
The lift 


acting on the body due to the pressure field of the wing 


is the wing chord, and ap the angle of attack. 


in this approximation is simply equal to the lift which 
would be carried by the region of the wing going through 
the body. By a geometrical calculation of the centroid 
of the surface area intercepted by the Mach planes 
from the leading and trailing edges, the location of this 
lift is found to be (7 /4)RV AJ? — 1 from the wing mid- 
chord. The moment of this lift about the wing mid- 
chord is therefore 
tacR T pu,’ 
Al, = — RV M? - 1 (217) 
. VM? —-14 : 

The total moment about the wing mid-chord is the 
sum of AJ, and AM». 

Schrenk” has reviewed the arguments leading to 
these expressions starting from slender-body theory, 
rather than from linearized supersonic theory. In 
this way, the extension to noncircular bodies and other 
interference problems is made more obvious. 

As the wing aspect ratio becomes smaller, the action 
of the wing-tip Mach cones on the body alters the 
simple considerations given above. The lift on the 
body can no longer be taken to be the lift which would 
occur on the extension of the wing passing through the 
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body. Fig. 11 shows the loadings produced on the 
body for a set of wings ranging in exposed aspect ratio 
from 0.5 to 4.29. The body loading curves all start in 
substantially the same way, but at the location where 
the tip effect is felt for each particular wing, the body 
loading begins to fall away from the AR 1.29 case. 
It is obviously necessary to modify the Ferrari pro 
cedure to account for the wing-tip flow field in these 
cases. The ratio of lift on the body to lift on the wing 
for the series of wings tested is shown in the inset in 
Fig. 11 and compared with calculated values obtained 
by application of the Trefftz-plane theorem, Eq. (22), 
on the assumption that the span lift distribution on the 


wings is given by Eq. (212). The result is 


Ce, Ryd - R/b) + 1/3[1 — (R/b)*)]) 918 
cc 4 (1 — (R/b)?] ‘Bas 


The agreement between the theoretical and experi 
mental values of this ratio is seen to be quite good, 
although it is well known that the correct lift distribu- 
tion is not that of Eq. (212). This is another indication 
of the insensitivity of this ratio to the loading which 
was demonstrated earlier for the subsonic case. 

The prediction of the moment of the body loading 
arising with low-aspect-ratio wings requires more care. 
It appears that the best approach would be to use the 
Ferrari method for the body loading with a wing of in- 
finite span, and to correct this loading by an increment 
applied everywhere downstream of the point where the 
wing-tip Mach cone intersects the body. This incre- 
ment might well be calculated on the basis of slender- 
body theory using the local angle of attack on the axis 
of the body computed from one or another of the sim- 
plified theories of wing downwash. To facilitate com- 
putation of the infinite wing-body loading, Rae* has 
computed a normalized loading curve for a semi-infinite 
plane wing centrally mounted on an infinite cylindrical 
body which is given in Fig. 12. By adding this loading 
imposed at the wing leading edge to the negative of the 
same loading imposed at the trailing edge, the solution 
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for a wing of finite chord may be obtained as illustrated 
in Fig. 12. 

Since the ratio of lift on the body to lift on the wing 
has been shown to correlate so well with theory, even 
for low-aspect-ratio configurations, it would be a simple 
matter to compute the lift of a wing-body combination 
if an acceptable method for computing the lift on the 
wings could be devised. Fig. 14, which is based on 
reference 25, shows the spanwise loading on the various 
wings corresponding to the body loadings of Fig. 11, 
with the two-dimensional strip theory values for the 
wing without upwash, and for the wing subjected to 
the Beskin upwash also given for comparative pur- 
poses. It is clear that the actual wing loadings bear 
no simple relation to either of these limiting cases of 
theoretical loading. As a check on the assumption of 
Beskin upwash for these cases, the upwash produced 
along the wings by the finite body was computed; the 
comparison of these values with the Beskin upwash 
shown in Fig. 14 indicates that while there are differ- 
ences they are generally small. For practical pur- 
poses, the calculation of wing lift can probably be car- 
ried out by using the Kirkby-Robinson formula, Eq. 
209), in the region outside the root and tip Mach lines, 
and applying simplified corrections in these regions. 
For example, examination of Fig. 9 shows that an as- 
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sumption that the average pressure in this region at 
any chordwise position is 90 per cent of its value on 
its value on the Mach line would give good results for 
total lift at least up to x = O.7RV JI? 2 
larly, simple results have been found for the experi 


Simi- 


mental results on low-aspect-ratio delta and diamond 
wings with supersonic edges reported by Dvye**; these 
wings have the same chord as the rectangular wings of 
Fig. 13. For these wings too (for the case of body and 
wings at a common angle of attack), it was found that 
the ratio of the measured pressures to the theoretical 
pressures on a delta wing alone did not vary by more 
than 10 per cent from the value at the Mach line inside 
the root Mach line. 

In the formulation of approximations for wing-body 
interference problems, it is often useful to have avail- 
able some general ideas as to the carry-over of lift from 
the wing to the body in the region of the wing rather 
than all along the body. Lagerstrom and Van Dyke 
have used what they have called ‘“‘lift preservation 
theorems” based on planar configurations to arrive at 
such ideas. The reverse-flow theorems for wing-body 
combinations may be applied to arrive at such results 
for nonplanar configurations. Consider a rectangular 
wing of large span mounted on a body of revolution. 


Assume also that the wing is at angle of attack and the 
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WING-BODY 


body at zero. Then any amount of lift which ts lost by 
the exposed wings (relative to the two-dimensional 
value) is recovered by a region on the surface of the body 
which is bounded by the intersection of the body with 
the plane Mach wave from the trailing edge at the rear 
and by the intersection of the body with the Mach cone 
from the wing root at the front. This region of ‘“‘lift 

The proof of this 


26). In 


preservation”’ is shown in Fig. 15. 
result may be obtained by the use of Eq. 
reverse flow, consider the wing at unit angle of attack, 
but without the body. 
ponents of velocity at the location of the surface of the 


This induces normal coim- 


body of amount equal to the z direction cosine of the 
normal to the body at every point. In direct flow, 
consider the exposed wings at unit angle of attack and 
The velocity 
The re 


with the body at zero angle of attack. 
normal to the body surface must be zero. 
verse-flow theoreim is applied to the region outside the 
exposed wings and body contour. Eq. (26 
directly to the result that the total lift on the be dy plus 
exposed wings in direct flow must equal the lift carried 
Contri- 


then leads 


by the exposed wings alone in reverse flow. 
butions to the integrals in Eq. (26) can only arise from 
the region of perturbed flow common to both the direct 
and reverse cases. This explains the particular shape 
of the bounding curves on the surface of the body. If 
unalogies to the planar case are employed, the entire 
region of the body influenced by the direct flow up to 
the trailing edge is the region ‘‘preserving’’ the wing 
lift. On the other hand, the nonplanar theorem shows 
that the lift on the body up to the wing-trailing edge is 
somewhat greater than that necessary to make up for 
iny loss in two-dimensional lift on the exposed wings. 
Lagerstrom and Van Dyke,”! and Morikawa” have 


used the planar analogs of the ‘lift’ preservation” 
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theorem as the basis for calculating the lift on a wing 
body combination without afterbodies by strip theory. 
They have then used slender-body theory to obtain an 
opposite limiting case for comparison, and have at 
tempted to infer conclusions for intermediate cases by 
interpolation between these limiting cases. This pro 
cedure provides a reasonable estimate of the total lift, 
but provides no information relative to lift on the indi 
vidual components. Further, the prediction of mo 
ments for configurations with long after bodies is not 
possible by such methods. For wing-body combina- 
tions with highly swept wings, slender-body theory 
probably gives good results both for total lift and for 
the components lifts. In this theory, however, the 
entire body lift is supposed to occur between the leading 
and trailing edges of the wings. For configurations of 
only moderate sweep, this may result in sizable errors 
in moment. A promising attack on the problem of 
‘‘not-so-slender’’ wing-body combinations is being made 
at Cornell University by means of an extension of the 


Adams-Sears second-order theory." 
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On the Stability of a Vortex Sheet in an 


Inviscid Compressible Fluid’ 


Suin-T Pari 
Vniverstty of Maryland 


SUMMARY 


rt bility of a vortex sheet in an inviscid compressible fluid 

h respect to two-dimensional small disturbances is investi 

g A new simple criterion for complete stability of this case 
fe l rhe significance of ‘‘supersonic,”’ “‘sonic,’’ and ‘‘sub 
disturbances is discussed. It is found that only the 

ral stable disturbance” with respect to the frec 


supersonic neut 


stream gives continuous eigen values. The determinantal equa 
tion for “‘supersonic’’ damped or self-excited disturbance with 
rt ct to the free stream here may be reduced to the same simpk 
form as those for other subsonic or sonic disturbance and the 
e1g¢ values so obtained are discrete 


INTRODUCTION 


y rHE stTupy of the stability of flow of a compressible 
fluid, one of the most interesting problems is that 
the laminar flow of very high Mach Number tends to 
be completely stable with respect to small disturbances. 
To show such a complete stability, it is customary to 
discuss either the smallest eigen value" ° of the deter 
minantal equation which in general has an infinite num- 
ber of eigen values, or the nonexistence of the ‘“‘subsonic 
neutral stable disturbance’ with respect to free stream 
velocity.! The following questions have never 
been fully clarified: (1) if the smallest eigen value shows 
a complete stability, what would be the effect of other 
higher eigen values? In other words, is it possible to 
show that all the eigen values of the determinantal 
equations give the possibility of complete stability? 
2) We know that for the supersonic neutral disturb- 
ances with respect to the free stream the eigen values 
are continuous and not discrete. But is there any pos- 
sibility of the existence of damped or self-excited super- 
with the 
If so, what will be the eigen values of these disturb- 


sonic disturbances respect to free stream? 


ances? In this paper, we study one of the most simple 
examples of the stability of compressible flow, i.e., the 
a vortex sheet or ideal jet in an inviscid 

fluid 


stability of 


compressible with respect to two-dimensional 


Received September 1, 1953 
his paper was written when the author visited the Aero 


nautical Engineering The State 
College during August of 1953, supported by ONR Contract 


Department of Pennsylvania 


Nonr 656(01 The author wishes to express his appreciation to 
Professor Martin Lessen for his interest in this problem and to 
Mr. John A 


Fox for his help in numerical calculations and the 


drawings 
+ Associate Research Professor, Institute for Fluid Dynamics 


ind Applied Mathematics 


w~ 


small disturbances in order to throw some light on the 
above questions 

Phe problem of the stability of a vortex sheet in an in- 
viscid compressible fluid has been studied by Landau‘ 
and Hatanaka® who both the 


bility of complete stability at Numbers. 


showed 
Mach 
Here we find a new simple criterion for the stability of 
Eq 


(Fig. 2). We also discuss the significance of the super 


Hiroshi 


] yOSS1- 


high 


26) | and 


the vortex sheet in a compressible fluid 


sonic disturbances with respect to the free stream for 
this problem which has not been discussed before. 


FUNDAMENTAL EQUATION AND SOLUTIONS 


We consider a straight vortex sheet between two 
parallel streams within the domain y y, and y Y: 
where y; and y. may be finite or infinite (Fig. 1 The 
velocity in the upper half plane (y > 0) is () and that 


0) is l The vortex sheet 
is situated along the x-axis at the beginning, f = 0. 
We assume that the flows on both sides of the vortex 


sheet are potential flows because it is known that no 


in the lower half plane (y - 


vorticity can be produced in the interior of an inviscid 
fluid which is initially in irrotational motion and where 
the stagnation enthalpy and entropy are constant; 
hence the vorticity remains concentrated in the vortex 
sheet. The velocity potentials of the flows on both 
sides of the vortex sheet may be written as 

®@, = l xT @ l 
where 7 = lor2. Subscript | refers to flow in the upper 
half plane (y > 0 


Q). ¢ 


and subscript 2, that in the lower 
half plane (y - are the perturbed velocity po 
tentials which are produced when the vortex sheet is 
slightly disturbed so that it has the shape 7 n(x,t). 
The differential equation which governs the small 


disturbances ¢; 1s 


OxOl 
Ox" T O-¢ 


O-¢ 
(Oo ¢ 


of?) + 21 
M ,? 


1 /a;* [(O°¢ 
Oy") (2 


(1 ~ 


where a; is the sound speed corresponding to the basic 
flow velocity U’; and J, = U’,;/a,, the Mach Number of 
the basic flow. 

Since the disturbance 7 = 7(x,t) is two-dimensional 
in nature, we may assume that 


x(x ci (2 


¢, = ¥, (ye 3 


where V,(y) is the disturbance amplitude function, a is 
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Fig! Vortex Sheet in Inviscid Compressible Fluid 


the disturbance wave number, C is the disturbance 


phase velocity which may be complex, 1.e., 
C=C,+7C (4 
Substituting Eq. (3) into Eq. (2), we have 
VY,” — a’o;7¥,; = 0 (5) 
where 
¢,;=1- [(U — C)/a;|? (6) 
and prime means differentiation with respect to y. 
The general solution of Eq. (5) is 
2 cosh ad V (7) 


V,(y) A, sinh ac;vy + A 


where Aj; and Aj» are arbitrary constants. Since Eq. 


2) is a linear equation, the method of superposition may 


be used here to find complicated solutions. 


SUBSONIC, SONIC, AND SUPERSONIC DISTURBANCES 


It is convenient in the study of the stability of flow of 
a compressible fluid to classify the disturbances as 
“supersonic” relative to a 


“subsonic,” “‘sonic,’’ and 
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certain local velocity when the x-component of the 


phase velocity of the disturbance relative to the local 
velocity is less than, equal to, or greater than the mean 
speed of sound at that point, respectively. In the pres- 
ent problem, such disturbances may be classified with 


respect to ( For such classified disturbances we have 


(1) C; = O neutral stable disturbance 


pure imaginary supersonic | 
: with respect 
\ to l 


a; = <zero—sonic 


real—subsonic 


(2) U7; = C 
o; = real, for all cases. 

(3) C; # 0 damped or self-excited disturbances 
a; = complex number, for all cases. 


The significant result is that only for supersonic 
neutral stable disturbances with respect to L’;, is a; pure 
imaginary which will give continuous eigen values. 
For all other disturbances, it will be shown later that 
the eigen values for the stability of a vortex sheet are 


discrete. 


BOUNDARY CONDITIONS 


The boundary conditions on the vortex sheets are: 
(1) the velocity components of the fluid perpendicular 
to the vortex sheet are equal on both sides and equal 
to the movement of the vortex sheet itself in that direc- 
tion. If the higher order terms are neglected, this con- 


dition may be written as: 
On/Ot = (O¢,;/Ov) — [U;(On/Ox S 
where the expression of 7 may be written as 
wet) = Be- si (9 


where B is a small constant. The expression of Eq. (9 
is assumed because of the expression of Eq. (3). 
(2) The other boundary condition is that the pres- 


sures on both sides of the vortex sheet are equal, i.e., 


¥1 E O¢) L =] = Y2 fe: Ov» L _ (10 
ay" Ox Of as" Ov Of | 
where y; = the ratio of specific heats (C, C,) on each 

side of the vortex sheet. 

As to the boundary conditions at y; and ye, the fol- 
lowing three different problems may be studied: 

(1) Vortex sheet between solid walls at y = y,; and 
y = yo. In this case, the disturbance velocity com- 
ponents perpendicular to the walls are zero, 1.e., 


O¢) oy = § at = yy = li 
Og./dy = 0 aty = w= —l (11) 


where /; and /, are positive real numbers. 

(2) Two-dimensional jet with y = y. as the axis of 
symmetry of the jet. It is well known that in this 
case the disturbances may be separated into antisym- 
metrical and symmetrical disturbances. These dis- 


turbances may be treated independently. For anti- 
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the conditions 


it 1 y; are the same as those of Eq. (11 


syminetrical disturbances, boundary 
while for a 
symmetrical disturbance, the value of Ov Oy must be 
zero at V ye, and then the boundary conditions at 


are 
0g, Ov = VU at y V1 l 
Og. / Ov" 0 at y y —/. 12 
EIGEN VALUE EQUATIONS 


into the 
boundary (12 
with the help of Eq. (9) the eigen value equation may 


If we substitute the general solution Eq. (7 


conditions Eqs. (8), (10), (11 or 


be obtained. For a vortex sheet and an antisymmet- 


rical disturbance of a two-dimensional jet, we have 


; 01 ¥i(1 ace o,7) |! tanh acy); = 


— }o2/[y(1 o2”)]} tanh aay /, (13) 
For a symmetrical disturbance of a two-dimensional jet 
we have 


male. = a1") |} 


\ oy tanh ag,/; = 


(14) 


= } a» [yo1 — a”) ]j coth aa»/. 


CASES FOR SEMI-INFINITE OR INFINITE DOMAIN 


It is usually desirable to study the stability of the 
vortex sheet itself without the influence of the solid 
wall y = yj. Then we may assume that /; tends to be 
infinite. We have to find the limiting value of tanh 
ag,l; as /; tends to be infinite in Eqs. (13) and (14). 


Now 


tanh ao,/; = (tanh ag;,/; + tanh tag;,/;) + 


(1 + tanh ao,,/; tanh taq;;/;) (15) 
where 
0; = 6; +10 (16) 
If o;, + 0, we have 
iim tanh ag,l; = +1 (17) 
If o;, = 0, we have 
lim tanh aol; = lim 7 tan ao;(l (18) 


lj— lj— 


Since tan © may be any value, the limiting value 
of Eq. (18) may be any value. 

From the analysis of the classification of disturbances, 
we know that only for the supersonic neutral disturb- 
ance with respect to U’;, is the value of o; pure imagi- 


nary. Only for such supersonic neutral disturbances, 


the eigen value a may be any value, i.e., both Eqs. (13) 
For all the 


other disturbances, including supersonic damped or 


and (14) give a continuous spectrum of a. 


self-excited disturbances with respect to l’; the eigen 


values are discrete and may be obtained from Eqs. (13) 


or (14) with the help of Eq. (17). 


A 


VORTEX SHEET 324 


VORTEX SHEET IN AN INFINITE DOMAIN OF FLUID 


For a vortex sheet in an infinite domain of fluid, both 


and tend to be infinite. With the exception of 


supersonic neutral disturbances with respect to L’,, Eq 
13) becomes 
Go 71 | [+ 1Q 
Now we assume that 7 y, and put 
AM; l C)/a,, M l C) a 20 
Eq. (19) becomes 
V7. Ve+ VM 2] 
where JJ; + Jo. 
We write 
C 
| 
l 
99 
( = C\ . (U, T U 2 \ e 
U = (U, — U2)/2 M140) Mo: 
then Eq. (21) becomes 
A‘ — [2 + (1/m,?) + (1/my,?) A: 2/(1/m 
(1 /mo?)|A + [1 — (1/m,7) — (1 /me* 0 23 


The complete stability condition is that all of the fou 
roots of Eq. (23) must be real. m,* > 
0, the sufficient conditions that all the roots of Eq. (23 


Since 1/7," + |] 


arerealare: (a) 


(24) 


= 
= 
V 


(1/m,*) + 


) 


and (b) the discriminant A of Eq. (23) is greater than 


zero,® 1.€., 


A>O 25) 

By Eq. (24) and the Descartes rule of sign, we see 
that Eq. (23) has at least two real roots, because with- 
out loss of generality, if we take 1/m,* 1/m2 > 0, 


there will be one positive real root since there is only one 


change in sign of Eq. (23) and one or three negative 
But A> 0O 


means that Eq. (23) has four distinct real roots or two 


roots since there are three changes in sign. 
pairs of complex roots. Combining these two condi 
tions, we have that Eq. (25 
The discriminant of Eq. (23) after simplification may 


has four real roots 


be written as follows: 


A = [64 — (48/m, 1S 71 + (12/m,4 
(S4/2,71127) + (12/23 1 my' - (3/m,4m") — 
(3/424) — (1/me®) |] 64/m,2m.? (26 


The boundary line of stability in 77, m2 plane* is 


* Dr. F. W. Geiger of Cornell Aeronautical Laboratory has in 
formed the author that he found by geometrical considerations a 
Even though it is diffi- 


criterion in place of Eqs. (26) and (27 


cult to reduce Eq. (26) into the formula of Dr. Geiger analytically, 
by numerical calculation (Fig. 2) the result of Eq. (26), A = 0, is 
identical to that given by the formula of Dr. Geiger whose result 


will be published by Cornell Aeronautical Laboratory 
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A = 6 (27) 


which is shown in Fig. 2. Since the condition of Eq. 


(24) always holds for the point just on the right side of 
A = O curve in Fig. 2, the region on the right side of 


A = Ocurve is the stable region, 1.e., in this region, the 
vortex sheet is completely neutrally stable with respect 
to small two-dimensional disturbances. Only neu 
trally stable disturbances can occur in this region with 
definite velocity of propagation for subsonic neutral 
disturbances with respect to Ul’; and with all kinds of 
velocity of propagation for supersonic neutral disturb- 
ances with respect to U;. The stability condition 
Eq. (27) also holds for (1/m,") + (1/m»?) < 1 which 


shows also that the four roots may be all real. 


When m, = mo = M, 1.e., Qj = do = a, we have the 


complete stability criterion as 
a> YZ (28) 
When either one of m, and my is infinity, the other 


for stability condition will be larger than one half, 
re. 


ms > 1/2 if mM, = (29) 


The interesting result in this section is that it is pos 
sible to show that all the eigen values of the deter 
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minantal equation show complete stability for flow of 


high Mach Numbers, 1.e., large value of m,. 

The above analysis gives the same order of magni- 
tude of 
(U,-U 


3 where the stability condition is (Uy — U2)/a > 2.5 


Mach Number for complete stability, ice. 
2)/a > 2V 2 as the analyses of references 2 and 
Here we do not use the unconvincing argument of the 
nonexistence of subsonic neutral stable disturbances 
with respect to free stream but solve the complete de- 
terminantal equation which includes the “supersonic” 
damped or self-excited disturbances with respect to the 
free stream if they exist. 
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Postbuckling Behavior of Axially 
Compressed Circular Cylindrical Shells’ 


JosEPH KEMPNER? 


Polytechnic Institute of Brooklyn 


SUMMARY 


an axially compress¢ d 


The 
thin-walled circular cylindrical shell loaded either by dead weights 


postbuckling characteristics of 


or by a rigid testing machine are determined It is shown that 
for either loading condition the minimum applied stress in the 
postbuckling region is 0.182(#t/R) and that the region of stable 
equilibrium corresponding to loading by the rigid testing machine 
includes and extends beyond that obtained with dead weight 
loading. The work here described is a continuation of work done 
earlier by von Karman and Tsien, by Michielsen, and by Leggett 


and Jones, 


INTRODUCTION 


— EXPLANATION OF THE discrepancies between ex 
perimentally determined failing stresses and those 
computed by the small deflection theory for axially 
compressed thin-walled circular cylindrical shells ap- 
pears to have been attempted first with the aid of a 
large deflection theory by L. H. Donnell! in 1934. In- 
cluding second order terms in the expressions for the 
Donnell 
equation 


median surface strains, derived a_ second- 


approximation compatibility which with 
energy considerations was applied to the determination 
of the postbuckling behavior of thin-walled circular 
cylindrical shells under axial compression. However, as 
pointed out by Donnell and Wan,’ the analysis was 
oversimplified and hence led to an unsatisfactory solu- 
tion. 

In 1941 Th. von Karman and H. S. Tsien* presented 
second approximation compatibility and equilibrium 
equations for the same problem. It is shown in refer- 
ence 4, and in the present text, that their equilibrium 
equation was incomplete. However, these investiga- 
tors did not use the equilibrium equation but, like 
Donnell, applied the principle of stationary potential 
energy. The results presented in reference 3 indicated 
that a buckled cylinder could be maintained in equilib- 
rium under an axial stress considerably below the 


critical stress computed from small deflection theory. 
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In the past ten years the work of von Karman and 
Tsien has been refined and extended by several investi 
gators.* In references 5 and 6 the compatibility 
equation, energy expressions, and assumed deflected 
the 


shape were accepted from reference 3 and prin 
ciple of stationary potential energy was again applied. 
However, in contrast to the analy sis of reference 3, in 
references 5 and 6 the energy was made stationary with 
respect to all arbitrary deflection parameters, whereas 
in reference 3 parameters characterizing the number of 
circumferential waves and the ratio of the circumfer 
ential and axial wave lengths were fixed. From the 
more rigorous analysis it was found that the smallest 
applied axial compressive stress o required to main 
tain the cylinder in a state of equilibrium with finite 
wave amplitude corresponded to oR Et 0.195. This 
value is considerably less than the critical value of the 
classical linear theory which is cR E:t = 0.605 

In reference 2, the analyses of references 5 and 6 were 
extended to include the effect of an additional arbi- 
trary parameter in the assumed deflected shape, as well 
as the effect of various types of initial imperfections in 
the cylinder. However, the solutions of the resulting 
set of equations for the deflection parameters were cal 
culated rigorously only for loads between 100 and 75 
per cent of the classical critical load; at lower loads cer- 
tain restrictions on the parameters were introduced. 
With these approximations the minimum postbuckling 
value of (cR Ft) for an initially perfect cylinder was 
found to be approximately 0.26 (Fig. S, reference 2) 
Since the investigators were primarily interested in the 
greatest axial load a cylinder with given imperfections 
could resist, the postbuckling region investigated was 
considerably smaller than that of references 5 and 6. 

The present investigation is an extension of Michiel 
sen’s and Leggett and Jones’ continuation of von Kar 
man and Tsien’s work.® 5 * The energy expressions 
of reference 3 are used to derive the corresponding 
The equations for equilibrium 
the 


equilibrium equations. 
stresses are then used in determi 
With the 


expression for the radial displacements at buckling ob 


of membrane 
nation of a suitable compatibility equation. 


tained from classical linear theory as a first approxi 
mation, the equation of equilibrium of radial forces 
and the compatibility equation are iterated to obtain 
This approximation 


includes the deflected shape assumed in reference 2 


a suitable second approximation. 


which in turn includes the shape assumed in refer- 
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ences 3, 5, and 6. In order to make tractable the com 
putations involved in the energy solution of the prob- 
lem, only that part of the second approximation for the 
deflected shape which corresponds to that assumed in 
reference 2 (in that region where the parameters were 
unrestricted) was considered. Thus in the present in 
vestigation one more arbitrary deflection parameter 
is introduced than in references 3, 5, and 6. 

The analysis shows that, for a cylinder axially com 
pressed by either a dead weight or a rigid testing ma 
chine, the minimum value of (oR £t) consistent with 
finite-deflection equilibrium configurations is 0.182. 
Furthermore, it is seen that the equilibrium configu 
rations are identical for the two types of loading. How 
ever, the results of a stability analysis indicate that the 
rigid testing machine affords a range of stable equilib- 
rium including and extending beyond that obtained 


with a dead-weight loading. 


SYMBOLS 


Ay, Aig, Agi = functions of wave-length ratio 

D = (1/12)Et?/(1 — v?) 

E = modulus of elasticity 

F(x,y) = airy stress function 

y = shell length 

R = mean shell radius 

U;, U2, U3 = nondimensional energy quantities 


extensional strain energy, bending strain 
potential of applied axial load, 


Wi, We, Ws, Ws 


energy, 
potential of applied lateral load, respec 
tively 
“ab 6, 2 deflection parameters 
dy, dx, doz = stress function parameters for first ap 
proximation 
fir, foo, fi ’ . 
: : = stress function parameters 
Sais Foo: Fon 
m,n = numbers of axial and _ circumferential 
waves, respectively 
p = inward radial pressure 
t = wall thickness 
u,v, u = axial, circumferential, and radial displace 


ments, respectively 
y = axial and circumferential coordinates on 
median surface of shell, respectively 
= median surface shear, axial, and circum 
ferential strains, respectively 


€ = unit end shortening 

n = (t/R)n? 

© = axial and circumferential half wave lengths, 
respectivels 

M = (A, A,) 

v = Poisson’s ratio (0.3) 

Fs = deflection parameter 


= axial, circumferential, and shear median 
surface stresses, respectively 


= applied average axial compressive stress 


o 
w = £7 
vi = (0'/Ox') + 2(04/dx? Oy?) + (04/dy! 


UNDERLYING ASSUMPTIONS 


In the ensuing analysis the thin-walled circular 
cylindrical shell of uniform thickness is assumed to be 
stressed elastically and to be of sufficient length such 
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that the effects of end restraints are negligible N 
Nojima and S. Kanemitsu have shown that the latter 
assumption is valid if the length is greater than 1.5 
The results of their investigations 
With only axial stresses 


times the radius. 
are described in reference 7. 
applied to the shell, omission of the effects of end re- 
straint leads to the condition of zero average circum- 
ferential and shear stresses. 

In addition to the foregoing assumptions, the changes 
of curvatures in the axial and circumferential directions 
and the twist of the median surface are assumed to be 
(O°w/Ox*), (O°w Oy"), and (O*w/Oxdy), respectively, 
where w is the inward radial displacement and x and y 
the coordinates of a point on the median surface of the 
shell (Fig. 1). Justification of the use of these simple 
expressions was given by Donnell in reference 8. 


BaAsIc RELATIONS 


The following relations upon which the present 


analysis is based are given in reference 3. 
Median Surface Strains 


With terms up to the second order included, the axial 
and circumferential median surface strains e, and €¢,, 
and the median surface shear strain y,, of a thin circular 
cylindrical shell of radius R are 


€- = (Ou/Ox) + (1/2) (Ow/Odx)? 
€, = (Ov/Oy) + (1/2) (Ow/dy)? — (w/R) (1) 
Yry = (Ou/Oy) + (Ov/Ox) + (Ow/Ox) (OW oy) 
in which , v, and w are displacements of a point of the 
median surface in the axial, circumferential, and inward 
radial directions, respectively (Fig. 1). 
Median Surface Stresses 

The axial, circumferential, and shear stresses in the 
median surface are from Egs. (1) and from the usual 
form of Hooke’s law for thin shells 


[E/(1 — v?)] } (Ou /dx) + 
(1/2) (Ow/ Ox)? + v[(Ov/Oyv) + 
(1/2) (Ow/Oy)? — (w/R)]j 
o, = [E/(1 — v?)]}(Ov/dy) + e 
(1/2) (Ow/dOy)? — (w/R) + 
v[(Ou Ox) + (1/2) (Ow /Ox)?} 
tr, = [K/2(1 + v)]} (Ou/Oy) + (Ov/Ox) + 
(Ow Ox) (Ow /Oy) { 


Or — 


in which o, and o, are positive in tension; / is the modu- 
lus of elasticity; and v is Poisson’s ratio (0.3 in subse- 
quent calculations). 
Extensional Stain Energy 
The elastic energy corresponding to the median sur 
face stresses can be expressed as 
"Ll  (°2eR 


W, = (t/2E) [((o, + o,)? — 


0 0 
211 + v) (a,0, — Try?) |dx dy (3) 
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POSTBUCKLING BEHAVIOR OF 
in which L is the length and ¢ the wall thickness of the 
shell 

Bending Strain Energy 


With the 
“Underlying Assumptions,’’ the energy of bending be- 


curvatures and twist mentioned under 


comes 


"I "2eR 
W D/2) / / ) [(O7zw ‘Ox 
1 — 


v) [(0°?w/Ox?) (07w /Oy") 
(O°w /Ox Oy)2]!dxdy (4 


(O-w Ov) |? 


in which D (1/12) k#t/(1 — vr? 
Potential of External Loads 
rhe potential energy of the axial stresses applied to 


the ends of the shell can be expressed as 


*2rR “2 
I’ —f / (or)r-1 AY / Ou, Ox) dx (5 
J 0 J 0 


The potential energy of an applied external radial 


pressure Pp 1s 


"Z *2rR 
W;, - / / pw dx dy (6 
se << J0 


If a variation of the total potential is contemplated 
for axial loads applied by a rigid testing machine, then 
SIT’ 0 since no axial displacements are permitted. 


EQUILIBRIUM AND COMPATIBILITY EQUATIONS 


Addition of Eqs. (3) to (6) and subsequent variation 
of the total potential yields the following three equilib 
rium equations, as well as boundary conditions which 


are not given here 


D t)Viw = (p/t) + (0/0x) (¢,0w/Ox) + 


(0/Oy) (7,,0w/Ox) + (O0/Ox) (7,,0w/Oy) + 
(0/Oy) (¢,0w/Oy) + (a,/R) (7 
(Og,/Ox) + (O7z,/Ov) = O (S 
(Or,,/Ox) + (O0,/Oy) = O (9 
in which V* = (01/dx‘) + 2(04/Ox? Oy”) + (04/Oy"). 
In the derivation of Eqs. (7) to (9), Eqs. (2) were 


used to relate the stresses and displacements. It may 
be noted that the energy term IV’; of Eq. (5) does not 
contribute to Eqs. (7) to (9), but rather enters into the 
expressions for the boundary conditions. 

Elimination of the stress derivatives obtained upon 
performance of the operations indicated in Eq. (7) 
with the aid of Eqs. (8) and (9) yields the following 


equilibrium equation :* 


D t)Viw = (p/t) + o,(0°w/Ox*?) + 


27,,(0°w/Oxdy) + a,(0°w/Oy?) + (co, R) (10) 


Eqs. (8S) and (9), which the conditions 
of equilibrium of forces in the median surface, are 
identically satisfied by the introduction of the Airy 


stress function F(x,y) defined by the relations 


express 
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COORDINATES AND COMPONENTS OF 
DISPLACEMENTS OF A POINT ON 
MEDIAN SURFACE OF AXIALLY 
COMPRESSED THIN-WALLED 
CYLINDER 


FIG. | 


(0?F /Ox dy), 
o O°F On 1] 


5 o°F dy’), r 


Combination of Eqs. (2), (8), (9), and (11) yields 


(O-w Ox Ow OY 
(1/R) (0°w /Ox? 


V'iF = E[(0°w dx Oy 


12 
This compatibility equation was first developed in refer- 
ence |. 

With the aid of Eqs. (11 
form of the equilibrium Eq. (10) is found to be 


and (12), an alternative 


(D tV'w + (E/R?*) (04w Ox) 


V[(p/t) + (O2F dy?) (O%w/ Ox 
2(0°F/Ox Oy) (O°w Ox Oy) (0°F/Ox") X 
(O°w Oy") | ER) (0? Ox?) [(0?w/Ox Oy 
(0°w/dx?) (O°w Oy?2)] (13 


Except for the last term on the right-hand side, Eq 
(13), with p = 0, is identical to von Karman and Tsien’s 
Eq. (12) of reference 3. It was indicated in reference 
! that omission of this term, which is of the same order 
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FIG. 2 VARIATION OF PARAMETERS a,b,e,¢ AND w 
WITH WAVE-LENGTH RATIO y 
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FIG. 3 VARIATION OF APPLIED STRESS, UNIT END 
SHORTENING AND WAVE NUMBER 
PARAMETER WITH WAVE-LENGTH RATIO 


of magnitude as those remaining, may lead to consider 
able error in the analysis of the postbuckling behavior 
of axially compressed cylinders. 

The form of the equilibrium equation used in the 
present analysis is from Eqs. (10) and (11) with p = 0 


(D/t)Viw = (0°F/ Oy") (O°w/ Ox") — 
2(0°F/Ox Oy) (0?w/Ox Oy) + 
(0°F Ox”) (O?w Oy") + (1/R) (0?F/Ox") (14) 


APPROXIMATE DEFLECTED SHAPE 


The well-known deflection function which satisfies 
the linearized form of Eq. (13) can be expressed as 


w, = &t[cos (wx/dy,) cos (ry/M1,) + e1] (15) 


in which &, @:, Az, and Ay, are independent of x and y, 
and the subscript 1 indicates a first approximation. 
Substitution of this approximation to w into the right- 
hand side of Eq. (12) and subsequent solution of the 
resulting equation yields 


l 


F, = dy cos (rx/4,) cos (ry/Ayy) + 
doy COS (27xN/Aiz) + 


dos COS (27ry My) (o:/2)y? (16) 


in which d is a parameter and go; is the average applied 
axial compressive stress. This expression for F satis- 
fies the condition of zero average circumferential and 
shear stress and hence, as discussed earlier, reflects the 
neglect of edge restraints. Introduction of the ex- 
pressions for w and F from Eqs. (15) and (16) into the 
right-hand side of Eq. (14) and determination of a par- 
ticular solution of the resulting equation lead to the 
following form for a second approximation to w: 


w= Et[cos (wx/Xzx) cos (ry/d,) + 
a cos (27x/dz) + bcos (2ry/d,) + 
c cos (rx/X;z) Cos (Sary/A,) + 


d cos (32x/X,) cos (rv/A,) +e] (17) 


in which é, a, }, c, d, and e are undetermined coefficients, 
and X, and X, are the half wave lengths in the axial and 


circumferential directions, respectively. That Eq. 
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(17) represents a suitable choice for the form of the de- 
flection function of an axially compressed thin-walled 
circular cylindrical shell is seen from reference 3 Che 
deflection function assumed in this reference which was 
arrived at from physical observations is equivalent to 
Eq. (17) ifc = d = 

For the present analysis Eq. (17) 


0 and a b. 

will be applied with 
( d = OQ, and thus the undetermined quantities are 
£,a,b,e,X,,andA,. The resulting simplified expression 
for w corresponds to that used in the analysis of refer 
ence 2. In the present paper all the deflection param- 
eters are determined from the principle of stationary 
potential energy whereas in much of the work presented 
in reference 2 1t was found necessary to assume fixed 


values of a, b, and (A,/X 


STRESS FUNCTION 


The expression for F obtained from the solution of 
Eq. (12) with w introduced from Eq. (17) is for « | = 


QO: 


F = —(Et?/pn*) [fu cos (4x/A,) cos (ary A 
fo. cos (27x/X,) cos (27d 
Jizs CoS (wx/A,) Cos (Sry XA 

fu cos (38ax/Xz) cos (rv /XA,) + 

foo COS (27x /X;,) + 


ufo. cos (2rry/dy)] — (o/2)v? (18 


in which o is the applied average axial compressive 


stress, 
fi — Aj,0[2w(a a b — 1] 
foo = Auw*ab 
Tis — 2A wb 
fa = 2A 1qw"a 
fo = (1/32)w(w — Sa) 
for = (1/32)w? 
and 
w = £m, w= (Ay/Az), n = (t/R)n 
(r?Rt/d,") 
Ay = p*/(1 7” ie)", Aig Le (] + Oy? E 


It may be noted that yu is the ratio of the circumfer- 
ential and axial wave lengths and » is the number of 
circumferential waves. The number of axial waves 
m is related to A, by the expression 777 L,/2hz. 


Unit END SHORTENING 


Since the operator of a standard type testing machine 
has control only of the distance between movable and 
fixed platens, the unit end shortening of an axially com- 
pressed shell is one of the principal parameters of the 


present problem. It is defined as*® 


71. 
e= —1/2) / (Ou /Ox) dx (19) 
J 0 
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and is the ratio of the total shortening in the axial direc- 
tion to the cylinder length. From the first two of Eqs. 
(2) and the first and third of Eqs. (11), 


(0u/Ox) = (1/E) [(0?F/dy?) — 
v(O?F/Ox?)] — (1/2) (Ow/Ox)? (20) 
Hence from Eqs. (17) to (20) forc = d = 0, 
(eR/t) = (cR/Et) + (1/8) (w7n?/n) (1 + 8a?) 
(21) 


UNIFORM INWARD DISPLACEMENT 


That component of the total displacement which 
corresponds to a uniform inward displacement is given 
Consideration of the condition that 
must be a periodic function of y 


in Eq. (17) as éte. 
the displacement v 
readily leads to the determination of e (reference 3). 


From Eqs. (2) and (11) 


(1/E) [(0?F/Ox?) — v(0?F/dy?)] — 
(1/2) (Ow/Oy)? + (w/R) 


(v/Ov) = 


(22) 


Eq. (22) together with Eqs. (17) and (18) and the 
periodicity condition yield (for ¢ = d = 0) 
(1/8) (w?/n) (1 + 8b?) — v(oR/Et) (23) 


fe = 
ceé 


in which € = w/yn. Thus in Eq. (17), e is not an inde- 


pendent parameter. 


DETERMINATION OF THE DEFLECTION PARAMETERS 


The deflection parameters of Eq. (17) can be ex- 
pressed in terms of a, b, n, py, and w. These parameters 
are determined from the condition that at equilibrium 
their values are such that the total potential energy 
of the system has a stationary value for any small vari- 
ation of the parameters. This condition can be en- 


forced for ¢ = constant (dead-weight loading) or € = 
constant (rigid testing machine loading, 6IV3; = 0). 
However, it can be readily verified that these two con- 
ditions lead to the same expressions for the deflection 


parameters. 
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FIG. 4 VARIATION OF AXIAL AND CIRCUMFERENTIAL 
HALF WAVE LENGTHS AND NUMBERS OF 
WAVES WITH WAVE-LENGTH RATIO 


(17), (18), and (20), 


the three energy quantities entering into the calcula- 


From Eqs. (3), (4), (5), (11), 


tions become 


U; —= (w? n”) (wfAn[(a + ))? + la b? | > om 
A,a* + Agb>? + (pi/128)} + 
Ay [(1 1) — wl(a + b) | + (1/128) (eo — Sa)?) "es 
(oR/Et)? (24) 
U, = (1 {(1 — vy?) |w?[(1 + py 
32(uta? + b?)] (25 
U3 = —2(cR/Et) [(cR Et) 4 


(] S) (w"u n) (1 
in which 
U, = (R/rEPL)W, ete. 

Expressed in terms of a, the following five equations 
are obtained from enforcing the condition that the sum 
of Uj, U2, and Ll’; has a stationary value with respect 
With 


these relations can be expressed in 


to a, b, n, wu, and w, respectively, for constant a. 
the aid of Eq. (21 
terms of the unit end shortening e rather than the aver 


age axial compressive stress o. 


a= (w/ 16) (SA, +1)-—- Ay wb | )w?[Ay (1 + 467) + Ajo] a 
[2/3(1 — v?) Jun? — 2u2n(oR/Et) + (1/2)} (27) 
b —= (Ay 2)w — Awa | »w?[An(1 — fa-) , a A | _ [2 3(1 a In?; 28) 
9 = (w?} An[(a + b i we $a*b*|] + Aya? + Agb? + (u4 128)} a 
Ay[(1 1) — w(a + b) | + (1/128) (w — Sa ) us (oR Et l Sd 29 
(oR/Et) = $8/{un(1 + 8a2)]} (wf [An/(1 + w2)] [(@ + 8)? + 40°B?] + 
(A 19” (1 + Qu") | a a [9A 9b? (9 + u”) | + (y* 128)| + [Ay 1 + gu“) | 1 4 -wat+ b + 
1/48(1 — v®)] (1 + uw? + 32u2a*)u*n?) = (30 
wAnl[(a + b)? + 4a7b?] + Aya? + And? + (1/128) (1 + wt)} — (3/32) [SAn(a + 6) + a] 4 
(1/8) }Au + 2a? + [1/1201 v?)] [((1 + pw)? + 32(uta? + 6?) |n? — w2n(ocR/Et) (1 + 8a); 0 (31) 
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Eqs. (29) and (31) can be combined to yield an equation independent of c. 


> 


n° = 3(1 — p*) ,4A1 +- Sa? — w[8A i1(a + hb) + al} (1 + a 2 32(p 4a? - b*) | 


tS 


Similarly Eqs. (30), (81), and (52) yield a relation which is independent of ¢ and 7. 


“/2 [(a + b)? + 4a°b?] + 9Aj9 ’a? + Ay 


| 3 
w*) An 


*h? + (1 128)} = w( An 


(a + Bb) + (4/8) 


(u2/16) }(1 + pw? + 32u2a?) [SAn(a + b) + a]/[(1 + w2)? + 32(uta? + b?)]} ) + 


(Ay*?/4) + (a?/2) — )(u?/4) (1 + w? + 3220?) (An + 2a") [(1 + w?)? + 32(u4a? + b2)]j = 0 (3: 


Hence, simultaneous solution of Eqs. (27), (28), (30), 
(32), and (33) for constant o or constant ¢ [Eq. (21)] 
determines the parameters a, 5, n, wu, and w such that 
the total potential is stationary. It may be noted that 
if a, b, n, wu, and w constitute a solution, then so does the 
set —a, —b, n, uw, and —w. However, the deflected 
shape of the shell [Eq. (17) with c = d = 0] is the same 


for the two solutions. 


STABILITY CONSIDERATIONS 


Although any solution of Eqs. (27), (28), (30), (32), 
and (33) corresponds to a state of equilibrium, it does 
not generally correspond to a state of stable equilib- 
rium. For such a solution to represent a stable equilib- 
rium state, it must correspond not only to the vanishing 
of the first variation of the total potential but also to a 
positive value of the second variation of this potential. 
Thus, since a solution of the aforementioned equations 
enforces the vanishing of the first variation of the po- 
tential, the condition that the solution correspond to a 
stable equilibrium configuration is that the second 
variation must be a positive definite quantity. The 
second variation is such a quantity provided that the 
determinant 


U,,| > 0, 1, 7 = a, 6, n, w, and w (34) 
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FIG. 5 VARIATION OF PARAMETERS a,b,e,w AND lu 
WITH UNIT END SHORTENING 
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and each of its four leading subdeterminants is likewise 


positive. The quantity U;; represents the total po- 


tential differentiated twice partially with respect to the 
Although the 
first variation of the total potential led to identical rela- 


parameters indicated by the subscripts. 


tions among the deflection parameters for ¢ = constant 
(iie., OU’ = 6U), + 6U2 + 6U3) and € = constant (i.e., 
6U = 6U, + 6U>) such is not the case for the second 
variation; thus, the set of five stability determinants is 
not the same for 0 = constant and « = constant. 
Hence, in general, the two methods of loading can be 


distinguished on the basis of stability considerations. 


NUMERICAL SOLUTION 


In order to determine the postbuckling characteristics 
of the axially compressed shell, it is necessary to solve 
for the deflection parameters a, b, n, u, and w from the 
set of nonlinear algebraic Eqs. (27), (28), (30), (32), 
and (33) for a range of the applied compressive stress 
o (or unit end shortening ¢). Such solutions can be 
obtained with a suitable numerical successive approxi- 
mation procedure. Since the wave-length ratio param- 
eter appears in these equations in a rather complicated 
manner, it is more expedient to fix uw and solve for 
(oR/Et) along with the remaining deflection param- 
eters than it is to fix (¢R/Et) and solve for all the de- 
flection parameters. Thus for a given value of y, first 
approximations are assumed for a and b, thereby re- 
ducing Eq. (33) to a quadratic in w. A value of w ob- 
tained from the solution of this equation together with 
the given uw and the assumed a and 6 permit the direct 
computation of the corresponding n and (¢R/Et) param- 
Substitu- 


tion of the calculated values of w, 7, and (oR/Et) into 


eters from Eqs. (32) and (30), respectively. 


Eqs. (27) and (28) determines new values for a and } 
In the 


procedure used, for each assumed value of a, several 


which in general differ from those first assumed. 


values of 6 are chosen until the calculated value of } 
[Eq. (28)]| is equal to the assumed value. Thus a set 
of values for a and 6} is determined such that the dif- 
ference Ab between the assumed and calculated values 
of b is zero. Similarly, for each assumed value of a, a 
value of } is also determined such that Aa is zero. In 
these calculations plots of Aa and Ad vs. 6 for each 
assumed a are helpful. To assist in the determination 
of a set of values a, 6 for which Aa and Ab are simul- 


taneously zero, Aa = 0 and Ab = 0-curves are plotted 
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The intersection of these curves deter- 
and 


in the ab-plane. 
determine w, 7, 
(30), (32), and 


mines a@ and #6, which in turn 


oR ft), 


33) are satisfied for the assumed value of yu. 


such that Eqs. (27), (28), 
The en- 
tire procedure must be repeated for a suitable number 
of u-values. By means of Eq. (21) the stress parameter 
sR /:t) can be related to the end shortening parameter 
eR ¢) once a solution of the simultaneous equations 1s 


found 


RESULTS AND DISCUSSION 


The results of the solution of Eqs. (27), (28), (30), 
32), and (33) are presented in Figs. 2 to 4 as functions 
of the wave length ratio uw, and in Figs. 5 to 7 as func- 
tions of the unit end shortening parameter (eR ¢). In 
Fig. 6 the deflection parameter € and the applied stress 
sare compared with the corresponding quantities de- 
termined by Michielsen in reference 6. In this refer 
ence, as well as in reference 5, the minimum value of the 
applied stress in the postbuckling region was found to 
be 0.195(Et/ R) at p 0.40. In the present analysis 
this minimum stress was found to be 0.1S82(/:t/ R) at 
u = 0.362. These results compare favorably, and both 
are considerably below the critical value of classical 
theory 0.605( Et R). 
out that the value of (oR/Et) reaches a maximum of 


0.77, whereas the present solution leads to no such 


In reference 6 Michielsen points 


maximum but rather indicates that for large values of 
eR t) or uw the value of (oR Et) approaches asymptot- 
ically the critical value of classical theory 0.605 (Fig. 5 
or 6). In Fig. 6 the straight line segment of the 
oR Et) vs. (eR/t)-curve drawn from the origin corre- 
sponds to the prebuckled state of the shell. It ter- 
minates at the classical buckling stress 0.605(/:t R). 
Application of the stability conditions discussed ear- 
lier showed that in the postbuckling region for a dead- 
weight loading the solution described corresponds to 
stable equilibrium configurations for the region of the 
curves defined by u« > 0.362. The lower bound of yu 
corresponds to the minimum value of applied stress 
0.1S2(Et/R). However, if the load is applied by a 
rigid testing machine, the region of stability increases 
tou > 0.25, the lower bound of which corresponds to the 
minimum value of the unit end shortening found to be 
0.307(¢/R). Thus from the (¢R/ Et) vs. (eR t)-curve 
of Fig. 6, the region of postbuckling stable equilibrium 
configurations is represented by the portion of the curve 
to the right of the point of minimum stress for dead- 
weight loading and to the right of the point of mini- 
mum unit end shortening for rigid testing machine 
loading. The remaining region of the postbuckling 
curve for each case of loading represents unstable con 
figurations. The results of the stability analysis for 
loading by means of a rigid testing machine are in 
complete agreement with the analysis of Tsien in refer- 
ence 9. In this reference it was found that a column 
supported laterally by nonlinear supports and loaded 
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FIG. 6 VARIATION OF DEFLECTION PARAMETER ¢ AND 
STRESS o0 WITH UNIT END SHORTENING 
AND COMPARISON WITH REF. 6 


axially by a rigid testing machine became unstable when 
de/do = 0. 
which the results of analyses of the type presented here 


Tsien has also proposed a method with 


can be used to determine conditions which correspond 


to the transition from an unbuckled to a buckled state. ! 
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A Guide to the Use of the M.I.T. Cone 
Tables 


RICHARD C. ROBERTS? ann JAMES D. RILEY? 
1.8. Naval Ordnance Laboratory 


SUMMARY 


M.1.T 


been found to be unsatisfactory in two respects 


Che second and third volumes of the cone tables have 
They have been 
criticized because of their inconvenient tabulation and because 
are based is inadequate near the cone 
fa 


Empirical evidence is presented to show that 


on which they 
The 


transformation 


the theory 
coordinate 


surface former is eliminated by means 


the latter may be ignored in practice The exact nature of cer- 


tain numerical errors in the table is also pointed out 


INTRODUCTION 


— THREE VOLUMES of the M.I.T. tables’ * * for 
supersonic flow over right circular cones with at- 
tached shocks were originally prepared to provide 
a handy catalog of this important type of conical flow. 
The first volume, containing the results for unyawed 
cones, is tabulated so as to give easy access to the vari- 
ables in any part of the flow field. However, several 
that the 
ing cones, contained in the latter two volumes, are not 


authors! have noted results for yaw- 
so satisfactory. 

The two reasons for the unsatisfactory character of 
these two volumes of tables may be stated briefly as 
follows: 

(1) The tables are not tabulated in the most con- 
venient coordinate system for practical use. 

(2) The method of solution used to compute the 
tables is not valid near the surface of the cone. 

The first of these difficulties may be easily remedied 
by a change of coordinates as noted by Van Dyke, 
Young, and Siska.* This brief note was expanded by 
Young and Siska’ into a report giving formulas for the 
flow variables in a coordinate system which is still not 
the most desirable one, since one of the angles is meas- 
ured from the wind axis. This present report com- 
pletes the the their 


method A, which assumes that the intermediate sur- 


transformation along lines of 
faces are yawed with the solid cone and remain circular 
cones. However, their method B, although more com- 
plicated, has merit in that it renders the approximation 
for the complete flow uniformly valid, somewhat in the 
spirit of The 
present method, as does their method A, has the dis- 


Lighthill's straining of coordinates." 


Received November 11, 1958. Revised and received Febru- 
ary 17, 1954. 

* The authors wish to express their appreciation to the re- 
viewing committee for the detailed constructive criticism. 
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advantage that it is not appropriate for part of the 
region just behind the shock. However, at the shock 
and for most of the flow field, the present method is 
entirely appropriate and is the most simple and direct. 
He 


© used 


The second difficulty was first noticed by Ferri.’ 
pointed out that the form of the Stone solution? 
to compute the tables was not valid near the surface of 
the cone. Briefly, Ferri showed that the entropy must 
be constant on the cone, while the perturbation solu- 
tion of Stone gives a cosine variation in the entropy to 
the first order in the angle of yaw. He also gave a 
method for correcting the Stone solution near the cone. 

The present authors have found that, although these 
corrections are logically necessary, they may be neg- 
lected in practice, as will be seen later. This is espe- 
cially desirable since these corrections are difficult to 
apply in an unambiguous manner in the first order, and 
the second order corrections have never been derived. 

This paper is intended to give a descriptive review of 
the problems mentioned above together with a collec- 
the tables. 


Any one who is merely interested in using the formulas 


tion of formulas that facilitate the use of 
may turn immediately to the last section. 


DESCRIPTION OF THE PROBLEM 


The Stone theory® ® has been developed in terms of 
spherical polar coordinates (7, 6, ¢), with the axis @ = 0 
parallel to the flow at infinity. (These coordinates 
will hereafter be referred to as the wind coordinates. 
The axis of the cone has the equation 6 = «, ¢ = 0 
In the nonyaw case (e€ = 0), the axis of 
This means that 


(see Fig. 1). 
the cone and the axis 6 = 0 coincide. 
the equations for the cone and the shock are 6 = @, and 
§ = 0,, respectively. However when ¢ + 0 (in our coor- 
dinate system e€ is negative for positive angle of attack), 
then the cone axis does not lie parallel to the wind axis, 
and the flow quantities at @ = 4, are not generally those 
on the surface of the cone. Likewise the flow quanti- 
ties at 6 = @, are not in general those at the shock. 
In spite of this the range of the M.I.T. tables is @, < 
o =. G. 

The range of @ given in the tables arises from the fact 
that Stone satisfied his boundary conditions at 6, and 
6,, instead of at the true positions of the cone and shock. 
This procedure is perfectly consistent with the pertur- 
bation method used by Stone, but can cause some con- 


fusion when the tables are to be used. The situation 





for 
bh | 


the 


the 
an¢ 


wil 


eq! 
ot 
alc 

wh 


CoC 


at 
th 
an 


th: 


t the 


shock 
1Od 1s 
lirect. 
He 
used 
ice of 
must 
solu- 
py to 
ive a 
DORE. 
these 
neg- 
“Spe- 
It to 
and 
d. 
w of 
llec- 
bles. 


ulas 


s of 


ates 


es. 


5 of 
hat 
ind 


Or- 


X1S, 
ose 


iti- 


lA 


ict 


nd 


IT- 
n- 
on 














A GUIDE TO THE USE OF 
wa 
“ > as 
x, 4: A 
. pA 6. y h < 
> ¢ + } 4 
: SOD AD aT TT 
$ 2 +. YA 
- ’ ? 
[ L 
/ 
/ - =— 
¥ 
Fic. 1 Wind and body coordinate systems 


at the surface of the cone may be remedied by trans 
forming the flow variables from wind coordinates to 
body coordinates (7, @,, y 
When this is done, the equation 
The equation for 


, Where the polar axis is now 
the axis of the cone. 
9 #, is the equation for the cone. 
the shock is somewhat more complicated since the cone 
and the shock are not inclined at the same angle. This 
will be treated in a later section. 

‘he perturbation method already referred to uses the 
angle of attack as the perturbing element. That is, 
the solution is expanded in powers of the angle of attack 
and the coefficients are calculated by substituting the 
expansions into the exact equations of motion and 
equating like powers of «. Because of the assumption 
of conical flow these coefficients are functions of 6 and @ 
alone. Reference 1 gives the zero order coefficients 
while references 2 and 3 give the first and second order 


coeflicients, respectively. 


THE COORDINATE TRANSFORMATION 


If the coordinates of a point in body coordinates 


r,0,, w) are given and the flow quantities are desired 
at this point, one could compute to first or second order 
the coordinates (7, 
and then compute the flow quantities by the Stone- 


6, @) with respect to the wind axis 
Kopal equations. This procedure has the disadvantage 
that certain regions of the flow will be outside the range 
of the tables;* in particular, that section of the surface 
of the cone where 0 < W < 7/2 will be out of the 
range. 

The following procedure is much more satisfactory 
and will be developed here. Transform the Stone- 
Kopal equations to body coordinates and then compute 
the flow quantities using the coordinates of the point 
r,@,, WY) as arguments directly. 

This method has the advantages that the entire sur- 
face of the cone may be reached and that the formulas 
may be specialized to the surface of the cone, where the 
coefficients may be looked up with no interpolation; 


rhe tables could have been extended to cover these regions. 
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except between the tables for Mach Number, if de 
sired. 

It should then be noted that the velocity components 
computed will be the components in the wind coordi- 
nate system. Only the arguments of the components 
have so far been changed. The velocity components 
must then be changed to components in the body co 
ordinate system. 

The first and second order expansions will be dis 


cussed separate ly 


FIRST ORDER THEORY 


In the first order theory the Stone-Kopal formulas 


associated with the tables are the following 


ul ii + ex COS @ la 
_ eV COS OD lb 

uw ez sin @ le 
po- €n COS @ ld 

p p + €€ COS m le 


where “, v, w are the velocity components, p the pres 
sure, and p the density. The coefficients %, x, 7, y, 
s, p, n, p, and & are functions of @ alone, and the bars 
indicate the nonyaw quantities. 

The above formulas refer to the wind coordinate 
system, and the velocity components are the com 
ponents in this system. The quantities of general 
interest are those in which the axis of the cone itself 
is the central reference line, as is pointed out by Van 
Dyke, Young, and Siska. 

To the first order in ¢, the transformation between 
the coordinate systems is / r, 6 6, + ecos yp, @ 
ware the velocity com 


Il’ the 


y—ecotésny. Nowil uv, 2, 
ponents with respect to the wind axis and uw, I’, 
components with respect to the body axis (the # com 


ponent is the same in both systems), then Eqs. (1) be 


come 
u(6,, y = Ur [x = % le cos y Da 
‘ / 
1(6,, ¥ 
I + cot @ 
7+ \ Wu jecos y 2D 
WV(6,, ¥ Z ese 6,) € sin ¥ X 
pi + 3 cot 6 
p é., y = p — Bus € COS y 2d 
; 
p cot 6 
pig ¥ =P Ts € cos y Ye 
a* - 
where @ is the nonyaw velocity of sound. Here the 


coefficients are to be evaluated from the tables, using 
for 6 the angle @, measured from the body axis. 

On the surface of the solid cone 6, = 4,, i 0, and 
y = 2%. In this case Eqs. (2) reduce to Eqs. (1) with 
y replaced by ¢ and (1b) replaced by 
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y= 0 (1b’ To second order the relation between the coordinate 

. a ; systems is as follows (see also reference 6): 

If the flow quantities are desired at the shock on the ~~ 

downstream side, a slight modification of the Eqs. (2), r = 7 

6= 6, + ecosp + !/2e? cot 6, sin? ¥ | p 

y — ecot 6, sin y + 
1/,e7(1 + 2 cot*@,) sin 2y | 


as pointed out in reference 4, is necessary. In each 
bracket, each term except the first should be multiplied ¢ = 
6/€, where 6 is the yaw of the shock cone, ex- 


by a = 
cept in the formula for 1’. For example, at the shock — ,, are ae P : 
Thus if fis any analytic function of 6, the value of f in 
wave, heats: d 
body coordinates may be computed, to second order, by 


u=u+ [x + ale cos y (3) the formula: 
f(0) = f(0,) + ef’(6,) cos W + 


with 6 = @,. Since 6 and « are of like sign, @ is a posi- 1/,€°[f'(0,) cot 6, sin? wy + f"(6,) cos*y] (Ga 


where the coefficients are to be taken from the tables 


tive quantity. IV is independent of a in the first order. . : 
peal: aa é f : : lhe following second order approximations also hold 
rhis coordinate transformation provides, in a sense, 


an extrapolation to certain regions neglected in the ecos@ = ecosy + e’ cot 6, sin’ y | 

computations. This neglect was due to the simplifica- € sin @ e sin Y — 1!/s€? cot 6, sin 2y 61 

tion, mentioned previously, of shifting the boundary e? cos 26 = e&* cos 2y " 
<* sin 2 = e* sin 2y 


conditions. 
The following formulas relate the velocity compo- 


SECOND ORDER THEORY ri , 
7 nents to second order. (The radial component is the 
When the yaw of the cone is larger, and second order same in both systems.) (See also reference 6 
approximations are desired, the formulas become con- , eee nia 
PI : ' V =v(1 — !/2€? csc? 6, sin? y) + 
siderably more complicated. In wind coordinates the . 
R iy = we csc @, sin y i 
formulas are as follows: “ae c 
I] = w — vie csc 6, sin y = | 
1 ‘»€* ese 6, cot 6, sin 2y) 


u= tu + ex cos @ + E(u + Ue COs 2p) 
v =od+ eycos¢ + e7(% + %% cos 26 


: : Note that since w is itself a first order quantity the 
w= essin @ + EW, sin 2¢ (4 ‘ : cee a ; 

® change in v is entirely second order and the first order 
p= pt+encos ¢ + €*(po + pe cos 29) . 
p = pt e€ cos gd + €(po + p2 cos 2g) 


change ts nil. 
Eqs. (4) transformed by Eqs. (6) and (7) are listed 


where 1, X, Mo, M2, ete., are functions of 6. below. 


u= t+ ex + 8) cosh + (my + 14a’ + 1/8 cot 0, + 1 /oxy + 1/ox cot 6,) + 
€°(u2 + 1/40" — 1/45 cot 0, + 1/2v — '/ox cot 6,) cos 2y (8 


v= i+ ey + 8’) cos + e(19 + 1/40" + 1/48’ cot 6, + '/oy’ + '/2y cot 6,) + 
€°(v2 + 1/40" — 1/40’ cot 0, + '/oy’ — '/sy cot 8,) cos 2y (9 


V= a+ ey + 0’) cos P + ev + 1/40" + 1/48’ cot 0, — 1/40 csc? 6, + 1/oy’ + 1/oy cot 0, + 
; 92 CSC 6, ) + €7(v2 + 1 io” —! Pid cot 6, a l 4v esc” 6, 4. 1 oy’ —! 2V cot 6, — 1 23 CSC 4 cos 2y (10 
w= essin W + e2(we + !/o2’ — 1/22 cot 6,) sin 2p (11 
W = es — dese O,) sin Y + €?{ we + 3/22’ — 1/22 cot 6, + '/2[8 cot 0, — (y + 8’)] csc 6,} sin 2y (12 


p/p = 1+ [(n/p) + (p'/p)] € cos Y + &[(po/ Pp) + (n/ 2p) cot 6, + (n'/2p) + (p’/ 4p) cot 6, + (p"/4p)] + 
e*[(po/p) — (n/2p) cot 0, + (n’/2p) — (p'/4p) cot 6, + (p"/4p)] cos 2y 
p= 1+ [(é/p) + (p'/p) Je cos W + &[(p0 p) + (E/2H) cot 0, + (E’/26) + (p'/4p) cot 6, + (p"/4p)] + 


€*[(po/p) — (E/2p) cot 6, + (£'/25) — (p'/4p) cot 6, + (p"/4p)] cos 2y (14 


where the quantities are to be evaluated from the tables using for @ the value 6, measured from the body axis, and 


(see also reference 9) 


o’ = —{[a°(a + d cot 6,)]/(a@ — w)! -— a (15) 


a” = —o'[cot 6, + (p’/yp)] + a[cot? 6, — 1 + (1/y)(p’/p)? + (p"/yp)] | 


p'/p = [yvi(a + é cot 6,)]/(a? — 8) 





for 


wl 


wl 


O1 
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p"/b = (a? — yo*)/[y(a*® — 3) ]} (b’/p)? + ('/8)- (a? + ®), (a d*)]-(p’/p) + 
[yi cot 6,(3’ — 3 cot 6,)|/(a? — 7 18 
y’ = —{ [cot 6, + (p’ yb) ly + 2x + 2csc 0, + (3/4 [(n’/ pb) — (n/p): (p'/p)]} 19 
n’ | Pe: ,n| Pp’ , p’ { 
= — Ee°- — 3 - yu +o) + yex + OL ¥ Cot 6, + + Zyv csc 6, 20 
p a? — ol p Pj p P 
it san = , ; ) Se ; 
'/p = (] ¥ ' (€/p) — (n/P)\ (PP /P) +r (nn /P)i; P/pP = P/YP 21) 
=—V - > 
"p= (1 - ¥)/Y¥7] (p’ Pp)? + (i/y) (p” P) (22 
s’ = (a°/7a)-(n/ p) ese 0, — [(a4/%) + cot 6,]z (0, # 0, 23 


These formulas are not strictly accurate, due to the change in coordinate system. However when used in the 
formulas for the flow variables all discrepancies fortunately become third order or higher. 
The following alternative formulas for y’ and #” also hold 


a” = —(6’A + 6B (16a) 
where 


A = cot 6, + A[2a + 32 cot 6, + (y + 1)Ad(a + 6 cot 8,)] 
B = 1 — cot? 6, + A[A(y — 1)a(% + 6 cot 6,) — 6 cot? 6,] 


A =6/( — 9°) 
y’ = —(Ay + Bx + Cd) (19a) 
where 
Cd = [(1 + do) /(sin? 6,)] (x + 2 sin 6,) 


On the surface of the solid cone 


6, = 6,6 = 0, y = 24, 6’ = — 2, p’ = 0 
p’ = QO, p” p= —2(u/a)*, n’ aoa 2pu- 
p" = —2pi?, o” = 20 cot 4, &' p = 20/4)? 
y’ = —2% cot 6, — 2x — zcsc O,, % = xX 
Vv. = ucoté, + zcscé, + x 
z’'* = —2u csc 0, — 2 cot 8, 
Therefore 
u(6,, ¥) = A+ excosy + e(uy + '/o% + '/ ox cot A.) + (uy + 1/0 — '/ox cot 8,) cos 2p (24) 
v(6,, W) = —e*s csc 8, sin? ¥ (25) 
V(6,,¥) = 0 (26) 
w0,, W) = essin Y + e(w. — a csc 0, — z cot 6,) sin 2y 27) 
W(0,, W) = es sin yy + e?(w. — a esc 0, — 2 cot 6,) sin 2y 28 


[p(6,, ¥)]/[p(@.)] = 1 + e(n/p) cos W + e?[(po/p) + (yit?/2a?) + (9/ 2p) cot 4.) + 
e[(pe/p) + (yit?/2a") — (n,/ 2p) cot @,] cos 2y (29 


[o(6,, w)]/[p(@.)] = 1 + &(E/p) cos PW + €[(po/p) + (a7/2a) + (&/ 2h) cot A,] + 
e*[( po p) + (a 2a-) — (& 2p) cot 6 | cos 2y 30 


Note that on the cone w = VW to second order. 


UsE AND VALIDITY OF THE THEORY trate this by indicating the methods of computing the 


The transformed formulas of the previous section PTESSUre coefficient C, on the surface of the cone, where 


may now be easily used to compute the flow variables C, = (2/yM?) [(p/)-(b/p:) — 1 


at points between the shock and the cone. We illus- 
In the above equation .\/ is the free stream Mach Num 


* Using L’Hospital’s rule, as in reference 9. ber and the subscript | refers to the free stream. 
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o EXPERIMENTAL DATA 
—— COMPUTED CURVE REF. 7 


e COMPUTED POINTS, 
PRESENT METHOD 


0, 7.5° 


M=1.6 


oO 
ro) 





PRESSURE 
COEFFICIENT 
° 
°o 


fe) 
180 150 120 90 60 30 ie) 
COORDINATE W 
€=-6° 





180 I50 120 90 60 30 (@) 
COORDINATE W 
e€=~-l2° 
Fic. 2. Pressure distribution around the cone at WM = 1.6. 


C, may be computed in either of two ways for both 
first and second order theory. We present both ways 
here and give only the recommended method in the 
next section. These methods consist of (a) computing 
pb, p to the first or second order using the tables and pre- 
vious formulas and computing /, p; exactly to give p/p, 
or (b) computing the square of the speed g° using the 
tables and formulas, and obtaining p, p; from 


b/fr = TU — @?)/0 -— ww). (31) 
where 
y+ 1\tMo- ir 1 
nl = € ) —— x 
: i — 
1 oy 
1 We) 
ee vy— ] 
(vir — : ) 
and 
M? = M? sin? 6, (33) 


It has been pointed out’ that when using (b) to first 
order better agreement with experimental data is ob- 
tained if the second order terms in g? = u*? + V? + W? 
are retained, i.e., if wu, 1’, and IW are computed to first 


order before squaring. 


MAY, 1954 


So far we have tacitly assumed that the Stone theory, 
and as a consequence the last two volumes of the 
M.I.T. tables, are mathematically correct. 
Ferri’ has shown that the entropy distribution given by 
Stone is not correct at the surface of the cone. In fact 
Ferri then gives 


However 


the entropy must be constant there. 
some corrections to the Stone theory which are de- 
signed to make it consistent with his general conclu- 
sions. 

The actual method of applying these corrections is 
not clearly stated, so that we have been unable to re- 
produce his results. It turns out that the present 
theory may be applied without these corrections to give 
good results, as we shall show. 

We now present some empirical evidence to support 
our thesis that the M.I.T. tables may still be used for 
practical results in spite of any theoretical objections 
to the Stone theory. In the following calculations we 
have used method (a) above, where p/p has been found 
to second order. 

Figs. 2 and 3 show plots of the pressure coefficient 
on the cone surface against y. Also appearing in the 
graphs are Ferri’s corrected results, which are first 
order, and some experimental data taken from refer- 
ences 7 and 8. The computations carried out by our 
method agree with the experimental data extremely 
well in all cases with the possible exception of the last 
graph. In making the computations for 1/ = 6.86 
it was necessary to use the fact that the pressure co- 
efficient for hypersonic flow is essentially independent 
of the Mach Number. The data in the third volume of 
the M.I.T. tables go up only to JJ = 4+, which was 
the Mach Number used in our calculations. This may 
account for some of the discrepancy in the last graph. 

It is really surprising that the theory gives such ex- 
cellent results for angles of yaw as large as 14°. This 
is certainly better than would be expected from any 
preliminary analysis. 

Our recommendation is that for a yaw less than 
about 6° method (b) be used to first order and that for 
a yaw greater than this the second order form of 
method (a) be used. Method (a) is more direct and is 
to be preferred for second order computation. How- 
ever, using method (b) for first order permits the reten- 
tion of the important second order terms (mentioned 
previously) without the use of any second order theory. 
No really sharp limits on the yaw angle can be set for the 
first order theory to hold since it will vary slightly with 
Mach Number, cone angle, etc. In all cases it would 
seem that the Ferri corrections may be neglected. 


TABULAR ERRORS 


Although many precautions were taken to avoid 
errors in computing the M.I.T. cone tables the second 
volume does contain a systematic error.'? The tabu- 


lated values of z are incorrect. The error comes from 
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© EXPERIMENTAL DATA FOR TWO 
DIFFERENT VALUES OF r 


——COMPUTED CURVE REF 7 


COMPUTED POINTS, 
PRESENT METHOD 


G,= 10° 


M = 6.86 


0.24, 
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COORDINATE y 
e€=-14° 
Fic. 3. Pressure distribution around the cone at J = 6.86. 


using the wrong sign for C in computing from Eq. (54 
of that report. The correct value of z is found as fol- 


lows: 


(correct z) = —(s from tables) — (2x, sin 6 


Ferri, Ness, and Kaplite’’ suggested that the error 
just mentioned might have carried over into the third 


volume of the tables. By performing an indirect check 


Sec ond Order. 
in the tables at 6 = @,. 


TO THE USE OF 


Use the formulas stated below together with formulas (15)—(2 
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it is easy to see that this is not the case.* Thus it ap 
pears that the only correction that need be made is one 
for 2. 

PROCEDURE FOR COMPUTATIONS 


(A) On the Surface of the Cone 


First Order. 
ties. Look up a, x, 2,9 pandé pata value of @ equal to 


Use Eqs. (la)-(le) for the flow quanti 
the semi-apex angle 6, of the solid cone. Use for € the 
value of the angle of yaw in radians. If the pressure 
for / , keeping 


Note also 


coefficient is desired, use formula (31 
the second order terms which appear in q 
thatv = 0. 

Second Order.—Use 
and (30) looking up all quantities at 6 


formulas (24), (26), (28), (29), 
#,in the tables. 
If the pressure coefficient is desired, use formula (29) to 


compute p, p, then compute pp; exactly, etc 


(B) Between the Shock Wave and the Solid Cone 


In general, the flow quantities will be desired at points 
around the cone at a uniform distance from the surface 
of the cone. Therefore bi dy coordinates will be used. 

First Order. 
quantities in the table using for @ the angle measured 


Use formulas (2a)—(2e) looking up all 


from the axis of the solid cone. If the pressure coeffi 
cient is desired, see the suggestion made for the surface 
of the cone. 

Second Order.—-Use formulas (S), (10), (12), (13), 
(14), and formulas (15)—(23) looking up all quantities 
in the tables using for 6 the angle measured from the 


axis of the solid cone. 


(C) At the Shock Wave on the Downstream Side 


First Order.—-Modify formulas (2a), (2b), (2d), and 
(2e) by multiplying each term in the coefficient of €, 
except the first, by a = 6 «. The value of a may be 
obtained from the tables. Then use the modified 
formulas to compute the flow quantities, with @, = @,, 
and looking up all quantities in the tables using for 
6 the angle @,, as given in the tables. Use formula (2c 
without the above modification. These formulas are, 
of course, the terms remaining in the second order 
formulas (stated below) after neglecting terms with 
F: 

* In a private communication Prof. Kopal has confirmed the 
fact that the correct values of s were used in computing the third 


volume of the tables 


(23) with all quantities looked up 


1 -/ 


u= i + ex + a@) cosy + € luo + '/ex cot 0, + '/2ay + (Bo + '/2@ cot 6,)% + * sad’ | 
e[u2 — ox COt A, + I/oay + (Bs — |} 2a cot 6,)3 + 1/40 é’] cos 2y 34 
"=a t+ ey + ad’) cosy + [mm + (1/2)y cot 0, + '/say’ + (Bo + '/2@ cot O,)5’ + 1/40" — 
9 csc? 0, + 1 /o2 esc Oy] + [ve — (1/2)y¥ cot 0. + 1/2 ay’ + (B2 — '/2@ cot 6 a’ + lM /ya*’e” + 
1/45 csc? 6, — 1/28 esc 6,| cos 2¥ 1) 
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1 + el(n/ Pp) + al(p’/p)| cosy + Ee} (po p) + [A 
[Bo + (a/2) cot 0] (b’/p) + (a? 


(a/2)-(n’/ p) + [pe 


B) + [a(p’,p)|jeos Y + €*} (p/p) + [! 
[Bo + (a2) cot 6.] (p'/ 5) + 


(a 2)-(é Pp) 
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Thermal and Creep Effects in 
W ork-Hardening Elastic-Plastic Solids’ 


ALEXANDER J. WANG? 


eesti PRINCIPLES GOVERNING the isothermal 
deformation of a work-hardening elastic-plastic 
solid have been given by Hodge and Prager! and 
Hill. In the present note it is shown how these prin- 
ciples can be extended to include thermal and creep 
effects. 

Using rectangular Cartesian coordinates x; (7 
|, 2, 3), denote the infinitesima! displacement from the 


standard state by u,;, the infinitesimal strain by e€;;, the 
stress by o;;, the mean normal stress by o, and the stress 


deviation by s The mean normal stress is defined as 


3) Ci; (1) 


where the usual summation convention regarding re- 
peated subscripts is used; the stress deviation is defined 
as 


Sig = Oy — 006i; (2) 


ij 


6,, being the Kronecker delta. The mean normal stress 
is a invariant of the stress tensor. 


useful in the theory of isotropic plastic solids are 


Further invariants 


Jo = (1/2) siti, J3 = (1/3) SijSjeSq i (3) 

To specify the mechanical behavior of the isotropic 
solid, consider an element that, at the generic instant 
', has the temperature @ and the strain ¢;; and is under 
the stress o;;.. In the interval between the instants 
tand ¢ + dt, let the temperature be changed by dé@ and 
the stress by do;; = ds;; + do6;). 

The corresponding change of strain, de;;, will then be 
assumed to consist of the following components: 


(1) The e/astic component 


de; = a(@) ds;; + B(A)do6; (4) 


where a(@) is one-half of the reciprocal of the temper- 
ature-dependent shear modulus, and @(@) is one-third 
of the reciprocal of the temperature-dependent bulk 
modulus. 

(2) The thermal component 


de:” = [a’(0@)s;; + B'(0)06,;| dd + y(9)6;; dO (5) 


where the prime denotes differentiation with respect to 
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n us rersily 


the temperature and y(@) is the coefficient of linear 
thermal expansion at zero stress. 


(3) The creep component 
de = P(6, Fe, Js (Og Od dt 6 


where ¢ = ¢(0,J2,J3). 
(4) The plastic component 


\0 if P(0,J2,J3) < h 
W(0,/s,/3) (Oy Oc; ) (dy _ dy if 


V(6,Jo,J3) = h (7) 


where the scalar / describes the state of hardening of the 
considered element at the time f. 

Eqs. (4) and (5) result from a generalized form of 
Hooke’s law in which the elastic constants are functions 
of the temperature. 

In Eq. (6) the derivative O¢ O¢;, must be evaluated 
for the state of stress existing at the time /. The ex- 
pression ®0¢ Ooa;; then represents the rate of secondary 
creep corresponding to this state of stress. This ex- 
pression is sufficiently general to comprise all isotropic 
laws of secondary creep that have so far been proposed 
in the literature on creep. 

In Eq. (7), the equation ¥(6,/2,/;) = h represents the 
temperature-dependent yield limit for the state of 
hardening achieved at the instant ¢. The first line of 
Eg. (7) therefore states that there will be no change in 
plastic strain where the state of stress is below the yield 
limit. The differential dy in the second line of Eq. (7) 
must be evaluated from the given temperature and 
stress at the instant ¢ and the changes in temperature 
and stress during the considered time interval. On 
account of the absolute value in the second line of Eq. 
(7), there will be no change in plastic strain even if the 
state of stress is at the yield limit provided that dy < 0. 
A special case of the combination of creep and plastic 
flow described by Eqs. (6) and (7) has been investi- 
gated by Odgvist.* 

The extremum principles that are to be established 
Con- 
sider a mass of work-hardening plastic material that 
has been deformed and, at the time /, occupies a region 
Suppose that the tem- 


concern the following boundary value problem. 


|’ bounded by the surface S. 
perature 6, the stress o;;, and the state of hardening h 
are known throughout lV. If the unit vector along the 
exterior normal of S is denoted by n,, the surface trac- 
tion in the considered state is 7; = oijn;. Prescribe 


now infinitesimal changes dé of the temperature 
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throughout J’, infinitesimal changes d7); of the surface 
traction on the portion S; of the surface and infinitesi- 
mal displacements dl’; on the remainder S, of the 
surface. What are the corresponding changes of stress 
do;; and the corresponding displacement du; through- 
out |"? 

The infinitesimal displacement du; causes the strain 


to change by 
de;, = (1/2) [(0/Ox,;) (du;) + (0/Ox;) (du,)] (8) 


The change of stress must satisfy the equation of equi- 
librium which, to within higher order terms, can be 


written as 


(O/Ox;) (de.;) = O (9) 


Finally, the following boundary conditions must be 


satisfied 
doin, = dT;,on S7 (10) 


du, = dU,on S; (11) 


The problem thus consists in determining the change 
of stress do;; and the displacement du; in such a manner 
that the boundary conditions, Eqs. (10) and (11), and 
the equation of equilibrium, Eq. (9), are satisfied, and 
that the strain change computed from Eq. (8) is related 
to the given values of 6, d@, o;; and h and to the sought 
stress change do;; by means of Eqs. (4) through (7). 
The first extremum principle compares the actual 
changes of stress and strain, do;; and de;;, to a fictitious 


* and the corresponding change of 


nf 


change of stress do; 


strain de,;*. The stress change do;;* is supposed to 


satisfy the equation of equilibrium and the boundary 
condition on S;, and the strain change de;;* is asso- 
ciated with 6, d@, o;;, h, and do;;* by means of Eqs. (4) 
through (7), but need not be derivable from a displace- 
ment field. The principle of virtual work then furnishes 


the equation 


S \(do * — dg;;) de;;| dV = 
ST? — dT;)dU;|dS¢_ (12) 


where d7;* = do;;* n,. 


The integrand of the left-hand side of Eq. (12) can be 
transformed as follows: 


2(da;;* — do;;) de;; = (do;;*de;;* — da;,deé;;) 
(da j;*(de,;* — dej;) + de; (do;; — do;;*)| (13) 


On the right-hand side of this equation, introduce the 
components, Eqs. (4) through (7), of the change of 


e “° a 4 « é 
strain. Since de,;*" = de," and de;;*° = de;;‘, 


2(do;;* — dai) de; = doj;*(de;;* + dei;*? + 
de,;*°) — doi(de,; + de,’ + de;;°) — 
[da j;*(dej;*° + de,;*? — de; — de;;?) + 
(dei;° + de;;”) (doi; — do;;*)] (14) 


The bracket in Eq. (14) involves only elastic and plastic 


changes of strain. When the case do;;* = do; is ex- 
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cluded, it can be shown that this bracket is positive. 
This is done in exactly the same manner as in the proof 
of the extremum principle for isothermal deformation 
of a work-hardening plastic material (see Hill,’ pp 
63-64). 
than the first two terms on the right-hand side unless 


Thus, the left-hand side of Eq. (14) is smaller 


do;;* = do;;.. When this result is introduced into Eq. 


(12), the following relation is obtained: 


(1/2) S (do *(de.,* + de,” + de,;* )J dV a 
S (dT *dU;,) dS; > (1/2) S |do (de;; + de 
de;,)|dV — S (dT dU ;) dSv 15 


where the equality sign holds only if do,;* =do;;.. The 
relation (15) establishes a minimum property of the 
actual change of state. 

Before the second extremum principle can be dis- 
cussed, it must be shown that Eqs. (4) through (7 
can be transformed so as to represent do;; as function 
of de;;, We note first that de,” and de 
diately from the data of the considered boundary value 


follow imme- 
problem. Since these two components of the change 
of strain are known, giving de;; is equivalent to giving 
de,’ + de,’ = de, — de; — de 
the sum de;;‘ 


The proof that 


+ de,” specifies a unique do;, then pro- 


ceeds exactly as in the case where de;,/ + de,” repre 
sents the entire change of strain (see Hill,” pp. 68-69). 

The second extremum principle compares the actual 
changes of strain and stress, de;; and do;,, to a fictitious 
* and the corresponding change of 
* 


change of strain de 
a is supposed to be 
* 


stress do; The strain change de 


derivable from a displacement field du;* that satisfies 
the boundary conditions on S,; the stress change do 

is associated with 6, d6, o;;, h, and de;;* by means of 
Eqs. (4) through (7), but need not satisfy the equation 
of equilibrium or the boundary condition on S;. The 
principle of virtual work then furnishes the equation 


S (do; (de;;* — de, \} dV = 
S (dT (du,;* — du;)| dS; (16) 
The integrand on the left-hand side of Eq. (16) equals 


da;;(de,;** + de,;*? — de; — de;;’ (17 


and does not involve the thermal or creep effects; it 
can therefore be transformed in exactly the same 
manner as in the case where the sum of elastic and plas- 
tic strain changes represents the total strain change 
(see Hill,” pp. 65-66). As a result of this transforma- 
tion, the following relation is obtained: 


S (dT du;*) dS; — (1 2) S (do *(de;,*° + 
de,;*?)|dV < Sf (dT idu;) dS;p — 
(1 2) S [do;;(de;,;' Py dV (1S) 


+ de 


where the equality sign holds only if de;;* = de;;. The 
relation (18) establishes a maximum property of the 
actual change of state. 


(Continued on page 360 
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RIEF REPORTS of investigations in the aeronautical sciences and discu ms Of papers pub 

lished in the JOURNAL are presented in this special 

completed as soon as possible after receipt of the materia The Editorial Committee d , 
itself responsible for the opinions expressed by the correspondent 


15) 
The ° 
the 
The Compressible Laminar Boundary Layer J is « particular solution of the energy equation, Eq. (4). Eq. (9 
dis- on a Rotating Body of Revolution follows 
1 (7 t (u?/2) + 7/2) + Rwri constant ie) 
. Sheng To Chu and A. N. Tifford 
aaa Department of Aeronautical Engineering, The Ohio State University For Eq. (9) to be valid at both edges of the boundary layer, 
ime- Columbus, Ohio u,2/2) + must equal [(1/2) + k|w?r? 4 ind be identically 
alue anuary 27, 1954 the integration constant In other words, Eq 4%) applies, in 
inge general, only when the surface enthalpy varies parabolicly with 
ers RELATIVELY GENERAL analysis of the ‘“‘constant property” Phe arbitrary constant, k, may be adjusted, however, so as to 
ei A laminar boundary layer on a rotating body of revolution in allow a variety of parabolic representations of the surface enthalpy 
that ‘xial motion has been given in reference 1 As a sequel, consti distribution to be covered Additional interesting points worth 
pro- tuting a portion of his doctoral research, the senior author is noting ar¢ 
‘pre undertaking the investigation of the detailed effect therein of the a Eq 9%) reduces to the well-known Crocco integral (refer 
69), variability of fluid properties. The initial result obtained is a ence 3) when the body does not rotat« 
tual Crocco-like energy relation applicable when the Prandtl Number b) The analogy of reference 1 between the frictional torque 
: is unity. By means of this relation, the compressible laminar and the surface rate of heat transfer is readily shown to apply 
_— boundary layer on a rotating body of revolution may be correlated despite the variability of the fluid properties 
e ol with a fictitious incompressible laminar boundary layer under ¢) An insulated rotating surface acquires the surface enthalpy 
» be more general circumstances than those reported by Illingworth distribution, = hy + (wr)?/2, 1e., k l Conversely, a 
sfies reference 2 Furthermore, the analogy between the frictional rotating body having a surface enthalpy (temperature) corr 
a torque and the surface rate of heat transfer of reference 1 is sponding to the free stream tot ul enthalpy temperature k 
: readily seen to be unaffected by the variability of fluid properties — 1/2—does experience heat transfer 
Ss ol , A ‘ : : ‘ : 
: The steady compressible boundary layer equations for a ro Aside from its general character, Eq. (9) may be used ad 
lion tating body of revolution are vantageously in obtaining a correlation between the compressible 
The ee ee ( boundary layer on a rotating body of revolution having a surface 
ion - 4 ‘ pressure gradient and an associated fictitious incompressibl 
plu(Ou/Ox) + v(Ou/dy) — (w?/r) (dr/dx = boundary layer by a method somewhat analogous to that em 
pili(du;/dx) + (0/Oy) [u(Ou/oy (2 ployed by Illingworth (reference 2 However, neither his 
16) pf uld(rw)/dx] + vf[d(rw)/dy]} = (2/dy) fula(rw)/dy]} (3) restrictive condition of constant enthalpy inside the boundary 
. . ‘ layer nor his limitation of the local free stream Mach number 
tals piu(dr/Ox) + v(O1/Oy tT Up ui(du/dx) = variation along the surface to 1J; = (ar + 8 1/2 is required any 
(1/o) (0/dy)[u(di/dy)] + p[(Qu/dy)? + (Qw/dy)?] (4) longer 
17 (y — 1)/yloi = p = fi = [Cy — 1)/ylou (5 Assuming that the coefficient of absolute viscosity is propor 


; ; ; ’ t ait tional to the absolute temperature, Eqs. (2) and (3) may be trans 
it Here x is the distance along the surface of the body in a meridian formed into Eqs. (10) and (11), 


plane, y is the distance measured normal to the surface, r is the 


me ; eget Peer oe ee sitet 
distance of the surface from the axis of symmetry measured 1/1 (du, /d 2s ‘ u(Ou/o 
las- normal to the latter, « and v are the velocity components along (w?/ mi? )r(dr/ds i? = u(0/Oy) [W( OU /O¥ (10 
nge the x and y directions, respectively, w is the rotational velocity and 
na- component, and the other quantities are denoted by conventional ae ‘ 2h 
on : ; ‘ Ol r-u Os = 0/OY)U|O( ru Oy; (11 
symbols. The energy equation in the free stream is 
where the stream function, y, satisfied the relations, 
1) + (1/2)u,2 = hy = constant (6 
rpu = Oy /Oy and rpv = —(Oy/Ox 12 
It is assumed that the Prandtl number, o, of the fluid is unity, : 13 
- » e u = u/u, w= W/owr,t = ii/i = » p o 
re and the enthalpy, 7, is a function of H alone, where wal 
) 
, m and 
H = (u?/2) + (w?/2) + kRwru (7) > 
. z 
rhe a { ] 
° . . . s = ie 0 4;aX 14) 
th w is the angular velocity of the rotating body, and & is an abri J ai 
1e " + . ‘ ‘ 
trary c ant. z ass ition and Exc (2), . . : 5 . : . 
_ constant. Combining this assumption c Ee Essentially this transformation is a combination of that of von 
(5), and ( i ay be Ww ‘ : ” , _ : . 
3), and (4), it may be shown that Mises (reference 4+) and that of Mangler (reference 5 From Eq. 
di/dH = —-1 (8) (8) we can write, 
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(a? — 7) = (a? — 1) j}1 + 1 /2(,2/i;)| } + 
1 /2(w?r?/1,) |? + (Rw?r?/i, od (15 
Substituting Eq. (15) into Eq. (10), we have 
(du,/ds) (%? — 1)| + f(O%/Os) — 
w?/u,2) [r(dr/ds)|w? + [(1 + b)(1/m,)| (du,/ds) X 
[1/2(w?r?/hy tt + (kw?r?/hy i] = 0(0/0Ww)[u( On /oY 


(1 + by 1/u 


{ 


where b = 17/21 The term dr/ds may be expressed as 


dr/ds = (dr/db)(db/ds) = (1 + b)?[(1/h,)(du,/ds 


Thus the above equation becomes 


(1 + b)(1/u;)| [(du,/ds)(a@? — 1)| + a(ou/dos) — 
(w?/hy) [r(dr/db)| w2(1 + b)?[(1/u1)(diy/ds)| + 
(1 + bY(1/m,)(duy /ds)| [1/2(w?r? /hy ee? + 


(kw?r?/hy it] = G(0/Ow)[a@(Oa/dy)| (16) 
and Eq. (11) retains its original form 
O( rh) /Os = (0/OP~)} a[O( 127) /OY] | (17) 
The boundary conditions are 
at y = 0, a = 0, we = 1) 
P (18) 
at y—> om, v=}, b= 0} 


It will be found possible to establish a correspondence between 
the compressible boundary layer problem represented by the 
differential system, Eqs. (16), (17), and (18), for the rotating 
body of revolution having a surface pressure gradient and an 
associated fictitious incompressible boundary layer problem. 

Consider the following system of equations 


oU/OX + OV/OY = 0 (19) 


U(O0U/OX) + V(OU/OY) — (W?/R)\(dR/dX) X 
[1 + (1/2)(U,/U,)?] + U\(dU,/dX)[((1/2)(W/U,)? + 


(RORW/U,?) — 1] = Nx(0?U/O¥?2) (20) 
and 
U[O(RW)/OX] + V[O(RW)/OY] = Nelo RW)/OY?] (21) 


Here the capitalized letters denote corresponding quantities in the 
fictitious incompressible boundary layer problem, U’, is a constant 
reference velocity, and Nx is the uniform kinematic viscosity. 
Note that when W = 0, these equations reduce to that of a two 
dimensional boundary layer flow of an impressible fluid about 
a cylinder. By changing to s and y coordinates, Eqs. (19), 
(20), and (21) may be transformed into Eqs. (22) and (23) 


(1/U,)[(dU,/ds)(a? — 1)]| + aouw/os) — 
(Q?2/U,2)R(dR/dB)w?[1 + 1/2(U,/U,)?|(1/U, (dU /ds) + 


(1/U,)(dU,/ds)}(QR/U,)2[(1/2)e? + kw)| { = a(0/dy) X 
[w(ouw/Oy)| (22) 
O(R27)/Os = (0/Ow)} a[O( R27) /dy] | (23) 


The fellowing definitions and nondimensional quantities have 
been used in obtaining the above differential equations: 


LPsU = dy/OoY, LPsV = — (d%/dX) (24) 


om 
s= L?M, Ps | UdX (25) 


where Ps is the uniform density; Ws, the uniform viscosity; 
and L, the characteristic length of the incompressible problem 
After comparing Eqs. (22) and (23) with Eqs. (16) and (17), we 


choose 
(1/r)(dr/ds) = (1/R)(dR/ds) (26) 


and 
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(1 + b)(1/u,)(du,/ds) = (1/U,)(dU,/d 97 


The integrals of Eqs. (26) and (27) are, respectively, 


rs) = aR a 
and 
(117/21, )(s) = U,%(s 9 
where c; and c are integration constants Furthermor we 
choose 
Q=oe ¢ = 1, ao =1/U,? = l/h 30 
The last relation gives b(s) = B(s). It is therefore obvious that 


with these choices the two sets of equations, Eqs. (16), (17) and 
Eqs. (22), (283) are now identical to each other Besides, th 
boundary conditions for both sets are the same, as given in Eq 
(18). 

The correlation mentioned previously is thus established 
That is, under the conditions given in Eqs. (28), (29), and (30), 
the nondimensional functions @ and # may be found from the 


incompressible flow Eqs. (19), (20), and (21) instead of the more 
complicated Eqs. (1), (2), and (3). 
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Contrclled Motion with Time Lag in 
Autopilot 


G. M. Andrew 
North American Aviation, Inc., Downey, Calif 
December 31, 1954 


_ EQUATIONS developed below apply to a second order 
system with arbitrary initial conditions, controlled by the 
autopilot, having an average time lag 7. The method illus- 
trated on one degree of freedom can be readily adapted for two 
or three degrees of freedom. Denote the natural frequency and 
the damping ratio of the uncontrolled system by w, and ¢ and 
the forcing function by 6(¢) then the equation of motion can be 


written. 
O(t) + 2erw,O(t) + w,20(t) = M5d(t l 


Assume that the control surface deflection having the time lag 7 


is proportional to displacement @(t) and rate @(t). Then the 
control surface defiection is represented by 
6(t) = 6) + KeAt — T) + K0(t —T7T 2 


where generally constants Ag and Kg are negative 
Substituting Eq. (2) into Eq. (1) gives an operational equation 
(in the time domain) of controlled system 


(D2 + 2%w,D + w,?)0(t) = 
Msi0 + MiKo + KeD)t — T) (3 


For the initial conditions (at ¢ = 0) 6); @) and 6, the Laplace 


Transformation of Eq. (3) yields 
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+ 2fu + wp7)0 = s0y + Oo + Zon 4+ VM 560)/s] + [.\M5(Ko + Kos )0 — \5K 9 { 
For technical purpose a very good approximation can be obtained expanding « into a series and taking the first two terms 
eT ~(1 — Ts 5 
Eq. 5 substituted into Eq. (4) gives 
fi 2fw,,s8 + w,,70 = shy + 2w,A + 0 4 V/ 56 +MsKo0 + 
VW 5K 6s0 — MsKo% — TM 5Kos0 — T\Ms5Kos*0 + T\MsK4s0 ha 
rranging 
1+ 7W3Ke6 + (2ew TMsKo — M5Ko)s + (w,2 — M5Ko)]0 = 
1+ TAM 5Ko6)s0 2 5A 60 V/ A tb 
Or 
| fa + TMsKo — M5Ko [= - ed — rw, — M5Ko) + 0 VW 
a 1+ TM3Ké \l + TMK en i + Tek TMs! . 
Tos fy the notations let 
2iw, + TM5)Ko — MsKo\|/(1 + TM5K6) = 26a 6 
wy? — M5Ko)/(1 + TM 5K6) = wn’ 7 
wher ind ¢- are equivalent natural frequency and damping ratio of the controlled system with the time lag 
Phe 
+ QWowne + wne2)AS) = S70) + [(2tw, — W5K6)0) + O)s/(1 + TM5K6) + (M56 1+ 7AMgh td 
Let us shorten the notation as follows 
2tw W5K6)0. + 6 A(1 + TM5Ko6 8 
VW 56 A(1 + TM5K6)) = ) 
2Few@ne = 10 
w =R 1] 
Then Eq. (4d) solved for @(s) gives 
As) = O}(s? + ays + a s(s? + bs + R)]t = Oo} + + + a i 3 1 
where the roots of the characteristic equation are 
§ => @a@ -- >) l 5 
The inverse Laplace transform of Eq. (4e) yields the time history of the motion (for f > 7) as 
6 = 0)}(a)/By?) + (1/88) [(a? — B? — aa + + B — 2a @! sin (Bt + W); 13 
where 
y = tan }|B(a,; — 2a a? — B? — aa + a)} — tan 3/-a 14 
Bo? = a* + B 15 
REFERENCES strength of the resulting Mach shock has been obtained correct 
Automatic and Manual ¢ papers contributed to the conference at to the second order of small quantities by Tan! and Chester.? 
Cranfield, 1951, London: Butterworths Scientific Publications, 1952 The same problem is considered in this note, but for a thick wedge 
Gardner, M. F., and Barnes, J. L., Transtents in Linear Systems, Vol. I and weak incident shock (an acoustic pulse In the latter case 
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The Shock Strength in the Regular Reflection 
of Weak Shock Waves 


William Chester 

Mathematics Department, Bristol University, England 

December 7, 1953 

Sigs A PLANE FRONTED SHOCK be incident upon an infinite 
wedge whose leading edge is yawed through an angle 8 rela- 

tive to the direction of propagation of the shock (see Fig. 1) 

If the faces of the wedge are only slightly inclined to the direction 


of propagation of the shock, Mach reflection takes place. The 


regular reflection takes place, and the boundary of the perturbed 


region consists of the Mach cone with vertex at the junction of 
the shock and the wedg« leading edge, and possibly by the tan 
gent plane to this cone from the intersection of the shock and the 
face of the wedge [this latter part of the boundary is the ordinary 
reflected shock which occurs in the elementary theory of reflec 
Is so in 


tion and does not appear when the face of the wedge 


clined that there is a region of shadow and, hence, a diffraction 
effect 


The pressure and sonic velocity in the undisturbed flow ahead 


see Figs. 2(a) and 2(b 


Corre 
the 


respectively 
flow 
shock, and in the region of perturbed flow, will be distinguished 


of the shock will be denoted by py and ap, 


sponding quantities in the region of uniform behind 
by the suffixes 1 and 2, respectively. The dependence on time 
is eliminated by taking the origin of a Cartesian system of co 
ordinates O(x’, y’, at the vertex of the Mach cone with the 
=’ axis along the leading edge of the wedge and the y’ axis per 
pendicular to the direction of propagation of the incident shock 
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0 < ao < w/2 and (b) — 2/2 < a, < 0,0 < a < 2/2 
On the assumption that (p; po) is small, the main flow velocity 


behind the shock is approximately do/sin 8 parallel to the leading 
edge, and the pressure in the region of perturbed flow satisfies 

0?p/ox'? + O7p/dy"? = cot? B(d2p/ds”? (1 
where 


p = (pp pi) /(p p (2 


Since there is no fundamental length in the data defining the 
problem, the pressure must be a function only of x and y, where 


2 


x = x’/(2’ tan B), y = y’/(s’ tan B (3 


If, further, the Busemann transformation be made, namely 


x = (2p cos 0)/(1 + p?), y = (2psin @)/(1 + p? (4 
Eq. (1) reduces to Laplace’s equation with (», 9) as polar coordi 
nates. 


Angular distances remain unchanged by the transformation 
Furthermore the Mach cone becomes the unit circle in both the 
(x, y) and the (p, @) planes. We thus require a harmonic function 
which satisfies certain boundary conditions on the unit circle and 
the lines 0 = am, 0 = 27 a2, Where a and a» are the angles 
which the faces of the wedge make with the x’ axis. Such a solu 
tion can be obtained by a method similar to that used by Light 


hill.* 


When the conformal transformations 


}p exp 12(0 a)] +” 2n - ie = Z (5 
Mtinmn=3 = (1/2)[Z + (1/Z (6 
are effected, that part of the boundary consisting of 6 = a 


(between the origin and the unit circle) in the (p, @) plane then 


becomes y; = 0, — x, < —lin the z,-plane. The unit circle 
between 0 = a, and @ = 27 a2 becomes y; = 0, i <# < I 
and finally @ = 27 — a» becomes y; = 0, 1 < x - 


In the 2-plane, the boundary conditions on the real axis may 


be written 
Op/Oy = 0, “i > 1; Op/Ox, = O w<i @ 


together with 


NCES 


* 
rar ye eee sin ast 
dx = : 
€ . OX SIN QeT 
where 
T T 2a a a 
We now define the complex function 
i = OP/Ov + U(Op/Oxr 9 


The boundary conditions, expressed in terms of I], are that o1 


the real axis Wis real for x 1 and wholly imaginary for 


Near = cos ait, Eq. (8) suggests that 
1 sin air l 
W~ 
T SIll a) T T COS @,T 
with an analoguos relation near = COS a»T 


A function which satisfies all these conditions, and is integrabk 


it all points in the upper half-plane except the specified singu- 


larities, is 


W = (1/7 


i(sin ayT (1 + COS a)T 


(1/7)i(sin aor l COS a2T 10 


It is, of course, possible to deduce the actual pressure field using 
Eq. (10) and the pressure distribution on the Mach cone Phis 


has, however, already been investigated by Keller and Blank 
q § . 


using a different technique, and Eq. (10) is sufficient for the pur 
pose of this note. 
The above theory predicts a singularity in the normal deriva 


tive of the pressure at the Mach cone. For if 0f/0mn denotes the 


pressure derivative along the inward normal to the circle p = 1 
in the (p, @) plane, and if (7, @) are polar coordinates in the 
(x, ¥) plane, we have near r = 1, 

p (P)rar™ (1 p)(Op/on) ~ [2(1 y (Op/on 1] 


and hence 

Op/or ~ (i= (Op/on 12 

Now Lighthill®’ has shown quite generally that the true behav- 
ior near the Mach cone depends on the sign of 0p/on. If 

p (po)p ~ pi V1 y Fie 13 

then the Mach cone is a linear approximation to a shock wherever 

F(@) > 0. The discontinuity in pressure across this shock is 


Ap» = (3/4)v(y + 1) pr Vi8/lai( Vi? — a2)}} F200) (14 


where y denotes the adiabatic index 


When F(@ 
followed immediately by a relatively stronger 


shock 


se cond order 


< 0 there is probably a weaker (third order 
expansion 
In the present case, it follows from Eq. (11) and (13) that 


. r ot 
F(0) = {12"*(p, — po] /oVi2}(Op/on 15 


f 


are, to the present order of approximation, a» and 
ind (15) show that, when Op/or 


Since a; and V 
ao/sin 8, respectively, Eqs. (14 


0) 
Ap. = 3(y + 1)/27 cos? 8 KA — P )2 Po|(Op/on . 16 
Now 
O; dz, Op 
= li 
on dp |dy, 
and, on the circle p= l, 
dz; E 
= rsin7(é — a 18 
dp 
x = —cos 7(6 — a 19 


by virtue of Eqs. (5) and (6). It follows, with the help of Eq 
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Fic. 3. Diffraction pattern produced in a corner for a; = a: = 
a and (a) r/2 <a < 27/3, (b) 27/3 < a < 34/4 
(10), that on p = 1 
of I j sinfa,;r/(27 — a a 
On 2r — ay a Icos o cos|a,;7/(2r a a 
sin[asa/(27 a a ( 
2) 
. 4 . ») 
COS @ + COS\|aer/(2r a a 
where @ = (0 — a )ar/(27 a a), so that 0 < @ S zwonthe 
Mach cone 
In the symmetrical case, for example, with a; = a. = a, Eq 
20 simplifies to 
sin |ra/(7 a 
Op/on = (21 
(3 a); cos 26 cCos| ta/(7r a 
When 0 a 7/2, Op/On is positive for ra/2(r a @ 
T ra/2( 4 a), that is for 2a 7] 2( + a) if terms of the 


third order are neglected. Thus if terms of the third order are 
neglected a shock is predicted only along that part of the Mach 
cone which contributes to the boundary of the disturbed region 
(the continuation of the plane reflected shocks—see Fig. 2 
The remaining part of the Mach cone probably represents a 


weaker shock followed immediately by an expansion 


When 7/2 < a < 27/3, Op/On is positive in the intervals a 
6< Ar a) and 2a 6 (27 — a These sections of the 


Mach cone still represent the continuations of the plane reflected 
shocks [see Fig. 3(a As a is further increased multiple re 
flections of the incident shock wave occur at the two walls before 
the actual corner is reached,‘ and the flow pattern in the neighbor 
hood of the corner is determined by the last pair of reflected 
for example, shows the configuration near the 
corner for 27/3 a 32/4. The factor (7 — a 
suggests that strong shocks will be produced by 


waves. Fig. 3(b 
in the denomi 
nator of Eq. (21 
reflection in corners of small angle 


Diffraction by a half-plane is also a simple special case. With 
—a = a = a, Eq. (2()) becomes 
Op/On = [2 sin a/2 cos (0/2 + a/2 
lx cos (6 + a COS a } (22 


In this case the reflected shock appears on one side only of the 


plane (see Fig. 2(b)), and is continued along the Mach cone as a 
weaker shock whose strength is deducible from Eqs. (16) and 
(22) as far as the point 6 = (7 — a Finally, in the region of 


shadow the Mach cone represents a continuation of the incident 


shock as far as the boundary 
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Allowable Stresses for Channels and Zees in 


Bending 


B. E. Gatewood and D. W. Breuer 

Mechanics Department, USAF Institute of Technology, Wright- 
Patterson Air Force Base, Ohio 

January 4, 1954 


5 Re NOTE CORRECTS AN ERROR in the paper by Gatewood and 
Williams! ‘‘Allowable Compressive Stresses in Aircraft 
Structures’’ and presents some test results on channels in bending 

Several readers have been so kind as to point out an error in 
the bending buckling coefficient (A) curves in Fig. 6 of the orig 
inal paper, the error being of the order of 15 per cent conserva- 
1 shows the corrected 
1.25 and 1.50 

V in Fig. 1 (Fig 


inalogy to channels in com- 


tive in the region of },/b, = 1 to 5. Fig 
curves, as well as additional curves for ¢,,/t, = 
The table for the maximum stress coefficient 
6 of original paper) was based on 
pression and on seven bending tests on 24S-T formed channels 
Recently, the authors have conducted a few more tests on stand 
(zero shear 


ard 24S-T channel extrusions. A constant moment 


was applied to two channels back to back, separated by spacers 


outside the test region. To obtain failure in the test region i 


stead of at the ends, the thickness of the compression leg was r¢ 
Fig. 2 shows the 


duced about 10 per cent in the test region 









































test results as compared to the calculated values K and M 
M | 
b, |M-formed | M-extrusions; 
{in "T"cond. | formed in "O"cond 
0-25 0.85 0.80 | 
40 0.75 067 
75 0 67 0.60 
=Compression Leg | 
120+ | 
1a ! | 
tw } 
*=1.50 
100+ ‘ 5 
<Moy be different in"b" | 
P “and “t" trom ¢ompression 
0.80tT 
| 
K | 
0.60 + t 
0.407 
0.20 
4 + n 2 
1@) 2 a 6 8 10 
Dw 
oe | 
Fic. 1 Vand A values for channels and zees in bending 
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Fic. 2. Test results on bending of 24S-T extruded channels 


of Fig. 1, together with the crippling stress curves of Fig. 4 in 
reference 1, were used to calculate the allowable crippling stresses. 
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Propagation of Turbulence from an 
Instantaneous Point Disturbance 


Morton Mitchner 

Research Fellow, Combustion Aerodynamics Laboratory 
Harvard University, Cambridge, Mass 

January 6, 1954 


Sagem has recently proposed a mechanism for the manner 
in which transition occurs from a laminar to a turbulent 
boundary layer. The mechanism assumes that unavoidable 
sources of disturbance carried within the laminar boundary layer 
may appear at random as “‘point’’ sources of turbulent spots. 
Once initiated, these individual turbulent regions will continue 
to grow in extent as they travel downstream until they coalesce 
with adjacent turbulent regions and eventually produce a fully 
developed turbulent boundary layer. The following contains a 
brief account of the flow distortion resulting from an artificial 
‘instantaneous, point’’ disturbance produced in two different 
types of shear flow. 

If the steady laminar motion of water flowing in an inclined 
rectangular channel is externally disturbed by allowing a drop 
of water to fall into the fluid, the development of the resulting 
region of turbulence may be observed visually by the roughness 
of the free surface associated with the turbulence. Observations 
were made at a constant mass flow rate of 0.81 Ib./sec. ft. (or at 
Reynolds Number of 1230 corresponding to a water temperature 
of 21.0°C.) for channel slopes of 1.18°, 2.98°, and 5.26°. The 
corresponding mean velocities were 17.5, 25.0, and 31.0 in./sec., 


respectively. In all cases the developed turbulent spots were 
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observed to be very closely geometrically similar in shape Each 
spot is characterized by a well-defined upstream boundary with 
a more vaguely defined downstream edge, a typical shape being 
indicated in Fig. | The half-angle of the spots with respect to 
their point of origin, a, was found to remain essentially the same 
The half-angle of the spots with respect to 
is 50°, 48°, and 45° 


and equal to 6.6 
their own vertex, 8, was found to vary 


respectively The velocity of the upstream edge of the spots, 


s’, varied as 14.4, 20.2, and 25.4 in./sec., respectively, while th 
velocity of the downstream edge, s”, varied as 15.9, 23.0, and 
28.8 in./sec., respectively. The experimental values for were 


found to agree moderately well with a theory treating the up 
stream edge as a particular form of hydraulic jump.? 

A low turbulence level wind tunnel was emploved to produce 
critical laminar boundary layer over the length of a 30-in. flat 
plate at a free stream velocity of 550 in./sec. Point disturbances 
having a time duration of less than sec. were produced peri 
odically 10 in. from the leading edge, and the development of the 
consequent turbulent flow patterns downstream was observed 
with a hot-wire. Disturbances were produced both by popping 
a small projection into the laminar boundary layer from the 
undersurface, and by producing a spark discharge in the bound- 
ary layer. The plan and side view shapes of an average spot are 
indicated in Fig. 2. At a distance of 14.5 in. from the disturb- 
ance source, the upstream and downstream edges of the turbulent 
spot were found to have velocities of 310 and 450 in./sec., 
respectively. The half-angle of the spot with respect to the 
disturbance origin, a, was found to be 8.6 It is of some 
interest to note that the vertex of the turbulent spot for the flat 
plate case points downstream, in exact contradistinction to the 
free surface case. The inference is that in the latter case the free 
surface, in addition to the turbulence, plays an essential role in 
establishing the spot shape 

The observation that turbulent spots may be created artifi- 
cially does not offer a complete verification of Emmons’ mecha- 


nism of transition by any means, and considerable work 
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remains to be done. The observation that spots appear to have 

heads projecting well above the flat plate surface suggests the 

interesting that with 

the ragged turbulent high spots in the intermittent region above 
laver. If this 

case there should exist a high degree of correlation between the 

hot-wires 


speculation these heads are associated 


a fully-developed turbulent boundary were the 


properly time-delayed turbulence readings of two 
placed in the intermittent region and separated along the direc 
tion of flow It is also suggested that an obstacle, such as a 
cylinder, placed in a uniform flow acts as an intermittent source 
of turbulent spots, and that the wake consists of a moving pat 
tern of such spots, the observed wake puffs being merely the 
passage of successive spot wings 

The author would like to express his thanks to Dr. A. A. Town 
send for both his invaluable advice and his permission to use the 


low turbulence level wind tunnel at Cambridge University 
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An Integral Equation Occurring in the Theory 
of a Slender Quasi-Axisymmetrical Body 


Byrd 

Aeronautical Research Scientist, Ames Aeronautical Laboratory, 
NACA, Moffett Field, Calif 

January 29, 1954 


ZZ, 


W T. LORD! in a recent report derives an integral equation 
* for a quasi-axisymmetrical, slender body having a spec- 


He states, however, that no pub- 


Paul F 


ified sidewash on the wing. 


lished solution has yet been given. It is the purpose of this note 


to present an inversion of the equation. 
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The integral equation is 


interval of integration and all variables 


/ 


where y lies outside the 


are real. This equation may also be written in the form 

1 x, a =e la 
where the prescribed function 7 is related to the sidewash, and 
H is 


the be dy 


an unknown function defining the boundary condition on 


(See reference 1 Eq. (la) can easily be inverted by 


methods of complex variables or potenti il theory, and lends itself 
also to direct treatment by expansion of the kernel We shall 
employ here the latter method 

kernel in I'chebycheff’s 


Formal expansion of the terms of 


polynomials? of the second kind yields 


l l ~ 
= > ] é » 
. = aa 4 V1—-—¢j;=0 
where the polynomials are given by 
l Ne Ud+NvVi—{ ST +ivi — gs > 
e ) 
cos @ sin 6 dé 4 
Assuming now that the unknown 
H(x, £) = Yo bmailx) Ungil (4 
0) 
we thus have for Eq. (la 
vi ] ‘ : y 2 
T(x, a) = p ® } ee 4 
m=0 0 Z ViIi- 
If we further assume that 
T(x = S U \ a) 
fo] 
it follows from Eqs. (5) and (6) that x/ = m4 Hence, from 


Eqs. (3) and (4) we have the expression 


H(x, ¢) = Y~£ YS (1+ m) engl / 4 
m= 0 


\ 1— ¢*? cos @ in 6 dé r 
The series occurring in this equation for /7 is summabk 
i | 
From Eq. (6 
—_ d 
S 1 - & d \ ao = 7 \ ao + ao 7 i 7 
— lea 
, 0 
and thus upon setting o = ¢+ 1V 1 — ¢* cos @, one can write for 
Eq. (7) the inversion formula in the closed form 
j é Vi a 
H oe & = 7 xv, @ T 7 7 rE @ lo ‘“ 
wT J eeivi lo 
or finally 
9 
H(x,¢) = -Im[(E + 72Y1—-2)7(yF+1V 1 -— & Sa 
T 


where the symbol Im means that the imaginary part is to be taken. 
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On Heat Transfer in a Compressible Laminar 
Boundary Layer 


Alfred H. Silver 
General Electric Company, Cincinnati, Ohio 


December 24, 1953 


NOTATION 


( specific heat at constant pressure, ft.lbs./slug °F 
h = film coefficient of heat transfer, ft.lbs./sec. ft.2 °F 
k = thermal conductivity, ft.lbs./sec. ft. °F 

p = pressure, lbs. per sq.ft 

f = function of x defined in text 

u = velocity in x direction, ft. per sec 


= velocity in y direction, ft. per sec 
= value at wall, subscript 
distance along airfoil measured from forward stagnation 


point, ft 


distance perpendicular to and measured from airfoil 
surface, ft. 
I = modified Bessel function of first kind 
Q = unit heat flux, ft.lbs./sec. ft.’ 
[ = total temperature, 7; + (u?/2c,), °R 
Du Ou 
Vomentum: pu + pi = 
on oy 
- 0 
Continuity: (pu) + 
; ed 
; oT, o7 dp 
Ene rgy: pc pit T pit = 0 
. Ox oy dx 


TICAL SCIERNCES MAY, 1954 
7, = static temperature, °R 
I = gamma function 
6 = function of x defined in text 
# = momentum thickness of boundary layer, ft 
u = absolute viscosity, Ibs. sec. /ft 
v = kinematic viscosity, ft.* per sec. = u/p 
p = mass density, slugs per cu.ft 
a Prandtl Number ywc,/k, dimensionless 
T = shear stress at wall 
y stream function, pu = Oy /Oy, pv = —Oy/Ox 


ANALYSIS 


o oe SOLUTION OF HEAT TRANSFER in a two-dimensional lami 
nar gas boundary layer, considering compressibility effects 
ind with arbitrary pressure and surface temperature distribu- 
tions, involves the simultaneous solution of the momentum, con- 
tinuity, and energy equations. The exact solution of the equa- 
tions has not, to the knowledge of the writer, been found as vet 
This note, which is a condensation of work done for General 
Electric in 1951, attempts to extend the approximations of Light- 
hill! and Tifford? to include the effects of compressibility 

For the analysis, c, is considered constant, o is considered 
constant and near unity, and the effects of curvature of the sur 
face and pressure gradient normal to the surface are neglected 


rhe equations to be solved are 


Cy oO ( oT, ) 
be 
ao OV oy 


Transformation of the equations into the x, Y coordinate system yields 


Ou dp 
pil == t 
Ox dx 
oO 7 dp plc re) 
puc, =u- t 
Ox ax o Oy 
Multiplication of Eq. (la) by w and adding yields 
or i124 m) u> || 
= pul T+(eo-— ] (4) 
ox a Oy | oy 2 1 


The boundary conditions on the total temperature are: 
T(O, ¥) = T(x = () T(x, 0) = To(x 


Since the analysis is concerned with heat transfer, which is a 
surface phenomenon, and since o is assumed near unity, the 
second term in the bracket is neglected in comparison with the 


oT 1 oO oT : 
= puu (4a) 
Ox a OY oy 
A final assumption about the nature of the gas properties is 
The applica 


first. This yields 


that the viscosity can be expressed as ph = C)/ 
bility of this assumption has been shown by Chapman and 


Rubesin This means that py = Cop and 
O(pu)/OY = C2(0p/dy) = O 
Therefore 


OT/Ox = (pp/o) (0/dW) [u(OT/oy) (4b) 


pu 


0 Ou 12) 
pul 1a) 
oY Ox 


O/ : Ou \* 9 
peu + pp-u" va 
oY OY 


Similarly, from Eq. (la 


Ou —ldp re) Ou 
= + pu u Ib 
Ov pu dx oy j 


The simultaneous solution of Eqs. (1b) and (4b) is a formidable 
task in the general case. Therefore, a power series approxima- 
tion to the velocity distribution which is valid close to the wall 


will be used. 


Let 
(/V2 = ala + vy Tt CX) + 
Since u(x, 0) = O 
r(x, O To(X = pw pulou/oy)|y=0 
a(x Oand b7(x) = ro(x)/pum 


Applying Eq. (1b) at Y = O yields 
dp/dx = pwr ,u(o/oy 
(dp/dx) = (3//2)b*% 


(1/pw7n 


This gives 
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T(x) 27/2 (” ) 1 dp. 
= y + y 4 
2 Pwht 3 T Pwhtw ax 


If the first term only is used, the velocity distribution of Light- 
This, 
in the direction of flow 


| is obtained. however, neglects the effect of pressure 


on the velocity distribution 





ested that the term involving dp/dx be accounted 
This 


Tifford sugg 
yme wav while retaining only the first power of y¥ 


one, in this case, by letting 


+ 6 dp dx y 


[he validity of the use of 6 and its exact nature will be deter 


ifter the analysis is complete by comparing the results of 


minea 


can be solved by Laplace Transformation obtaining a Be 


Ea. (¢ 


its power series form, the following is obtained 


T(t) = Tet) 21 - ae [ ss 
( (4/3) L Jo (f - 


0 due to the change in 


Since T has a finite discontinuity at ¢ = 
surface temperature, the coefficient of y ‘* in Eq. (7) can be expressed 
9/2 
2/3 
I'(4/3) Jo 
rherefore 
] 2/3 i 
T(x) = T(x l— — 
{ r(4/3) Jo 
where /; is the value of fat which x = & But 
1 x 
0 ¢ 
so that 
9/9 "% “a 
T(x) = To(x 1 — ——— \ 
t r(4/3) Jo , 
The heat transfer from the surface, Q(x), may be determined from 
o7 
Q(x) = —k — 
oy 0 
or 
k V puwk lo 
Ax = 
(4/3 Qu 4 J/0 
Since o {, the effective and total temperatures are nearly equal. 


Q(x k (‘*) V pwF 
mee Te | (UNSDT AS Qu 


For the average heat-transfer coefficient of a fluid of constant 
properties on an isothermal surface, Eq. (10) simplifies to 


*) 


Vv F(z) dz 


lifford compared Eq. (22) with the solution obtained by nu- 


il integration of the differential equations for isothermal 
He determined that 


meric 
wedges with widely varving wedge angles 
an excellent correlation (+1.8 per cent) could be obtained if the 


value of 6 be chosen as 
3)Ox)/a (4 


6 = (4 


ind the resulting value of F be multiplied by 0.980°- 


FORUM 353 


the equation obtained with those from the special cases 


have been solved exactly 


Letting F(x t) + 6(dp/dx) and remembering that py is 
independent of ¥, Eq th) yields 
o D7 2) 17, Of 5 
V 2puF Ox dy (¥ a - 
Letting 
. wang 
J0 
O7'/Ot O/oOy) ly OT /Oy¥ 6 
subject to 
T(0, ¥ T(t, 0, TUt,0) =7 


solution If the Bessel function is written in 


ssel function (/ 


temperature from the free-stream value to the stagnation point 


as the Stieltjes integral 


“é dT\(t 
dT ; { 
Lee ( 
Vv 2 ppk d 
2puF(z) dz dT\(é |. X)y 7 { S 
O7 ; O7 
= — k,.pu 
oy 0 OY y=0 
~» 4 
\ 2ouF dz dT): 4 
Therefore 
"* a a 
V 2puF(z) d d7 10 
0 Jt 


It is therefore recommended that Eqs 


, 1/3 / 
F(x) = 0.98e0°-°% T + 
gs’ ax 
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Buckling of Thin Cylindrical Shells Subject to In 1934 Donnell* presented a second-approximation shel] 
Hydrostatic Pressure theory for finite deflections in which two of the displacement 
components were eliminated and the following relation between 
the stress function F of the membrane stresses and the radia] 

William A. Nash displacement w was developed: 
Associate Professor, Department of Engineering Mechanics, te A | —_—* ( my dew 24 Ed% 


University of Florida, Gainesville, Fla ey 
Ox? Os? Oxo oc? s? R On 
January 19, 1954 7 ‘ ° ° ™ 
Here x and s are axial and tangential coordinates in the middk 
surface of the shell, R denotes the mean radius of the shel] 
— ANALYSES have appeared in the literature for the and K involves the ratio between the initial deviation from per 


elastic buckling of a thin cylindrical shell subject to hydro fect shape w) and the displacement during loading w. Th 
static pressure.'~® These have all been based upon the classical — second-order terms on the right-hand side correspond to mem- 
small deformation theory of thin shells presented by Love,’ and brane stresses due to finite displacements, while the last term 
have assumed a geometrically perfect cylinder prior to loading corresponds to membrane stresses due to the curvature. For 
The procedure adopted in most of these papers has been that of — the case of small displacements of an initially perfect shell, this 
minimization of the total potential, which leads to a stability equation reduces to 
determinant for calculating the buckling pressure. The labor nee R Rs Sick : : 
of computing the elements of such a determinant is considerable Kor /ox") + (0"/0s°)/ °F = —(E/R): (O*w/ax* é 
and it is the purpose of this note to present approximate formulas Thus, if w is assumed, then F can be readily obtained. Knowing 
for the elastic buckling pressure for several boundary conditions wand F the internal strain energy of the shell is given by* 


a Eh . _— j (ow O*w \? o*w \? ow Ow |t 
E= + 4- O(] — vp) — . > dx ds + 
241 — »?) Jo So (\ox? © Os? oxds Ox? Os? If 
h . 2rR (/a2F Ot F\? OF\? wF oF) 
; + + 2(1 + v) _ . dxds (3 
2E Jy . (\ox? © ds? Oxds Ox? Os? f 


where / is the shell thickness and vy represents Poisson's ratio. 
The potential energy of the forces acting upon the shell has previously been demonstrated to be 


a! 727R i — so 3 
W, = pR Sa St os ae as { 
a J | | 2 ‘OR? 4 Jews 


where L is the length of the shell. Consequently, the total potential is given by the sum of /, W), and a constant term which 
corresponds to the mean axial and circumferential stresses 
Let us first consider the case of a clamped-end cylinder subject to hydrostatic pressure. The shell is not reinforced. If the 


deflection of a generator is represented by the function 


w = ah cos (ms/R)(1 — cos 26x) (5) 
where a is an arbitrary parameter, m is the number of waves in the circumferential direction, and 6 = 2/L, then from Eq. (2 
we have 
; $E62ah ms 
F= —- — i cos 26x (6 
[462 + (m?/R?)|?R R 


Minimizing the total potential with respect to the parameter a and solving the resulting equation for the buckling pressure p we 


obtain 


Eh’ j 3m4 45°m?| $Eh6* 
4 86! + - + 
24(1 — v?) { 2R3 R? {© [462 + (m?/R?®)]?R? . 
0, = ‘ 
f 3m? R8? 34 
(4k 2 1R§ 


Because of the simplified expressions used for the bending strain energy, the above equation for buckling pressure is accurate only if 
the square of the number of circumferential waves m is large compared to unity. 

Secondly, the case of a cylindrical shell having simply supported ends and subject to hydrostatic pressure may be treated by 
this technique. This is the problem considered by von Mises! and if the deflection of the generator is taken to be 


w = ah sin (ms/R) sin 6x 8 
as he assumed, then after minimization of the total potential, the buckling pressure is found to be 


m?|? Ehé* 
p ot 
R24 $[5? + (m?/R2)|2R? 
m? R86? ] 
jm? Re 1) 
l44R° 8 4Rf 


$8(1 — v?) 
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TABLE 1 
ee HN, R, d B B, D, I 9 } 
Mod In In In Ref. 1 Eq. 9 Eq. 7 b h Ref. 4 Ref. 6 
1 32 0.03 Ss 7.0(6 7.2(6 10.4(7 0.12 0.40 7 9(6 GQ R16 
2 32 0.050 8 23.1(5 24.1(5 35.2 (6 0.12 0.40 29.7 (6 32.5 (¢ 
16 0.053 Ss 51.2(7 53.2 (7 74.5(8 0.12 0.40 65 5(8 68 4(8 
{ 32 0.109 8 152.2 (4 162.1 (4 235.2(5 0.12 0.40 197.7 (5 213 .5(5 
5 Ss 0.0635 Ss 164.7(9 166 .6(9 246.5(10 0.12 ) 40 296 0(10 210.91 
{ 8 0.0901 8 384.9(8 390.6(8 538.7 (9 0.12 0.40 163.7(9 $92 5(9 
7 5.08 0.065 13.38 37.6(15 38.1(15 189.9(17 0.1875 0.719 175.1(17 179 O(16 
S 6.686 0.1324 13.41 6138.2(11 617.5(11 ) 852.8(12) 0.23 1.0824 781.0(10) 823.2(10 
This is equivalent to the result von Mises obtained using three mind that all of these pressures are for a geometrically perfect 


differential equations of equilibrium. It was also presented by 
Salerno and Levine? as a special case of their analysis of a cylinder 
reinforced by elastic circumferential rings 

[It is of interest to compare the numerical results indicated by 
Eqs. (7) and (9) with those found from existing theories derived 
from minimum energy considerations in which three orthogonal 
components of displacement were considered and the potential 
minimized with respect to at least three parameters. In Table 1 
ire listed the geometries of eight cylindrical shells together with 
A is 


for simply- 


buckling pressures given by various theories. In column 
listed the buckling pressure indicated by von Mises! 
supported cylinders, and in column B appears the value given by 
Eq. (9 
clamped ends* and column D lists the buckling pressure for this 
Eq. (7). 
buckling pressures indicated by the theories of references 4 and 


Column C lists the buckling pressure for a sheli having 


same shell as given by In columns E and F appear 
6, respectively, for multiple-bay cylindrical shells, each bay 
having the geometry indicated in Table 1 but reinforced by 
rectangular circumferential stiffening rings having the width 5 
and overall ring-shell depth iy) indicated. The pressures given in 
C and D could, of course, be regarded as buckling pressures for 
i multiple-bay shell having infinitely rigid reinforcing rings, each 
bay being of length Z. The number in parentheses after each 
buckling pressure is the number of circumferential waves which 


corresponds to the minimum buckling load. It is to be borne in 


An Empirical Method for Calculation of the 
Growth of a Turbulent Boundary Layer* 


Donald Rosst and J. M. Robertsont 
Ordnance Research Laboratory, The Pennsylvania State College 
November 4, 1953 


SUMMARY 
A new, simplified method for calculating the momentum thickness of a 


mmetrical turbulent boundary layer has been developed empirically 
In regions of strong adverse pressure gradient, a simple algebraic expression 
may be used involving an exponent which is a function of the Reynolds 
Number Limitations of the method are presented 


INTRODUCTION 


; pe USUAL EQUATION for the growth of a two-dimensional 
boundary layer is the von Karman integral momentum equa- 
tion for the rate of change of the momentum thickness @ as a 
function of the wall-shear-stress coefficient c; and the pressure 
gradient 


* Work partially supported by the Office of Naval Research under Project 
NR 042-139 
T Associate Professor of Engineering Research. 


+ Professor of Engineering Research 


shell having no initial imperfections and being free from initial 


stress. The presence of initial imperfections is known to reduce 
the buckling pressures from those given by classical treatments 
of idealized shells 
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do/dx = c;/2+ }(2 + H)[(0/pu;? (dp/dx)}; (1 


where // is the shape parameter and u, the velocity outside of the 
boundary layer. This equation was originally derived by von 
Karm*n! in 1921 for laminar boundary layers and later applied 
to turbulent boundary layers by Buri? and Gruschwitz Al 
though the equation has generally been verified for turbulent 
flows with zero pressure gradient, and by Clauser‘ for equilib 
rium flows of all types, it has been found faulty for nonequilib 
rium turbulent boundary layers both with and without adverse 
pressure gradients. The cause of the discrepancies was dis 
covered in 1949 by Newman and Wallis,® who pointed out that 
the turbulent normal stresses are not negligible and must be in 
cluded in the differential equation for the growth of the mo 


mentum thickness 


The revised equation derived by Goldschmied,’ Bidwell,” 
Van Le,’ and Ross!’ may be written 
fy) 
dé Cc 6 dp ] 3) 
= +(2+H - + u’? — v’*)dy 2 
dx 2 pu,;? dx u;? Ox 


- 


where p refers to the wall static pressure. The term involving 


the anisotropy of the turbulence has been related to the shear 
stress and thence to mean flow quantities by Ross,'®!! resulting 


in a new approximate differential equation; while Rubert and 


Persh'? have obtained an empirical correction term of similar 
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form. In any case, it is now agreed that some correction of the 
original von Karman equation is required for turbulent flows. 
We may write the revised equation 


do/dx = c;/2 — \(2 + H)[(0/u)(du;/dx)|{ + [Co(dH/dx)] (3 


where the term involving the wall pressure is replaced by one 
with the outer velocity, and the correction term includes the 
fact that the static pressure at the outer edge of the boundary 
layer is not equal to that at the wall. 

Eq. 3 reveals why Clauser* found good agreement with Eq. 1 
for equilibrium boundary layers even in fairly strong adverse 
For such flows, // is practically constant and 


pressure gradients 
The equation also shows why 


the correction term is negligible. 
some cases of redistributing flows have not agreed with Eq. 1, 
even though there was practically zero pressure gradient. For 
such flows, the correction is negative and the momentum thick 
ness grows less rapidly than as predicted by the simple von 
Karman equation. 
Integrated forms of Eq. 1 have been developed by several 
authors, often based on erroneous assumptions concerning the 
None of the 


equations is particularly simple to use and all of them are in- 


behavior of the wall-shear-stress coefficient. 
accurate near separation. The solution presented in the present 
paper is quite simple mathematically, and it is usually more 


accurate than the other methods 


EQUILIBRIUM BOUNDARY LAYERS 


For conditions close to equilibrium, // is practically constant 
and the only major problem is the determination of a suitable 
equation for the wall-shear-stress coefficient. The recently de 
veloped equation of Ludwieg and Tillmann! 


Cy = 0.246/( 10.6754 Rg?-25) (4) 


is applicable at any value of the shape parameter. Expressing 


Eq. 3 as a difference equation: 


AO = (c;/2)Ax — (2 + HW)6Alnu, (5) 


and using this with Eq. 4, step-by-step integration is possible 
This is the usual method for the solution of the von Kérman 
equation. However, because of the relative constancy of c¢ 
and H and the usual small pressure gradient, large intervals may 


be used and the solution requires only a few steps 


STRONG ADVERSE PRESSURE GRADIENTS 


For comparatively large pressure gradients, such as those on 
the suction surface of a lifting airfoil or in a diffuser, the step-by- 
step solution requires many more steps, as well as simultaneous 
solution of an equation for //. Here a new solution is presented. 
The complete momentum equation (Eq. 3) may be written: 


d0/0 = — duy/u, i\2 + H — [c;/2(0/u,) (du,/dx)| — 
Cu(dH/du,)| (6 


All of the functions in the brackets are variable with distance; 
but it has been observed in a large number of cases that their 
sum is approximately constant for any single experiment. If 
we make this assumption and replace the quantity in brackets by 


2 + G, then the solution of Eq. 6 is simply: 


0/0; = [u;/m]?2 7° (7) 


where the subscript 7 refers to some initial station. This solution 
of the integral momentum equation is very simple and allows the 
replacement of the usual laborious step-by-step solution by an 
algebraic equation. It remains to confirm this solution experi- 
mentally, determine its limitations, and evaluate the exponent 


CONFIRMATION FOR TWoO-DIMENSIONAL FLOWS 


To verify a solution of the momentum equation of the form of 
Eq. 7, the data from more than twenty two-dimensional experi- 


> 


ments were plotted in the form of @ as a function of “,? on log-log 
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paper. These experiments were chosen to cover a wide variety 
of two-dimensional configurations. The result of the »g-log 


plot was encouraging, as a general linear trend was found for 
most of the data. Closer examination revealed consistent dis 
crepancies for weak pressure gradients, and also close to separa 
tion. By trial and error it was found that one criterion for th¢ 
applicability of Eq. 7 is that the pressure gradient term dominate 
This is expressed approximately by the inequality 


— [(0/u,) (du,/dx > Cy/4 & 


Restricting the data to those satisfying this inequality, a second 
plot was made to confirm Eq. 7. The result is presented i 

1, the key to the experimental points being given by Tab 
The general linear trend is found for most of the experiments, 
with an average value of G of about 2.8 fitting most of the data 


1 


— OG 


Particularly impressive is the verification for the diffuser tested 
by Hall' for which particular care was exercised in obtaining 
two-dimensional flow. The only divergent set of data is that of 
Wieghardt and Tillmann”® in a channel 
criticized elsewhere, and Tillmann has attributed the divergence 


These data have beer 
to observed secondary flows. Ignoring this one test, the rest 
of the data give a reasonable confirmation of Eq. 7 

The explanation of the success of Eq. 7 lies in the approximate 
constancy of the relation 


G = H — [c;/2(0/u;) (du,/dx)| — Cu,(dH/du 9 


all terms of which are usually positive in an adverse pressure 
gradient. Apparently, as the boundary layer progresses, the 
decrease in the skin friction term closely equals the increase in 
both #7 and the correction term, thus keeping G approximately 
constant. Near separation, the sum of #7 and the correction 
term is large, and the momentum thickness may be expected to 
exceed the value given by the simple equation. However, for 
redistributing flows, such as a reattached boundary layer follow 
ing separation, the correction term is negative and so Eq. 7 is 


valid even at high values of H close to separation 


AXISYMMETRIC BOUNDARY LAYERS 


For flows with rotational symmetry the momentum equation 
is quite similar to the two-dimensional one, and so a similar 
solution is to be expected. Unfortunately, bodies of revolution 
exhibit only relatively weak pressure gradients and such data are 
not applicable. However, strong adverse pressure gradients are 
found in conical diffusers. For flow in conical diffusers, the 
momentum integral equation becomes 
2 _ dH 

— Cu 10 


d(@r) du, , / 
(0/u,) (du, /dx du 


= i 
Or uy 
where r is the local pipe radius and the three-dimensional momen 


tum thickness is defined by 


The corresponding solution for the momentum thickness in 


strong adverse pressure gradients is 
(Or/O:r;) = [(ayy/u 2 12 


This is the form in which the present solution of the momentum 
equation was originally derived by the authors’ for application 
to a diffuser experiment 

To confirm Eq. 12, data from several conical diffuser experi 
ments have been used. Fig. 2 shows the results of a study of 
three diffuser angles and several entrance boundary layer thick- 
nesses carried out by the authors" preliminary to the design of a 
» 


large water tunnel. These data are seen to confirm Eq. 12 


VARIATION WITH REYNOLDS NUMBER 


Examination of Eq. 9 for G would lead one to expect a Reynolds 
number effect. At lower Reynolds numbers, both // and ¢ 
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riety 
| Pn TABLE | SYMBOLS FOR’ EXPERIMENTS 
dis 
Pini SYMBOLS INVESTIGATORS TEST OBJECT EXPERIMENT 
sie IDENTIFICATION 
; ©  FAGE, FALKNER (1930) JOUKOWSKI AIRFOIL, © R.= 125 Xx 10° 
= Q@ FAGE, FALKNER (1930) JOUKOWSKI AIRFOIL, O Re= 1.67 X 10° 
in Fig. © PRESTON, SWEETING (1943) JOUKOWSKI AIRFOIL, ©, R.= 4.2 xX 10° 
tble 1 9 PRESTON, SWEETING (1943) JOUKOWSKI AIRFOIL, 6 R,=4.2 X 10° 
ceRES, O SCHUBAUER, KLEBANOFF (1950) SIMULATED AIRFOIL, R.= 15.3 X 10° 
pee: & GRUSCHWITZ (193!) GOTTINGEN 387 AIRFOIL, 12 R,= 8.4 X 10° 
ee < _ vonDOENHOFF, TETERVIN (1942)  66,2-216 AIRFOIL, 10 R.= 2.6 X 10 
ray V von DOENHOFF, TETERVIN (1943) 65(216)-222 AIRFOIL, 8.1° Re= 0.9 X 10° 
e been A von DOENHOFF, TETERVIN (1943) 65(216)-222 AIRFOIL, 8.1° Re= 2.7 X 10° 
rgence > vONDOENHOFF, TETERVIN (1943) 65(216)-222 AIRFOIL, 10.1° Re=2.6 &X 10° 
warn oO KEHL (1943) CHANNEL AK IIT 
— B KEHL (1943) CHANNEL K3 
N KEHL (1943) CHANNEL K7Ta 
. a KEHL (1943) CHANNEL K7b 
| . 0 WIEGHARDT, TILLMANN (1944) CHANNEL IZ 208 FT. PER SEC 
en OQ HALL (1946) CHANNEL B 10.2° DIFFUSER 
vile ©  BURSNALL, LOFTIN (1951!) 66,-O18 AIRFOIL, O° Re= 2.4 X 10° 
nately A NEWMAN (1949) AIRFOIL WITH SPOILER, 0° Re= 3.3 X 100 
ection Vv NEWMAN (1949) AIRFOIL WITH SPOILER, 10° Re= 33 X 10 
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1c Fic. 1. Verification of relation for momentum thickness for two-dimensional flows in strong adverse pressure gradients symbols are 
defined in Table 1 
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Fic. 2. Verification of relation for momentum thickness in 


conical diffusers 


tend to be higher, so that G should be somewhat higher.  Ex- 
amination of the data in Fig. 1 shows consistent deviations from 
the single average curve. To test for Reynolds number effects, 
the values of G which best fit each individual set of data were 
determined, and these were plotted, together with the diffuser 
results, as a function of Reynolds number in Fig. 8. The more 
reliable data do show the expected trend. It is also noteworthy 
that the results of the conical diffuser experiments agree well with 
the two-dimensional data. It should be remarked that the en- 
trance momentum-thickness Reynolds number is used as the 
independent variable. The point at the highest Reynolds num- 
ber was obtained from measurements made at only one station 
in the diffuser of the 48-inch water tunnel’ at The Pennsylvania 















State College. Also included are data obtained by Uram!* 
in a recent study of conical diffuser boundary layers 
5.6 
lo+6 
5.4F 
< % CONICAL DIFFUSERS (I7) 
5.2 > ® CONICAL DIFFUSERS(I9) 
sok + 48-INCH WATER TUNNEL 
48r 
4.6r 
44 
4.2r 
40r a 
386r 
iil 1 L 1 Roi 
L i L 5 
s0O0. (10>) 2xilo®) —s sxtl®—i‘ét* ss 2x10°* 5xi0* = 10° 25°98 


Fic. 3. Variation in exponent with Reynolds number (except as 


noted, symbols are defined in Table 1). i. 
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CONCLUSION 


A solution of the integral momentum equation has bx pr 
sented which is particularly simple for flows in adverse pressur 
The solution is a power relation, the exponent being 
Reynolds The 
method, both at low pressure gradients and near separation, hay; 


gradients 
function of the number limitations of thy 
been discussed and defined 

To solve a typical boundary layer problem, the step-by-ste 
solution given by Eqs. 4 and 5 should be used from transition to 
Calling this 
point the origin (subscript 7) Eq. 7 may then be used directly for 


the point where the inequality of Eq. & is satisfied 


any subsequent station. In using this equation, the value of ( 


may be taken from Fig. 3, or a mean value of G = 2.8 may be 
assumed if only a rough approximation is required 

For the cases for which the present method has been compared 
with experiment, the results are in each case at least as good as 
any other published method. The simplicity, when compared 
to the step-by-step solution of simultaneous differential equa 
tions required by most methods, should appeal to most workers 
in this field. 

The authors wish to thank Mr. E. M 
unpublished data on conical diffuser flow 


Uram for the use of his 


REFERENCES 


von Karman, Th., Uber laminare und turbulente Reibung, Zeitschrift fiir 


angewandte Mathematik und Mechanik, Vol. 1, pp. 235-252, 1921 
2 Buri, A., 


beschleunigter und 1 


Eine Berechnungsgrundlage fiir die turbulente Grenzschicht be 
Thesis 652, Federal Technical Col 
British Ministry of Aircraft Pro 


erzégerter Strémung 
lege, Zurich, 1931; R.T.P. Trans. 2073, 
duction 

Gruschwitz, E., Die turbulente Reibungsschicht in ebener Strémung bie 
Druckabfall und Druckanstieg, 


* Clauser 


Ingenieur-Archiv, Vol. 2, pp. 321-346, 1931 
Adverse Pressure 


Department of Aeronautics, Con 


Francis H., Turbulent Boundary Layers in 
Gradients, The Johns Hopkins University 
tract AF 33(038)9862, June, 1953 

> Newman, B. G., The Re-Attachment of a 
Behind a Spoiler 
Australia, October, 1949 

6 Wallis, R. A., 
Layer Flow, 
tories, Australia, November, 1949 


Turbulent Boundary Laye 


Report A64 of the Aeronautical Research Laboratories 


Turbulent Energy Considerations in Turbulent Boundar 
Aerodynamics Note 86 of the Aeronautical Research Labora 


Goldschmied, Fabio R., Skin Friction of Incompressible Turbulent Bound 


ary Layers Under Adverse Pressure Gradients, NACA TN No. 2431, August 
1951 
8 Bidwell, Jerold M., Application of the von Karman Momentum Theoren 


to Turbulent Boundary Layers, NACA TN No. 2571 
9 Van Le, Nguyen, The von Kdérmdn Integral Method as 


December, 1951 
Applied toa Turbu 
lent Boundary Layer, Readers’ Forum, Journal of the Aeronautical Sciences 
Vol. 19, No. 9, pp. 647-648, September, 1952 
Ross, Donald, Integration of the Reynolds Equations for Incom pressible 

Turbulent Boundary Layers, The Pennsylvania State College, Ordnance Re 
search Laboratory Report NOrd 7958-268, ONR Project NR 062-139, May 
1953 

1! Ross, Donald, Evaluation of the Momentum Integral Equation for Turbu 
lent Boundary Layers, Readers’ Forum, Journal of the Aeronautical Sciences 
Vol. 20, No. 7, p. 502, July, 1953 

12 Rubert, Kennedy F., and Persh, Jerome, Procedure for Calculating the 
Development of Turbulent Boundary Layers Under the Influence of Adverse 
Pressure Gradients, NACA TN No. 2478, September, 1951 

'S Ludwieg, H., and Tillmann, W 
spannung in turbulenten Reibungsschichten, Ingenieur-Archiv, Vol. 17, pp 
288-299, 1949; Trans. NACA TM No. 1285, May, 1950 

14 Hall, A. H., /nvestigation of Separation in Two-Dimensional Flow, Ph.D 
Thesis, Purdue University, June, 1946; Data published by R. C. Binder 
Journal of Applied Mechanics, Vol. 14, pp. 213-216, 1947 

‘% Wieghardt, K., and Tillmann, W., Zur turbulenten Reibungsschicht bet 
Druckanstieg, ZWB U and M 6617, 1944; Trans. NACA TM 
No. 1314, October, 1951 

6 Ross, Donald, and Robertson, J. M., Water Tunnel Diffuser Flow Stud 
ies, Part 3—Analytical Research, The Pennsylvania State College, Ordnance 
Research Laboratory External Report NOrd 7958-230, March, 1952 

” Robertson, J. M., and Ross, Donald, Effect of 
Diffuser Flow, American Society of Civil Eagineers 
No. 141, July, 1952 

18 Ross, Donald, Robertson, J. M., and Power, R. B., Hydrodynami 
Design of the 48-Inch Water Tunnel at The Pennsylvania State College, Trans 
actions of the Society of Naval Architects and Marine Engineers, Vol. 56 
pp. 5-29, 1948 


} 


Untersuchungen iiber die Wandschu 


November 


Entrance Conditions on 


Proceedings Separate 





for | 
entry 
ump! 
it wi 
sepal 
more 
than 
vant 
the ¢ 
the t 
is, ol 


ind 

mom 
the 

negle 
deter 
tary 
i firs 
Is CO 
two 
be I 
cedu 
rotat 
first 
one-' 
show 
inert 
lor € 
foun 
The 
seen 
spon 


in th 


pre 
ressure 
being a 
or the 
n, have 


by ste] 
ition to 
ing this 
ctly for 
1e of G 


nay be 


npared 
rood as 
npared 
equa- 
orkers 


of his 


rift fiir 


Laye 


tories 


indar 


abora 


sound 


ugust 


lorem 


ences 


ssible 
e Re 
May 


urbu 


nces 


h.D 
ider 


t bet 
™ 





READERS’ 


Ur Kk. M., Investigation of the Growth of an Axisymmetric Turbulent 
Layer in an Adverse Pressure Gradient, MS Thesis, The Pennsy] 
ania State College, Dept. of Engineering Mechanics, June, 1953 


Corrections and Additions to ‘‘A Modification 
of the Holzer Method for Computing 
Uncoupled Torsion and Bending Modes’’! 


Henry E. Fettis 

Mathematics Branch, Aeronautical Research Laboratory, 
Patterson Air Force Base, Ohio 

January 29, 1954 


Wright 


=< THE PUBLICATION of the subject paper,' experience has 
indicated that some revisions in the procedure originally 
recommended would render the method more useful. Also, two 


typographical errors have been noted which are hereby corrected: 


(1) Eq. 10(c) should read 


Mut (= a (**") My (= ? 
ee ek i EI, | \ 2 


2) The heading of column 14 of Table 2 should read: 
(2/Ax,) 


In the original paper it was suggested that the computations 
for bending frequencies be carried out in one table, with each 
entry consisting of two portions, one depending on the initial 
umplitude and the other on the initial slope. In actual practice 
it was found preferable to carry these two calculations out in 
separate tables. It was also found that in some problems it is 
more suitable to use the wing root as the starting point rather 
than the tip as originally suggested. This is particularly ad- 
vantageous when the tip mass is subject to variation, since then 
the conditions at the free end can be used to directly determine 
the tip mass as a function of the frequency. The disadvantage 
is, of course, that the symmetric and antisymmetric cases must 
be treated separately. For symmetric modes, the solution is 
started with zero slope and known shear at the root, while for 
intisymmetric notes, the solution begins with zero deflection 
ind a known moment. At the tip, the conditions of shear and 
moment are applied to determine the tip mass corresponding to 
the assumed frequency. If rotary inertia of the tip mass Is 
neglected, then the moment is set equal to zero and the mass 
If ro- 
tary inertia is included, then it is necessary to assume (at least as 
a first approximation ) that the radius of gyration of the tip mass 
The conditions of shear and moment then result in 


determined from the resulting value of the shear force. 


is constant 
two simultaneous equations from which the unknown mass may 
be found. When computations are carried out by this pro- 
cedure, an interesting result is obtained concerning the effect of 
rotary inertia which might not otherwise be recognized. It is 
first noted that if rotary inertia is neglected, the tip mass is a 
one-valued function of frequency, the relation being somewhat as 
However, when rotary 
That is, 


shown in Fig. 2 of the original paper 
inertia is included, the tip mass becomes two-valued. 
for each value of the frequency, two values of the tip mass can be 
found which satisfy the end conditions of shear and moment.* 
The result is a plot of the type shown in Fig. 1. It can be 
seen that values on the second branch of the curve, which corre- 
spond to the second mode may actually lie below other values 
in the first branch. In other words, for sufficiently large masses, 
the frequency of the second mode may be less than that of the 
first mode for smaller tip masses, a result which could not be ob- 
tained from the theory if rotary inertia were not considered. The 


* One of these values may sometimes be negative 
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Manner of variation of tip mass with frequency when 
rotary inertia is not neglected 


Fic. 1 


reason for this is that rotary inertia is more critical in higher 
modes as pointed out previously in a note in this JOURNAL.? 
Finally the following extensions to cases not considered in the 


original paper are noted: 
I BEAMS WITH DISCONTINUITIES AT AN INBOARD POINT 


To illustrate the treatment of this general type of problem by 
the present method, we consider the case of transverse vibrations 
and a cantilever beam, simply supported at some inboard point 
If the coordinates of the root, support, and free end are, re- 
spectively, x = 0, x = x,,x = L, the problem may be formulated 


as follows: 


For all x, 


(d2/dx?)EI(d2h/dx?) = wh 
h=0 
sities (dh/dx) = 0 
; EI(d*h/dx?) = 0 
saetieain (d/dx)EI(d*h/dx?) = 0 
Atx = 4M, h=0 


In addition to the above five conditions, we require that (dh/dx) 
and ./ be continuous at x = x However the shear force at this 
point will suffer a discontinuity corresponding to the reaction of 
Numerically, the problem is handled by consider- 
Starting at 


the support 
ing separately the cases 0 <x < x, andx,<x<L 
the point x = 0, and using an assumed value of w, the tabular 
method is employed to obtain two linearly independent solutions 
of the equation which satisfy the initial conditions of zero de- 
flection and slope. The process is carried up to the point x = x, 
at which point the condition of zero deflection is applied which 
determines the ratio of the initial values of moment and shear 
Similarly, starting at x = L and proceeding inward to the point 


x,;, two solutions satisfying the free end conditions are ob- 


x= 
tained. These are likewise combined to satisfy the condition 
of zero deflection at x = x,. The two portions of the solutions 


are thus uniquely determined with the exception of numerical 
factors, which now are chosen in the proper ratio necessary to 
Had the proper value of w been 
Thus by 


render the slope continuous 
used, the moments would then a priori be continuous 
plotting the moment discontinuity against the assumed values of 


frequency, the correct values of w may be located, 








Materials with Strain Hardening, J. Math. Phys., Vol. 27, pp. 
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Il. EXTENSION TO COUPLED MODES (d/dx)GJ(d6/dx) = —w?(Sh + 10 
The differential equations for coupled bending-torsion oscilla where S is the mass static moment per unit span about the toy 
tions are sion axis. The resulting algebraic equations for h,4; and @ 


in terms of previous values of deflection, slope, moment 
(d?/dx?)EI(d*h/dx?) = w*(uh + SO shear are* 


if) 


Ax\? Ax\‘ Ax\? 
hess { ot (**) So | toss | Glas } (=) Io | = in| -w# (*)'s,] 4 
; ‘ Ax \* Ax 
6, | Gr - «(*) 1 + 00! (SGI, + GJ 


which may be solved simultaneously for /,,4;, 0,4; in terms of 
BC On’: Dav Ma, one F, Advance values of 6, p, 7, and Fare 


then found as in the uncoupled case. Since three initial values : Fettis Henry E., A Modif Gro ay tNe Holzer Method of Computing Un- 
oupled Torsion and Bending Modes, Journal of the Aeronautical Sciences 
Volume 16, No. 10, October, 1949 


REFERENCES 


are not specified, it is necessary in this case to construct three 


inearly independent solutions for an arbitrarily chosen value . 
linearly i I u s for irbitrarily chos luc 2 Fettis, Henry E., Effect of Rotary Inertia on Higher Modes of Vibration 


of w. At the opposite end, application of the three remaining Journal of the Aeronautical Sciences, Volume 16, No. 7, July, 1949 


conditions will result in three homogeneous equations which will 


only be consistent if the correct value of w had been chosen. * Ax is assumed constant 


Errata on Alternative Forms for the Basic exponent in the numerator should have been instead of ? 
Equations of Transonic Flow Theory* so that the ordinate of Fig. 3 should read 

cy + IVa /t/e)” 

John R. Spreiter . “d ei 


National Advisory Committee for Aeronautics, Ames Aeronautical instead of (y + 1)%*cq,/(t/c) 
Laboratory, Moffett Field, Calif 


February 26, 1954 and of Fig. 4, 
(y + 1)M,?] /8cq,/(t/c) 
i HAVE DISCOVERED that the ordinates of Figs. 8 and 4 in my 


. — . n- rather than (y + 1)MMo?]/*ca,/(t/c 
recent note in the Readers’ Forum are mislabeled. The “ 
* Readers’ Forum, Journal of the Aeronautical Sciences, Vol. 21, No. 1 All the calculations and discussions are correct as are the curves 
pp. 70-72, January, 1954 themselves when relabeled. 





Thermal and Creep Effects in Work-Hardening Elastic-Plastic Solids 


(Continued from page 344) 


REFERENCES * Hill, R., Zhe Mathematical Theory of Plasticity; Clarendon 
Press, Oxford, 1950. 

3 Odgvist, Folke K. G., Influence of Primary Creep on Stresses 
im Structural Parts, Transactions, Royal Institute of Technology, 
10, 1948. Stockholm, Sweden, No. 66 (1953). 


' Hodge, P., and Prager, W., A Variational Principle for Plastic 
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